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PREFACE 


This report gives the results, in detail, for a stagnation-line 
analysis of the radiative heating of a phenolic-nylon ablator. The 
analysis includes flowfield ecupUng with the ablator surface, equlli- 
briu» ehemstry. a step-function diffusion model and a coupled line 
and continuum radiation calculation. This report serves as the 
documentation. 1. e. users manual and operating instructions for the 
computer programs listed in the report. Copies of the decks have been 
transferred to Mr. James K. Moss, grant monitor, of the Langley Research 
Center, and can be obtained from him or from the authors. 

This report also served as Carl D. Engel's dissertation require- 
ment in obtaining a Doctor of rhllosophy degree in chemical engineering. 
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Ja-nerq The Apollo 11 

u ^ iiiq qurroundrngs . i 

• Koq-! rally curious about his sun 
"" „as in a concrete »ay increased 

landing of man on the moon i curiosity 

•tarlnde the solar system. This basi 

tnan's surroundings t planetary 

.. • ^ -Itcplf to assure that tnci. 

: .,„rtfic o.:ectives for missions to t„e pianets 

• .. The main scienti 

missions. tne ma „f the Lunar and 

11 icif.R chairman oi mu 

. 1 1 ct-pied by Findlay, 

has been concise a x. C 1 1 ) "To learn more of the 

Hoard of NASA, (Ref- 

nautetai-j - explore the possi- 

^ svste^n ana lo 

and hiatorP of tUe 

--- :: :;:: ianetar. expiorntion nidi reRoire 

.. P,..ieni reeiitntion .e 

vinnlncrical advances, ihe g 

„any techno „ technological area 

. ..to contribute to the advancement of 
herein is t nlanetary mission. 

1 1 nd for a successful manned planeta y 

vhich must be eve op- successful 

t inctiral developments required lor a 
Ph.„y of the technological 

. oore discusaed and documented at 

manned planetary mission - - 

Technology for lianned Planetary Missions held 

. 1 2) discussed the limitations 

---^^“^""^Tilc lt^and the undeveioped potent 

-- -^-'\:::;::::s from the perspee 

° onts must he made in structural design and 

tequlrements. Improvements mus 
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synthesis (Rot. 1.3). Developments In hardnerc, software, m.an- 
machine integration .and other sub .areas of gnidance and control 
technology .are recinlred (Ref. 1.4). The reliability and length of 
operation of present life support systems most be increased (Ref. 1.5). 
Advanced con™„ieations hardware most be developed to accomodate 
high data rates transmitted from interplanetary distances (Ref. 1.6). 
The technologies for man-machine integration and experimental design 
among others must be developed. Furthermore, planetary entry and 
landing technology must be improved (Ref. 1.7). The vehicle 
heating aspect of planetary entry is the technological area of 

concern of the present research. 

The technology which must be developed for a successful manned 

earth entry from interplanetary tretve. c... e-eue-e... , --- -.i 

interp lane tear y trajectory, atmospheric braking trajectory and 
vehicle shape used. These trajectories determine the type of thermal 
environment, which in turn, determines the type and amount of heat 
shielding required to adequately protect the entry vehicles. 

Interplanetary trajectories are primarily determined by minimum 
energy and transit time constraints. For a manned mission to Kars 
two types of Interpl.inetary trajectories are under consideration. 
These trajectories consist of either a direct flight to Mars or a 
Venus swing-by as Illustrated in Fig. 1.1 from Ref. 1.8. The swing- 
by mission has two major advantages. First, additional scientific 
information could be gained by a manned pass near Venus. Secondly, 
the entry velocity upon earth arrival is significantly reduced due 
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, . with tho Ctavltauional field o£ Venus. This is 
to tho encounter with the gr 

illustrated in Fis- " 

to 1„, period. Barth entr. speeds are displayed for hoth direct 

end Venus s»ins-hy traiec tories . The left ed.e of each har rn 

Pi, 1.. indicates velocities for traiectcries optimised assu.rnP 

l^rs atmospheric hrahin« for capture while the ri.ht ed.e indicates 

velocities for traieetories optimised assuming propulsive ratmo. 

n Itinr entry velocities for several planetary missions are 

The resulting euLi-> 

presented in Fig. 1.3. It is noted from Fig. 1.3 that the cart 
entry velocity is hetween 11 and 19 Wsec for essentially a 

missions considered, 

Phe earth arrival velocity from interplanetary traiectorres 

eevr-t- entry trajectories. 

provides initial conditions for eai ... - 

These traieetories are limited hy heating, aerodynamic and docc r 

constraints. Due to hum.sn limitations a loading limrt of nomrna y 

10 g's (i.e. 10 times the earth gravitational force) is place on 

For a direct entry this trajectory 
the possible entry traieetories. For 

/n C 1 q^ The aerodynamics ot 
13 called the undershoot boundary (Ref. 1.9). 

the entry vehicle, determined primarily by body shape and weight, 

eatahlishes the sice of tho usable entry corridor. Furthermore, the 

hody Shape has a significant effect on the surface heating. Thus an 

h-cle shape is a compromise between maximising manuverabali 
optimum vehicle shape 

and minimising tho heat input to the vehicle. 

The magnitude of the hinetlc energy which is dissipated during 

atmospheric brahing is proportional to the velocity s,uared . Since 
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vexocu. .0... U. .3 aXso voU.«3.U no„ 

per unit area of streaxntuba toward the vehicle, it follows that 

. ■! P flux, toward the vehicle is 

energy flow per unrt area, i.e. tiux, 

,„ro«Xo„aX 30 X.XcaX e„x. V.03..3 

X. a. X.3 a. aoxX. .. X. .OX. 

T , n Im/sec it is evident that the amount 
velocity is approxxmately 

energy £Xng Is snbstantlaXXy higher £or gXanetary 

of energy and energy 

•lies Such velocities, which arc greater than 
return velocrtrcs. 

, ,ra caXX. hyperboXic since intcrpXanctary 

earth escape veXocrty. a 

rraxectcry £or hinetic energies o£ this .agnit.e result in hyper- 

bolic shaped trajectories. The main concern, of co , 

aeter.ine the iraction o£ this hinetic energy .hich is transXerrea 

rn the uehicle's sur£ace during at.espheric deceXeratron. 

-ihilitics of experimentally simulating the 
The present capabilit 

- anticipated during hypetbeXlc entry is iXXustratc 
flight conditions antic ipc-cu 

u x-iicr present facilities are not 
in Fig. 1.^* This figure shows tha P 

1 txlv producing both high enthalpy and flow 

capable of simultaneously produci g 

Ki nf building a test facility which piovides 
energies. The problem of buildi g 

required energy £Xuk and enthalpy is a a,aJor one (Ref. X.7). 

.XternateXy, £ree-£Ught .odeX er.perX.ents eonXd he performed, but 

1 riiff-tcult to scale. Hence, 
-.-.raTipc; arc very expensive, and diffi 
such experiments arc veiy r 

our hnowlcdgc would not be rapidly Increased even if the bag, cos 
„ero accepted (Re£. 1.10). These are the reasons why there has een 
and will continne to be a considerable reliance placed upon 
analytical techniques to predict atisosphorlc entry heating and the 

resulting surface material response. 
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THE VI8C0US hyperso:jtc shock t.ayer prof.lkm 

In this work the interaction of the stagnation region flovj- 
field and an ablative protection system of vehicles which experience 
a hyperbolic atmospheric encounter is analysed. The analysis is 
aimed at resolving major uncertainties in the current state of 
knowledge by systematically examining the effects of radiation, 
viscous coupling and ablator coupling. These effects are measured 
in terms of the heating rate to the ablator surface. 

The processes wliich govern the heat transfer rate to a blunt 


vehicle in a hypersonic flov7 are, for the most part, contained in 
3 layer adjacent to the vehicle. This shock layer formed by a 
blunt body in a hypersonic stream is shovm in Fig. 1.5. The 
radiating snock laycj- (white region) is the thinnest in the nose 
region where the highest heating rates are experienced. For 
protection from the high heating rates encountered during entry an 
ablator will be used. An ablator is a surface material which absorbs 
heating loads by changing phase ai\d mass loss thus reducing the 
transmission by conduction to the interior of the protected material. 
In terms of weight efficiency for entry deceleration an ablator 
protection system requires 10 to 50 times less entry vehicle weight 
than would be required by either nuclear or chemical propulsive 
system (Ref. 1.11). Vehicle weight minimization is quite iinportant 
since one pound of re-entry spacecraft weight requires approximately 
300 to 1000 pounds of launcli vehicle vjeight (Ref. 1.11). Many types 
of ablative materials are available which could be used for vehicle 


Fig. 1.5 (Next page this location) 

Photograph of a blunt body in a multiple arc 
jet wind tunnel showing the ionized shock 
layer about the body. (Mach no. prf 7) 
Courtesy of: T. A. Barr, Jr., U. S. Army 

Missile Command, Redstone Arsenal, Alabama, 

1969. 


Fig. 1.6 (Next page this location) 
Photograph of a cross section of a phenolic 
nylon ablator. Courtesy of: C. W. Stroud, 

NASA TN D-4817, 1968. 
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protection. Typical of the better neiterlals is the charring phenolic- 
nylon ablator (Ref. 1.12). Fig. 1.6 prc.sonts a photograph of a cross- 
section of a phenolic-nylon ablative composite ohich has been 
exposed to a severe heating environment and nhich decomposed to a 
char of porous carbon and loo molecular weight gases. This type of 
ablator protects the vehicle not only by subliming the char, but 
also by acting as a heat sink, by transpiring high energy gases 
which reduce convective heating, by reradiating from the char 
surface, by reacting exothermically and by blocking radiation from 
the shock layer. These ablator processes are intimately coupled 
with the shock layer processes. Thus to determine the amount and 

type of ablator to use for entry protection a quantitative under- 

, .so T-ro QTion is nGcdcci. 

standing oi both tne stiocr. • 

The major mode of energy transfer to the surface for entry 
velocities above 11 km/soc Is by radiation as illustrated in Figure 
1.7. This figure shows that, for typical flight conditions of 
interest, the convective heat transfer coefficient rs essentially 
Independent of velocity whereas the radiative coefficient incieases 
rapidly dwarfing the convective coefficient for the higher velocities 
The principle reason for the high radiative transfer is the elevated 
temperatures experienced in the shock layer (i.e.. on the order of 

15000*^K) for typical flight conditions. 

Figure 1.8 presents a schematic of the important regions 
ablation coupled shock layer. The outer region of the shock layer 
primarily consists of high temperature radiating air which was heated 
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P«cceain. .hoc. wave, So»o of fhc va.iaUon fs o.atcea 

,.„cu. awa, foo. Che ho., an., ff noc ahsovho. fs loaf U.vou.h 

t,.o shoe. wave. The va.iatloo loss re.ucos Che ho^pooatooe so the 

Shoe. ha.et an. thus te.oces the ta.tatWe flow. This .esctlptioo 

excludes the consideration of radiation absorptron by the 

v.irrher the radiation 

strean, gas called precursor radratton. Further. 

emitted fro™ the hot air region is partiall. ahsorhed h. the coo or 

air and ablation species and the remainder is absorbed by the ablator 

r-1 r enrfoce is reradiated into the 

o ^ fh-ic; energy to the suric^ce xb 
surface. Some ot tnis 

1 bv the ablator surface, 

shock layer and the reminder is abso 

(■ o 90 '’^ of the free- stream flow rate) 

For large ablation rates (i. • '« . . • j 

u ,1 b.Pll as the near shock region is inviscid. 

the region near the body as well as the n 

The convecrive heating rate is much smaller .h.n indicated in Fig. 

•pp exists between these two inviscid regions m 
A viscous region exis 

which sovor.nl Important processes occur. In tho vis-ous rog 
energy is transferred hy viscous dissipation and gnsoous conductron 

in addition to radiation. Moreover, species concentrations change 

A hhe laree temperature change in 
quite rapidly due to diffusion and too large P 

T7*rr 1 R Is fot a casc of large 
the viscous region. The schematic in Fig. . 

however as the ablation rate is lowered the viscous 
ablation rates; however as 

tegion becomes attached to the surface. In this situation convective 
beating to the surface hecomes a more significant heat transfer mode 
to the ablator. Furthermore, the effects of air species diffusing 
to the surface and a consequential inhom.ogenuous surface reaction 
become probable. In addition to these processes, 

03 a particle flows away from tho a.ts of symenetry (i.e. the 
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reduced and 

• it is accelerated, its temperature r. 

stagnation ) surface 

decrease.. Consequently the heating rate 
the pressure dccrea.eo. 

. line e.porianccs tl« most severe 
The reclon near the stasnation Ime .p 

.fieri region in a deslsn analysts. 

A thus is the inost critical o 

eonserv.rtle„ la„s nhieh .eserihe the ...nation 

:: ;:.io: are more Simpie than those nhioh 

.fa fore their solution provides a succrnct an 
shoch layer. Therefore, their 

f the entire shock layer for a 

cprvative characterization of the 

conservative „„^phers have mathema- 

f o' flisht conditions. Numerous research 
given set oi fl. 

rlcally modeled the near stagnation region p 

degrees of completeness during the past ten to fifteen years. 

. _ c 1.7, 

kt'CGUciy three icvi<-*'^ 

lushed uhich diseuss the profuse amonnt Of norh 

f t-hc most current v;ork v-hic 
done. Consequently only some of 

u tT.ll be reviewed. 

pertinent to this research «U 

CO r 1 15 1.16. and 1.17) present 

Several recent papers (Ref. . . 

f the blunt body flow Problem whleb include line an 
solutions to the 

eontinuum air radiation coupling and assume ^ 

et al (Ref. present stagnation poi 

inviscid. rage et.al. layer 

radiative and convective heating ■ 

, „e The effects of ablation products 

pressures and thicknesses. The 

n, f d (Pef 1 16) calculated radiative heating 
included. Olstad (Ret. 

• r an inverse method (i.e. specifying the s )OC 
distributions using an 

, h. this method describes non- 
cor^-yrntinc the body shape), this men 
shape and computing im 

. r. only The shock layer was assume 

stagnation line regions onl> . 
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consist o£ an inviscld air and Inviscld ablation layer, altbongh air 
propcrtica were uacd in tbc ablation layer. Tbe radiation fodel 
used included line and continuum mechanisms and was coupled 
flow field solotion. Chin (Ret. 1.16) asaerr,ed the stagnation line 
shock layer flow could be divided into a inviscld air and an inviscid 
nblation layer. Computationally each of these layers were divided 
into six sublayers. The coupled radiative heating rates cere 
eompwted with a model which contained line and continuum mechanisms 
of air and ablation producis. A carbon-pbenolic ablator was coupled 
to the flow-field analysis by assuming equilibrium sublimation. 

In order to determine the shock layer processes more precisely, 
many researchers (Ref. 1.18 to 1.24) have Included viscous effects 
in their £iov;-iield models. ..a-iSw... v 

line solutions of fully coupled viscous, radiating shock layers 
including specified ablation injection. The momentum equation was 
aolved by finite differences using two methods. One method is valid 
for small mass injection rales whereas the second Is valid for large 
rates. The effect of radiation blocking by carbon-phenolic ablation 
products was studied. Shock layer elemental compositions wore 
determined using a single species equation for ablation products 
diffusing into air. The radiation model used contained line and 

1 -omc fnr C 11 0 and N atoms. A limited 

continuum radiative mechanisms for C, H, 

number of solutions were presented. Rigdon et. al. (Ret. l.W) 
presented several viscous, radiation coupled stagnation line solutions 
which included cases tor specified carbon phenolic ablation rates. 
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The viscous stnEHUtlon Hue mo^uu™ equation vas solved 

T„ltlal value technique uhleh „as slatted at the stasnatlon point 

then carried out toward the body and toward the 
The integration v;as then carriea 

shoch- The radiative model used included Une and co 
a«chanls.,s o£ air and ablator species. A species equation 

solved uslns an effective binary diffusion coefficient for 

1 rfci The effects of precursor radiation 
of air into ablation products. 

1 rRof 1 20) reported solutions for a 
.•.nrirl Smith et. al. (Ret. 

were examined. omiiu 

quasi-steady response of a phenolic-nylon ablator to a flow- re 

solution which was broben into two regions. The outer air legion 

f'i-T'Tn tnethod . 

was analysed as an inviscid region using a one o 

The inner ablation layer was analysed using two techniques. o 

^ used; v?heio.ac :or 

small ablaciou rates a boundary .aye-., 

•niPfiril method was used. No radiation 
large ablation ratos an integral t.e 

coupling between air and ablation Inyer was oonsidere . 

^ael used included various Une and contlnuur. reohanls.s for ato„s , 
ions and molecules excluding line weehanism for C and 11 atosts. 

Heating rates for tbe stagnation point and around tbe body along a ^ 
trajectory were presented. Engel and Spradley (Ref. 1.21) prdsen e 
stagnation point radiativo boating rates for a typical hyperboll 

a radiation model which contained only 

entry trajectories using 

uontinnum mechanisms of air. The viscous momentum equation was 
solved using ibe integral technique of UoshfcaRl et. al. (Ref. 1. 
which IS limited to low ablation rates. Tbe ablator response and 
resulting estimated beat shield weights were computed using the calcu- 
lated cold wall heating rates m an uncoupled ma 
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other vlscoes hnelyscs »hlch endohvor to define the 

effects of finite rate chemstry and culticcwonent diffusion In 
the shoeh layer, have not inclnded radiation (Ref. 1.23 to 1.26). 

Adaos et. al. (Ref. 1.23) present results for chemical „onet,uiUbrtu. 
fuviseid and laminar viscous floo over spherically blunted cone 
Oeo^trics. The calculations were »ade for flight velocities and 
altitudes where radiative transport is negligible. The chc.istry 
was restricted to air species and injected species of argon, heliu. 
or carbon dioxide. Multicoraponent diffusion was discussed; however 
„„ results were presented, for the cases studies the influence of a 
noncatalytic wall with .ass injection significantly reduced the 
convective heat transfer. Pavy et. al. (Ref. 1.24) presented 
solutions of the iiiul ticompouent , reacting, stagna.-c. p— ..t 
layer with choc, leal eguiUhriu..,. A comparison with binary solutions 
was given. A nitrogen or air external stream with injection of 
hydrogen or a miKture of hydrogen and argon were the chemical systems 
studied. The results presented indicate that the and H species 
compositions computed with multicomponent diffusion were the only 
species that significantly deviated from their binary solution. Mu 
(Ref. 1.2A) studied the finite rate chemistry effects m injecting 
hydrogen into air at an axisymmetric stagnation point. The hydrogen- 

1 -to™ M-iq restricted to low stagnation temperatures 
air chemistry system was rescricuc 

where ionisation does not occur. Mu demonstrated the difficulties 
associated with the chemistry of these flows hut did not propose a 
general solution. Blottncr (Ref. 1.26) investigated a finlte- 
difforcnco method and a nonlinear ovcrrolaxation method 
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the stagnation line viscous blunt body problon.. A seven species air 
rodel nas used to study nonequiUbtium chemstry eifccts. Results 
„ere presented for flieht conditions ..'here radiative transfer is not 
inrortant. Numerical studies indicated that the finite difference 
procedure converged nore rapidly than the nonlinear overrelaaation 

method . 

The studies reviewed represent the best currently available, 
unfortunately no single analysis includes all of the important 
effects of radiation and ablator coupling although many of the analyses 
adequately account for certain shcch l.nyer processes. The inviscld 
analyses lacR the generality of being applicable to small ablation 
rates. «itb the exception of Chin (Ref. 1 . 16 ) the inviscid analyses 
have not adequately accounted for ablatio.-, product affects rn the 
shoch layer. Unfortunately Chin's results do not include all of 
the important line radiation mechanisms. The viscous analyses, 
which are by and large the best, ate limited by either nusmrlcal 
difficulties or computation time in addition to Incomplete radiation, 
flow-field or ablation models. The main limitations of Wilson's 
(Ref. 1 . 18 ) analysis are the lack of molecul.ar radiation and 
numerical difficulties with the momentum equation. The analysis of 
Rlgdon et. al. (Ret. l.W) Includes the most detailed radiation 
model which unfortunately consumes excessive computation tine and 
this analysis has not been used to study .ablator coupling. The other 
viscous analyses reviewed do not include as many Important ablator- 
shock layer processes as the two just discussed. The third group of 
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papers reviewed did not include radiation transport effects which is 
the main heating mode being studied. These papers are helpful in 
deciding how the particular processes examined might be handled 
numerically and in understanding the role that finite rate chemistry 
or multicomponent diffusion might play in the shock layer. 

In summary, the studies reviewed are inadequate in varying 
degrees. Further, the analyses vjhich are complete enough to define 
some shock layer processes have not been used successfully to define 
many cause and effect relationships. Consequently many questions 
remain unresolved. The most important of these questions are; 

1. Under what conditions, if any, can the shock layer be 
treated as inviscid? 

Is nolccu’' ir radiation important and if so what molecules 
contribute? 


3. How effective is ablation in reducii'i:^ the heating rate to 
the surface? 

4 What arc the ablation rates v.-hich correspond to 
estimated shock layer heating for various flight condition; 

5 Vlhat is the error introduced into analyses by assumptions 
made in the flow- field model? For examiile, shock shapes. 

6. Can the shock layer heating process be correlated in a 
simple manner vjith any shock layer parameters? 

In addition to the uncertainties remaining in current heating 
analyses, there is a significant uncertainty in the basic data used 
in these analyses. Several major areas v;hcre basic data is not 


precisely known are: 

1. Gaseous radiation data 

2. Surface cmissivities 
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3. Sublimation temperatures 

4. High temperature transport properties 

Using the 1968 state of tlie art as a basis, the estimated heat 
shield weight would account for 12 to 24 percent of the total 
vehicle weight for a 50,000 ft/sec entry velocity trajectory 
(Ref. 1.27). The 12 percent range of uncertainty in the total 
weight is due directly to uncertainties associated with defining 
the flow- field and ablator behavior for various vehicle shapes. A 
detailed study of uncertainties due to incomplete analyses and a 
discussion of the effects of uncertainties in available data were 
made in an attempt to reduce them. This document describes this 
study. 

PRE8]::NT KhShAKCH e)hJhCTIVi:S 

The overall objective of this research is to develop tlie 
capability to accurately predict the performance of ablative thermal 
protection systems when exposed to aero- thermal environments such 
as those encountered by planetary atmospheric probes and return 
vehicles from interplanetary missions. Emphasis is placed on 
diminishing some of the uncertainties presently existing in ablative 
thermal protection design. Specific research objectives are: 

1. Investigate the interaction of the stagnation region and 
around the body flow- fie Id and ablation protection system. 

(a) Develop the governing equations from a general 
property balance and systematically point out 
the assumptions made in obtaining the equations 
to be solved. 
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(b) Describe inathematically the various levels of 
approximation which may be introduced into the 
equations describing the shock layer flov;. 

2. Calculate the stagnation flow- field with tlie use of a 
radiation transport model vdiich does not impair the 
accuracy of the best available radiation property data. 
Emphasis is placed on making the resulting program an 
engineering tool . 

(a) Develop a numerical method to solve the momentum 
equation v^hich is stable and equally valid for 
largo and small ablation rates. 

(b) Determine a procedure for evaluating the relative 
contribution of molecular species radiation and include 
the ones necessary for an accurate x^adiation calculation. 

(c) Examine numerical simplicat ions for the species 
equation. 

3. Perform parametric studies on entry velocities, altitudes 
and vehicle shapes, and then determine if a more simple 
empirical model could be used to represent the r*esults 

oI lliose detailed calculation^.-, 

(a) Examine the effects of ablation at various rates 
to determine cause and effect relations, 

(b) Obtain coupled ablator- shock layer solutions for 
various flight conditions. 

(c) Exc^minc parametrically the effects of stagnation 
line shock bluntness on surface Ideating rates. 

(d) Determine the importance of molecular species 
radiation in effecting the surface heating rate. 

In general, the research reported in this dissertation is a 
study of the uncer taint ics which exists in entry heating analyses. 
Furthermore, it is a definitive study of the thermal radiation 
which occurs during such entries. 
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CHAPTER 2 


DEVELOP^^::NT OF GOVERNING EQUATIONS 

As with most physical problems, the flow about a blunt body 
entering the earth's atmosphere obt-ys the conservation laws of mass, 
momentum and energy. The equations representing these laws for a 
multicomponent, radiating, chemically reacting fluid in an inertial 
reference system are derived using a general property balance approach 
in Appendix A. Following the derivation, the conservation equations 
are written in time independent vector form from w’hich they arc v?ritten 
in general orthogonal coordinates. In the third section of Appendix A 
the conservation equations are written in orthogonal body oriented 
coordinates for apnlication to the blunt body flow problem. In this 
chapter these equations are simplified using physical arguments and 
order of magnitude assessments. Care is tahen to indicate the 
approximations made throughout the development. 

THIN VISCOUS SHOCU LAYER EQUATIONS 

In order to determine the proper mathematical model to describe 
the flow-field developed by a blunt body moving at hypersonic 
velocities, one must assess the behavior of the gas that the vehicle 
will encounter. Fig. 2.1 based on the work of Ref. 2.1 presents 
the flight regimes which are encountered by a body during atmos- 
pheric entry. The regimes can be grouped into two gasdynamic 
domains - continuum and noncontinuum. Hayes and Probstein, Ref. 2.2, 
demonstrate the continuum domain can be divided into five regiiiies: 
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(:i^) ufcj-i 


CM 


Kegimes (,Bas(jd on Ref. 2.1 an> 
Modified Fron Ref. 2.27) 
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( 1 ) classical boundary layer, (2) vorticity interaction, (3) fully 
viscous, (4) incipient merged layer, and (5) fully merged layer. 

The beliavior of the gas flowing over a body in tlie five continuum 
regimes can be described using the equations developed in Appendix A. 
Let us consider further tlic characteristics of fluid flow in the 
five continuum regimes, 

1. Boundary layer regime: The classical boundary layer 

equations are a valid approximation of the viscous effects 
for high Reynolds numbers corresponding to lower altitudes. 
Viscous effects dominate near the wall in a region which 
is small compared to the shock layer thickness. Vorticity 
generated by shock curvature is therefore negligible thus 
not affecting the boundary layer flov;. 

2. Vorticity interaction becomes important at lower Reynolds 
numbers where shock generated vorticity becomes significant 
in respect to viscous effcct5j near tl'ie body. Here the outer 
region of tlic shock layer, usually considered the inviscid 
layer, becomes coupled through momentum transport to tlie 
higher shear region near the body, usually thought of as 

the boundary layer. The high shear region near the body is 
also larger than that experienced at higher Reynolds numbers. 

3. Viscous layer Regime: Viscous effects from the body inter- 

action are spread throughout the shock layer (i,e., the 
boundary layer and shock layer thicknesses are of the same 
order). This occurs at lower Reynolds numbers and 
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correspondingly higlier altitudes than docs the vorticity 
interaction regime. Viscous dissipation at the shock is 
still small m comparison to dissipation at the body. This 
condition is true so long as tiie ratio of the mean free 
path behind the shock over the shock layer thickness is 
much smaller than the square root of the density ratio 
across the shock wave. Ref. 2.2. This implies that the 
Rankine-Hugoniot shock wave equations are valid for deter- 
mination of shoe!; layer boundary conditions. 

4. Incipient merged layer regime: The incipient merged layer 

begins when dissipative effects at the shock are significant 
The shock wave is thin relative to the shock layer thickness 
but the Rarkinc-l!ugcniot relations must be modified to 
account for viscous effects at the shock boundary. 

5. Fully merged layer regime: At higher altitudes and low 

Reynolds numbers a distinct sliock does not exist. The free 
stream mean free path over the major body radius is approxi- 
mately one or less. The flow behaves continuously from the 
free stream to the body. Above this altitude range continuum 
concepts are no longer applicable and the flow goes through 

a transition to free molecular flow. 


The foregoing discussion of the five continuum flow regimes 
follows in part the reasoning of Hayes and Probstein, Ref, 2.2. This 
reasoning was based upon tlie assumption that radiative energy 


transport and ablative mass 


injection were negligible. 


In the present 
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development these two effects arc the primary flow field-body inter- 
action mechanisms which are to be assessed. Fig. 2.1 shows the flight 
regimes where radiative heating to a one foot splierical body becomes 
significant. For the most part, significant ablation rates are also 
encountered in these regimes when using present day charring ablators 
such as carbon phenolic or phenolic nylon. Therefore, let us make 
additional observations about the flow characteristics in these 
continuum flight regimes where the effects of ablation and radiative 
transfc). in the shock layer arc important. In proceeding, 
our attention will be restricted to the first three flight regimes, 
where the heating rates to a vehicle's surface are the most significant. 

Significant radiative energy transfer has several important 
effects on the shock layer behavior. First, radiative transfer 
couples the energy equation and thus the thermal boundary layer over 
the entire shock layer. Tliis is apparent by recalling that the flux 
divergence term in the energy equation is evaluated by an integration 
over all space in the shock layer. This effect has been demonstrated 
by several auUiors including Ref. 2.3 and 2.4. Further, the thermal 
boundary layer exists from the shock to the body for all three fliglit 
regimes in the radiative coupled domain. Secondly, radiative energy 
transfer produces nonadiabatic or energy loss effects. Principally, 
radiant energy is lost through the transparent shock wave. Thirdly, 
the effect of radiative transfer in the shock wave is coupled through 
the energy equation to the momentum equation. Although this coupling 
effect 1 .S not altogetlier negligible, it does not change the conclusions 
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obtained about momentuin transfer in the shock layer in the first 
three flight regimes. Therefore, even though the viscous effects 
may be approximated through boundary layer concepts with possible 
modifications of edge conditions in the vorticity layer regime, the 
energy transport occurs over the entire shock layer. In the viscous 
layer regime both viscous and energy transport arc significant over 
the entire shock layer. 

Appreciable mass injection rates of ablation products results 
in additional effects on energy and momentum transfer within the 
shock layer. High mass addition rates tends to enlarge the region of 
shear dominated flov? near the body. Libby, Ref. 2.5 showed experi- 
mentally and theoretically that in the boundary layer regime, 
boundar^^ layer concepTs could bo applied \rhcn mass injection or 
suction rates were quite large. This study did not include the 
effects of radiation, but since energy transport does not change 
the character of momentum transport these conclusions are also valid 
insofar as momentum transfer is concerned for radiative coupled 
shock layers. Mass injection has other effects such as reduction of 
shear at the wall. Ref. 2.4, and reduction of heat transfer at the 
wall. Ref. 2.4, 2.5 and many others. These effects although of great 
importance do not change the basic characteristics of momentum or 
energy transfer in the shock layer. 

We may conclude that for flight conditions in the radiative 
coupled domain where ablation rates are also significant, tlie character 
of the momentum transfer is essentially the same as without these 
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effects. „o„evet, tl.c Catactctistic, „t enet.y ttensfet arc 
stEntftcantly ,t«ctc„t t„ t„at the c.ttte choc. la,cr .c 

considered in all three fli(th^ 

foregoins stateraents as background tlie problem ohicli 
we Wish to solve can be stated. The basic conservation equations 

stated in .ppendia A are appropriate to describe the flo„ oi a 
continuum reactinp 

trng and radiating gas mixture over a blunted suriace 

when thermodynamic equilibrium e.vists For the 

tor the present work, we will 

determine the reduced set oi equations which desetibe the How i„ , 
sHock layet ovet a blunt body when the outer boundary ei the shock 
layer rs a shock wave described by the Ra„ki„e-„„g„ni„t equations 
Thus the equations governing the flow in the shock layer will be 

applicable to tbe three higher keynoids number regimes both in a„d 

out of the radiation coupled domain Th^ 

cicmain. The prime concern and motivation 

-T Obtaining this set oi equations is to describe the heat transfer 

condrtrons can be predicted by numerical calculation. 

^ll£^.o lMa.onatude An.nltrom. 

Tn order to determine tbe appropriate sot of equations which 

realistically approximate the flow sitnat' 

situation just described, an order 

magnitude assessment of the terms in the basic co 

Dasic conservative 

equations is needed. This iq i 

. . . “ ‘==■"100 out by first uondimen- 

sionali..ing the conservation equations. The folli ■ 

inc following nondimens ional 


variTbles are introduced. 
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It should be noted that the equations in the Appendix A are in 
dimensional form. In this chapter a superscript * will denote 
dimensional variables unless it is explicitly stated otherwise. 

The dimensional global continuity equation is: 


.y . 0 
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( 2 . 2 ) 
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Using the dimensionless variables stated in Eq. 2.1 the above equation 
may be written as 
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Dividing by p U yields the dimensionless form 
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(2.3) 


From Eq , A. 7 8 the d iri'.ens ional species continuity equation is: 
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Units of the terms in Eq. 2.5 are mass of specie i per mass of 
mixture . 

Substituting Eqs. A. 86, A. 88 and A. 89 into Eq, A. 84 yields 


the dimensional x-momentum equation: 
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Proceeding as before the din’iOnsionless variables are introduced. 
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and cancelling dimensional terms yields the dimensionless 
^-momentum equation. 
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The y-momentum equation can be stated in dimensional form from 
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Introduction of noad ir.icMis ional variables into the y- momentum 
equations follows the same procedure and pattern as in the x-momentum 
equation. The resulting nondimensional y-momentum equation is: 
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Using Eqs. A. 90 with A. 86, A. 87 and A. 89 the energy equation 
can be written in dimensional form: 
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A 

wIiGiTc the diffusions! snd rsdistivc flux divergence terms j > nre 
defined in Eq. 2.1. Substitution of the nondiinensional ratios from 
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Introducing the Re^^nolds niunber the uondiinens ionnL energy 
equation can be v;ritten: 
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Having stated the nond inensional conservation equations we are 
confronted with the problem of estimating the relative magnitude of 
the terms in each equation. 
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According to the results of Hayes and Probstein, Ref . 2.2, the 
gas behind a bow shock of a hypervelocity vehicle is a continuum 
for freestream Reynolds numbers greater than 100 based on principle 
body radius. Further, the standoff distance nondimensionalized by 
body radius for flight Reynolds numbers greater than 100 has been 
shown, Ref. 2.2 and others, to be approximated ly equal to the 
density ratio across the bow shock. The density ratio for hyper- 
sonic Mach numbers is of the order of one tenth and less for 
dissociating gases. These stated relationships can be expressed 


as follows: 
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Since »e nrc concerned cith a thin layer »ith respect to the body 
radius, Prandtl's concepts for the relative order o£ Basnitude of 
terms in the conservation equations can be enployed. Following the 
procedure given by Schl ichtlng, Ref. 2.6, the relationships for the 
relative order of nondinensionaliaod terms may be written. 
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Using the above estimates the relative order of magnitude of 
the terms in the four conservation equations have been determined. 
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Energy: 
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Using the lower limit on Reynolds number v;e observe 


J_ -2 1 

R « P Ri 2^00 
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(2.19) 
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At this lower liniit on Reynolds number, let us drop all terms of 
o 

order p and higher in all equations except the y-momcntum equation. 

—2 

In the y-momentum equation terms of order p are retained for a 
specific reason. Along the stagnation line, ^ = 0, the u component 
of velocity is zero. Thus the y~momentum equation is of one order 

lov;er at g = 0. It is appropriate in this case to consider terms of 

— -2 

two orders of magnitude in this equation namely p and p . The 
resulting conservation equations are presented in Tab, 2.1, 

At § = 0 the y- morucntum equation, Eq. 2.23, has terms v^hich are 
of order p and p^“. Two terms which can be directly eliminated from 
this equation when u ~ 0 at C = 0 are indicated by arrov;s. 

The simplified set of conservation equations, Eqs, 2.15, 2.20 - 
2.24 form a set of partial differential equations (neglecting the 
radiative terms) v;hich are valid for Reynolds numbers greater than 
100. It is obvious that the terms v/hich have been dropped due to 
order of magnitude reasoning become less significant as the Reynolds 
number is increased. Tiiese *’thin shock layer" equations are the same 
as second order boundary layer equations with curvature terms and 
arc valid for continuum flov; of the viscous, vorcity, and classical 
boundary layer regime. 

To this point little has been said about the bulk viscosity v/hich 
appears in the \ term in the momentum and energy equations. This has 
been done for the sake of generality. Hov/ever , to interpret the 
pressure in our equations as the local thermodynamic pressure Stokes* 
postulate 



TABLE 2 . 1 

listing of conservation equation with 

ORDER ASSESSMENT RESUT.TS 
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Global continuity; 


0 [ 1 ] 0 [ 1 ] 

||(Au) + = 0 


Species continuity: 

0 [ 1 ] 0 [ 1 ] 
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0[1] 
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A 5u 1 




(2.15) 


( 2 . 20 ) 


( 2 . 21 ) 
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TABLE 2.1 (Cont.) 
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momentum: ( 0 [p] and larger terms) 
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(2.23) 


y. 
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TABLE 2.1 (Cont.) 
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momentum: ( 0 [p] and larger terms) 
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TABLE 2.1 (Cent.) 
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2^1* + 3X* = 0 (2.25) 

must be accepted. The bulk viscosity "ja is a direct indication of the 
departure of the ir.can pressure from the thermodynamic pressure 
expressed by the equation of state Ref. 2.7* Further, Laitone Ref. 
2.7 points out that by accepting Stokes postulate for compressible 
flows v;e are at best restricted to monatomic gases. This appears to 
be a rather stringent assumption until one examines the type of 
behavior a polyatomic gas must exhibit to significantly deviate from 
monatomic behavior. To a first approximation the bulk viscosity 
"Ji characterizes the dependence of pressure on the rate of change of 
density Ref. 2.8. Gases which exhibit sliov7ly excited internal 
dc£ 7 :ees of freedom (i.e. rotational or vibrational) in flows which 
have rapid changes in the state of the fluid, the pressure cannot 
follow the changes in density and differs from its thermodynamic 
equilibrium value. Thus, acceptance of Stokes* postulate fox“ bulk 
viscosit}^ is consistent v’ith our basic assumption of local thermo- 
dynamic equilibrium used throughout this development. Henceforth, 
we will use 


•k 2 

X = - 3 p, 


(2.26) 


in our equations. In thin shock layer equations Stokes' relation 
is needed only for the y-momentum equation. The order analysis has 
eliminated all terms containing X in botli the x- momentum and energy 
equation. 
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We note that dimensional equations will be used throughout the 
rest of this section without the superscript notation unless the 
superscript is needed for clarity. 
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Let us consider liie radiative transfer Eq. A. 28. 


1 

7 7~ + Qt • VI = O' (B - I ) 
c ot 1 V V V 


Follovjing the steady state assumption let 


1 


SI 

V 


c Bt 


0 


Therefore our transfer equation can be written 


Q-i 


Vl 


-ff (I - B ) 

V V V 


(2.27) 


By imposing tlic one-diniensional approximation, the radiative transfer 
equation for the y-dircction may be v/ritten 


(( j/ 




dl 


dy ) dy 


— (I 

V V 


- L ) 


(2.28) 


where j h cos Y 

For tlie one-dimens ional problem the absorption and emission charac- 
teristics vary only in one direction^ y. This fact is sufficient 
inforrution to solve Eq, 2.28 for the specific intensity by integration 
in y. We will see later that although the specific intensity is 
evaluated one-dimensionally the radiative flux and flux divergence 
must be evaluated over all space. Consequently the flux divergence 
is integrated over an infinite plane slab v/hich has the same intensity 
variation across the slab at any station down the slab. 

In order to clarify the solution of Eq. 2.28, Fig. 2.2 is 
presented. From Fig. 2.2 we observe 
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dy = cos Y ds = fds 


(2.29) 


By defining the optical depth as 


- J “v 


(2.30) 


and using Eq. 2.30 the radiative transfer equation can be 
rewritten 


dl 

^ = I - B 
dj V V 


(2.31) 


The radiative transfer Eq, 2,31 can be solved formally by using the 
variable coefficient method: 


Substitution of the above relation into Eq. 2.31 and solving for 
the function C(t^) yields 






(2.32) 


Thus the general expression for the specific intensity is 

Tv 2 

'^v,l 

Splitting the integration into tv-’o parts and evaluating boundary 
conditions yields 


I - I + I 

V V V 


where 
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\ exp (- (t^- t^//) 



T 


V ,W 


I 

V 


+ I ^(t ) exp (- (t - t )//) 

vv,w ‘ v,w V 


+ [ B exp (- 

J V 

T 

V 





dT 

~1 


(2.34a) 


(2.34b) 


+ I 

V 


(t ) exp (- 

V, E 


v.s 


\)//) 


The above equations describe the radiation field in terms of 
temperature througli Planck's function for a nonscattering gas. 
The quantities I "(t ) and I ”^(t ) are boundary conditions and 

^ V V , s V V , w 

the exponentials represent attenuation over optical path length. 

Using Eqs. 2.34 for the specific intensity, the radiative flux 
and flux divergence may be evaluated. Recalling from Appendix A, 
Eq. A. 22, radiative flux term can be expressed as 

„4rr 


qR<r) 




^0 


I n, dQ dv 
V 1 


(2.35) 


For the geometry under consideration the unit vector can be 
replaced by the direction cosine f . From Fig. 2.3 we note that 


and 


dQ = sin Y dY d0 


j = cos Y 
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Therefore 

dQ = - d/ d0 

Substitution of Eq. 2.36 into 2.35 yields 


iR,y 


,oa pi 

I d0/ df dv 

0-1^0 


(2.36) 


(2.37) 


si»i>iifyi"e £»y £1'=' »"= *>5' i”££8£‘>££“" £" « 

r“ 

q = _ 2rr f / 

^R,y ^0 -1 

It is convenient to split the integration in Eq. 2.38. 


(2.38) 


% = 2n I l//d/ 

^R»y>v Jq 

o_ = 2 tt ! I /df 

Thus the monocrorriiitic heat flux is the 


(2.39) 


^R,y,v ^R,y,v ‘^R.y.v 
Substituting Eq. 2.34 into 2.39 yields 


,.T 


qZ = - 2n 

%,y,v 


’^VjW 


R.y.v 






(2.40a) 


(2.40b) 
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where the direction cosine, f, dependence is expressed in terms of the 
exponential integral function of order n. 


E 

n 



(-t//)df 


(2.41) 


Let us examine the radiative flux equation given in the third 
section of Appendix A. 


q 


R,y 


|.y(rp 

\(^o) 



Differentiating v?ith respect to 


I (r)dnjdv dy 

J 0 V ' 

y we obtain 


(2.42) 


9q 




ay 




v;h i c h is 1 1 ic r a J i a t i v e f 1 u 
energy equation, Eq, 2.24, 


R>y Sy 


’R.y' 



divergence in tlie y direction. 
v;e have the term 


(2.43) 


In our 


Due to the one-dinensi onal planar slab approxina tion this term will be 
represented by 


5 ^ A 




yj 


A ^'•r. 




+ q 






As a result of this approximation, an evaluation of Eq. 2.43 is 
sufficient to describe the radiative transfer influence in the energy 
equation. 

In order to evaluate Eq. 2.43, the intensity at a fixed point y 
and in a direction define:! by 1 and f is integrated over all solid 
angles. Substituting for the solid angle, the integration for a 
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one- d iincns i ona 1. plnnc slab can be readily carried out. 


Sq 


9y 


a 


0 


0"'’ 


where the inner integral is 

.-1 .0 . 


“ -T ), 


v,w 




r J 3^ __^a/ 

0 T 


(2.44) 


(2.45) 


r” ^ ' "\r ^ '\t'\ 

* J„ I— If 


(t„ -t). 


Eq. 2.45 can be simplified by interchanging the order of integration 
as substituting the exponential integral function. 


p-1 

1, V^ = J 


B E,(^ - T )dr 

V 1 V V V 

, w 


-I (t )E„(t -t ) 
V v,w 2 v,w V 


+ E (t -t )dT 

J V 1 V V 

''v 

-I (t )E„(t -t ) 
V v,s 2 V v,s 


(2.46) 


Substituting Eq. 2.46 into 2.44 provides an expression for the 
radiative flux divergence in a one-dimensional slab. 
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-5q 
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Sy 


2tt 




vL J 


v,w 




(2.47) 


r V,: 

J 'v^r'v 'v' 


+ J • E, (t -T„)dT 
T, 


+ V‘"v.s>'=2<\.s-'v> - '“vJ"" 

where the exponential Integral function I!„ has the following 

characteristics: 


E (t) = E (-t) for n = 1, 3, 5, 7 , ... 
n n 


£ (t) = - 


(2.48) 


Eq. 2.47 is valid for arbitrary boundary conditions with the cxceptron 
that only one boundary reflection of a photon packet rs allowed. In 
practice, for a shock layer solution, the subscript "w" is interpreted 
as conditions at the wall or body and ”s" as conditions at the shock. 

Under this interpretation = 0 barring precursor radiation 

. ._L .,11.- — n Further, for the case 

and the optical depth at the wall 0. 

of a perfect absorbing wall 1*^(0) = 0. These simplifications are 
the usual ones made in describing radiation transport in a radiating 
shock layer. Raking these simplifications reduces Eq. 2.47 to 

Eq. B. 31 of Ref. 2.10. 
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In order to numerically compute the flux and flux divergence 
when considering line radiation it is advantageous to solve the 
radiative transfer equation in physical space rather than in 
optical depth space. To this end the following development is given. 

The radiative transfer equation is 


- ; 


dl 

V 

dy 






(2.31) 


The formal solution of Eq. (2.31) follows the same procedure stated 
before with tlie exception that f is assumed to be an average value 
of 1/2 (Ref. 2.12) 


I + = + 

V 




'o 


/ 


dy 


(2.49) 


O' B 

r v V e 


I = + . 

V Jy / 


1| a- 

-a a dy 

/Jy V 


dy 


These two equations are the counter parts of Eqs. 2.34 a and b with 
the boundary conditions set to zero. From Eq. 2.39 the radiative 
flux equations inay be written 


qn'' = 4n ^ 

jo V 

(2.50) 

V " Jq 

Finally using Eq. 2.44 the flux divergence equation may be v.^ritten as 


By 



a (2B 

V V 


I 

V 


)dv 


(2.51) 
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The one-dimensional planar slab approximations which result in 
Eq. 2,51 have important ramifications to our shock layer problem. 
Radiation calculations can be made using Eq. 2.51 at each § 
location independent of other ^ locations. This makes the thin 
shock layer equations a set of parabolic partial integro-dif ferential 
equations which can be solved using marching schemes which are used 
for solving classical boundary layer equations. 

An observation concerning the planar slab approximation is in 
order at this point. This approximation eliminates all curvature 
effects from the radiation calculation. A more appropriate approxi- 
mation for m,ost axisyn^^etrically blunted vehicles would be a concentric 
sphere approximation for the boundaries of the shock layer as pro- 
posed by Viskant.., Ref. 2.13. For a two- dimensional body the 
corresponding approximation is quite obviously concentric cylinder 
boundaries. However, as pointed out by Viskanta, Ref. 2.13, 
comparatively little attention has been given to radiative transfer 
in curvilinear systems. The paper by Viskanta analyzed the steady 
state radiative transfer betv.een two concentric, gray, opaque sphere., 
separated by a gray absorbing and emitting medium which generated 
heat uniformily. He concluded, for constant absorption coefficients, 
that curvature effects were evident for concentric sphere radii 
ratios as high as .99. This corresponds approximately to a shock 
standoff distance of 6/R .01. Nominal hypersonic standoff distances 

are .OA 6/R ^ .10- Viskanta' s vjork we arc led to expect that 
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curvature effects may be significant for both radiative flux and flux 
divergence in a typical shock layer. The actual magnitude of these 
effects aie difficult to assess because of the constant absorption 
coefficient assumption and differences in boundary conditions for 
the problem under consideration. Thus as far as is known an accurate 
assessment of curvature effects on shock layer radiative transport 
is absent today. It is felt that using a concentric splicre model is 
analogus to including both first and second order effects whereas 
the infinite parallel plate model includes only first order effects. 
However, for the present we will use the infinite parallel plate 
model in our development. 


Statement of Shock L a e r E ci u a 1 1 c j p. s 


i.s a resuj 


a s EUU'P 1 . ion , and tlie planar radiative triinsfer model the second order 
thin shock layer equations may be v;ritten in a more usable form. In 
addition, a relation for the y-component of the heat flux vector from 


Appendix A was used to yield the second order equations given in 
Tab, 2.2. 


Let us now examine the simplifications v;hich are needed to 
obtain the first order shock layer equations and classical boundary 
layer equations from the equations stated above. First let us drop 
all terms of order p or smaller. The resulting first order shock 
layer equations are given in Tab. 2.3. Additional simplifications 
can be made by assuming the boundary layer thickness to be small in 
comparison to the local body radius. This implies 
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TABLE 2.2 

SECOND ORDER SHOCK LAYER EQUATIONS 


Global continuity: 


(r%.) + I; Cr%v) = 0 


Species continuity: 

i (r%C,u) + I (~r\c^v) = 


Sy 


+ Ki: O)^ 


X - Momentum: 
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, S A f^Uv 5r g. 

+ ~ \K'^ P 
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y - Momentum: (0[pl larger terms) 


A Sv , ^ A 5v ^ 

pr u — + pnr V — - p^r u - - nr 
+ 3^ " 3 by cx ^ + 3 by ^ by^ 


2 b . A , ^ br' 

- 3 („r + «*v 


■) 


(2.52) 


(2.53) 


(2.54) 


(2.55) 
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TAIiLE 2.2 (Cent.) 
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(2.56) 


(2.57) 
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^ ^ A A 

K ® > K "♦ 1 ^ 


(2.61) 


Using these limits in the first order shock layer equations results in 
the boundary layer equations given in Tab. 2.4. 

Equations 2.62 through 2.66 arc essentially the same as the 
boundary layer (B.I..) equations v;hich are given by Fay and Riddell 

Ref. 2.14, Dorance Ref. 2.15 and others. The boundary layer equations 

A . 

for a flat plate are obtained by simply noting that r^ is not a 
function of x. We can conclude from the foregoing simplifications 
of the thin shock layer equations that the classical Prandtl type 
boundary equations contain only first order terras which exhibit no 
normal component curvatux^e effects. 

Boundary Condil:ions 

Having stated the thin shoeV; layer and boundary layer equations, 
the appropriate boundary conditions for the two sets of equations can 
now be discussed. Figure 2.4 presents a sketch of the various 
regions and boundaries of particular inteiest in the t.iin shock 
layer problem. We note that, in addition to the shock layer region 
the char layer and decomposition zone (sec Fig. 1.6) are important in 
our problem. These regions are important because the momentum, 
energy and mass transfer in the char and decomposition regions are 
intimately coupled to the transfer in the shock layer. Theoretically 
we could consider all the processes which take place between the shock 
wave and the virgin plastic of the body and attempt to solve the 
governing equations for this boundary value problem. However, it is 
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TABI.K 2.4 

BOUI'IUARY layer equations 


Global contimiity: 




(2.62) 


Species continuity; 
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(2.63) 


(2.64) 


(2.65) 


( 2 . 66 ) 
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iDore practical to divide the solution of this general problem into 
a shock layer and a material response problem and iterate on the 
boundary conditions at the material surface. Therefore, it is 
important to realise what information is available from the material 
response solution which could be used for boundary conditions of 
the thin shock layer equations. This is accomplislied by integrating 
the conservation equations across the chai‘-gas interface and reducing 
the spatial increment to zero to yield surface balance equations. 

With tills perspective of the general problem in mind, the nature of 
the thin shock layer equations and boundary conditions V7ill be 
discussed . 

The parabolic nature of the thin shock layer equations mathe- 
matically requires initial conditions ar well as boundary conditions 
in order to obtain a solution. The entry veliicles axis of symmetry 
is the appropriate location of the starting line for zero angle of 
attack problems. The determination of conditions along this line, 
called the stagnation line, is a major and important problem in 
itself. Consequently, development of the method to obtain these 
initial conditions (i.e. stagnation line solutions) is delayed until 
after the boundary conditions are established. 

As discussed in the next section, the thin shock equations are 
a set of parabolic integro-differential equations with initial values 
given along x = 0, the stagnation line. Because the shock Xv’ave 
location is not known before hand, the blunt body problem is mathe- 
matically referred to as a free boundary problem. Given initial 
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conditions along the stagnation line and boundary conditions along 
the body, the thin shock layer equations can theoretically be solved 
with a simultaneous development of the shock geometry and corresponding 
shock boundary conditions. The shock geometry (sec Fig. A. 3) can be 
obtained by carrying out the following integration. 


pX * * 

= (1 + K 6 ) tan edx + 6^ 


(2.67) 


In practice another technique has been used to determine the 
shock geometry Kef. 2.3, 2.10 and others. The shock geometry is 
assumed and specified in terms of de/dx. Iterations are made around 
the body until the input and output shock geometry coincide. 

If the shock geometry is known, the Rankine-Hugoniot equations 
can be used to obtain the shock boundary conditions. The development 
of these equations in curvilinear coordinates follows directly from 
Ref. 2.10. The dimensional Rankinc-Hugonoit equations vjritten in 
rectangular coordinates are: 


Continuity: 


* * * ^ 
0 V = P V ^ 

^cD cOjll SjH 


( 2 . 68 ) 


Momentum: 
(norma I) 

(tangential) 




(2.69) 

(2.70) 
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Energy: 


% V ^ + hM ^ 
^ 00, t <» 


‘2 + h * = V + % V + h * (2.71) 

» ^ s,t ^ s,n s 


Using Fig. 2.5 the above equations can be written in body oriented 
coordinates. From geometry we have 


* * * 

V = V sin e - V cos e 

s s , t s , n 


(2.72) 


* 


* 


u = V ^ cos e - V sin e 

s s , t s , n 


(2.73) 


where 


V = U cos ci) 
o>, n ^ ' 


* — * 

V = 0 u cos rn 

s ,n ^ fo 


CO 


V _ = V t- = U sin cp 

S , t CO, t CO 


Substituting for V V “ and ^ Eqs. 2.72 and 2.73 yield 


V * = U sin CD sin c - p U cos cp cos e 

S CO ^ 


(2.74) 


* 


u = U sin CD cos p. - p U cos cp 

r' rr, ' CO 


Sin e 


(2.75) 


The pressure behind the shock can be obtained by using the normal 

ic 

momentum equation and substituting for and • 

(2.76) 


* * 2 * * — * . 2 * 
p (U cos co) + P ~ Po ^s 

•CO m ' CO S ^ ^ 


By 


substituting normal and tangential velocities the energy e.quation 


can be v;ritten 
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h 


* 


s 



CO 


2 



+ V 


*2 


) + h 


(2.77) 


The tangential velocity is negligible near the stagnation line and 
thus reducing 2.77 to 


* 


*2 


(1 - 


P 


cos 




+ h 


(2.78) 


Nondimensionalizing Eqs. 2.74 through 2.78 and dropping and h^ 
which are order (p ^) yields the following shock boundary conditions 


V == sin CD sin c - p cos cp cos c 
s ^ 


u sin CD cos G + p cos cp sin e 
s ^ 


P = (1 •• p) cos cp 

s 


-2 2 

h = (1 - p ) cos cp 
s 


(2.79) 

(2.80) 
(2.81) 
(2.82) 


or 


h 

s 


2 2 

1 - (u ^ + V ^) 

s s 


(2.83) 


It is important to realize that the Pvankine-llugoniot relations are 
valid only if strong precursor radiation effects do not become 
important. The shock conditions can be more adequately described 
for the strong precursor radiation problem with modified Rankine- 
Hugoniot relations presented by Zeldovich and Raezer (Ref. 2.16). 
However, significant precursor radiation effects are not experienced 
in air below flight velocities of approximately 60,000 to 65,000 
ft. /sec. as demonstrated by Lasher and Wilson (Ref. 2.17). Therefore 
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the Rankine-Rugonlot relations provide satisfactory boundary condi- 
tions for the outer edge of the thin shock layer equations for many 
problems of current interest in atmospheric entry. Let us now 
write the shock boundary conditions at y = 6. 


u = u 

s 


V = V 

s 


p == p 

^ ( 2 . 84 ) 


h = h 

s 


C = C (P h ) (Assuming chemical equilibrium) 
i is s s 


The Rankine-Hugoniot equations provide expressions for u^, v^ , P^ , 

and h . The equation of state and free stream elemental mass fraction 
s 

provides the additional information needed to determine the post 
shock species mass fractions assuming chemical equilibrium. The 
specific intensity coming through the shock towards the body is 
specified as zero. We note that in total four boundary conditions 
are needed for the energy equation because of its integro- 
diffcrential nature. Thus two boundary conditions, enthalpy and 
specific intensity, have been specified at the shock. 

The corresponding body surface boundary conditions can be 


written for y - 0: 
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u = 

0 

V = 

V 


w 

p = 

p 


w 

h = 

h 


w 

C. ^ 

= c. 

1 

1 


(t 

I 

1 ^ 

V, 

1 

T 


V, 

V? 


The boundary conditions specified in 2.84 and 2.85 are sufficient to 
solve the thin shock layer equations. However, substitution of 
equivalent boundary conditions for some surface conditions is found 
to be practical. For example the normal velocity at the wall is 
usually replaced by (pv)^. Of greater practical importance is the 
wall boundary condition on pressure. This pressure is not known a 
priori . An equivalent boundary condition is then needed. There are 
at least two suitable boundary conditions which might be used in 
lieu of pressure. These are the normal pressure gradient at the 
shock or the normal pressure gradie t at the body. The normal 
pressure gradient at the shock could be specified by evaluating the 
inviscid y - momentum equation at the shock using the R.ankine- 
Hugoniot equations. The normal pressure gradient at the body could 
be set equal to zero from boundary layer theory. Each of these 
conditions would involve some degree of approximation. 


The effect 
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of these approximations are not evaluated here. To evaluate the 
pressure gradient at the shock an approximate form of the continuity 
equation is needed. Correspondingly the zero normal pressure 
gradient assumption at the wall neglects the wall velocity head 
at the body vdiich would push the true stagnation pressure point off 
the body- An additional complicating factor arises when one 
observes what boundary condition is needed in the material response 
analysis- The pressure at the outer wall is usually specified as a 
boundary condition Ref. 2.18. Ideally one would like to knov; and 
specify the pressure boundary condition for both problems. This 
would eliminate iterating on this variable between the two solutions. 

The surface boundar^^ conditions can be derived by integrating 
the conservation equations across the boundary and taking the limit 
as the spatial increment approaches zero. This method assures 
inclusion of all the effects accounted for in the flow- field 
equations • 

A photograph of a section of charring ablator is shown in 
Fig. 1,5 in vdiich the important zones are indicated. The ablative 
corr.pos ite * s response during entry may be analyzed in two v;ays . One 
is a transient analysis which gives the response of the material as 
a time function of its heating environment. The other is a quasi- 
steady analysis which predicts a constant, history independent, rate 
of decomposition for a given heating environment. Experimental 
evidence indicates the conditions under which the quasi-steady 
behavior exist. As the material is lieated, the surface is removed 
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by chemical reactions, sublimation and erosion; as a result, the 
total thickness of the material decreases as shown by the data in 
Fig. 2.6. Concurrently decomposition in depth and a char build up 
occur. When an equilibrium situation exists such that a constant 
char thickness and a constant surface recession velocity are 
maintained, a quasi-steady response would physically exist. 

Additional evidence that a quasi-steady ablator response would occur 
during hyperbolic entry has been presented by Ref. 2.29. Theoretical 
calculations using a transient ablator analysis and a quasi-steady 
analysis were shown to yield essentially identical results over tne 
peak heating portion of a typical trajectory. This portion of the 
trajectory is the conditions of current interest. 

As indicated in Fig. 2.6 the char depth is of the order 0.3 
inches for lunar entry conditions. This thickness should be nearly 
independent of heating rate and thus applicable to hyperbolic entry 
conditions. The flow through the char can be considered one dimensional 
unless the ratio of the char thickness to the local surface curvature is 
somewhat smaller than .05 (Ref. 2.19). For ratios of the order .05 
or greater the flow-field pressure variations may cause multidimensional 
flow through the char. Since the surface radii under consideration 
are of the order of 1 to 15 feet, the one dimensional flow approxi- 
mation should be quite good. Moreover, since the porous char and 
decomposition zones are quite thin with respect to the body radius, 
the pressure through the zones can be assumed constant. 

The quasi-steady surface balance equations were derived by Esch 
In Ref. (2.20) using the preceeding assumptions applied to equations 



Total Thickness 





Time (sec) 


2.6 


Thickness of 
as a Function 


the Char and Virgin 
of Time (Ref. 2.28) 


Plastic 
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2.52 to 2.55 and 2.57. From the development in Ref. (2.20), the 
spGciGs boundciTy conditions sr6. 

Species boundary conditions: 


(pvc.)" + Jf + (Ri + spci - €p) = (pvc.) + J. 


( 2 . 86 ) 


where = mass rate of formation of species i by heterogeneous 

reactions 

s = mass rate of sublimation of species i by homogeneous 
i 

reactions 

e = char porosity (volume of voids per unit volume) 

P 

and where the superscript - and + means evaluated on the char side 

and flow- field side of the surface respectively. 

^ ^ ‘ e he obtained by multipli." 

The elemental bounoaiy ca.. dc uul... . j 


cation of Eq. 2.86 by 


e . . 
i-J 


A. .m. 
1-1 I 

M. 


(2.87) 


where 

A == moles of element j per mole of compound i 

ij 

m = atomic weight of element j 

j 

M = molecular weight of compound i 

i 


and summing 

over 

all compounds i. 



n 

n 

n 

(pv) 

1 

e . .C . + ) e J . 

ij 1 Zj ij 1 



i=i 

ii 

i 



n 

n 



(pv)-^ I + 

Y e.. V 

-L. IJ 1 



i=l 

i=l 


( 2 . 88 ) 
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The above equation can be expressed in terms or 
fractions, C., and the elemental mass fluxes. J.. 
Elemental species boundary condition: 


cd in terms of the elemental mass 


(pvC.r +X" = (pvC.) +Jj 

The momentum boundary condition was obtained using the y - 


( 2 . 89 ) 


■'omenturo equation of order p . 

Momentum boundary condition: 


P+ = p- - (pv 


4 / dvx"^ it 

“3 r dyj 

The surface energy balance equation uas derived noting that no 
significant radiative transfer occurs ulthin the char. 

Energy boundary condition: 

(l+a)V f 1' - pvl'',(':i''-'=i') + >'->1 

1 j^i " ^ 

dT ! ” , + 

+ (1 + k)\ dy 1 

„ 1 C the sum of the radiation flux to the 
The radiation term, , 3.S tne s 

surface from the flow- field minus the reradiated energy flux. 
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These surface boundary conditions express the coupling relation- 
ships which exist between the flow- field and the ablator. If 
diffusion is significant near the surface, the species balance is 
a boundary condition of the third kind. The moipentuin balance 
indicates a complicated coupling of the surface pressure to both the 
internal and external normal velocity profile if terms of order p 
are retained. The energy balance is a boundary condition of the 
second kind in temperature with additional convective, mass and 
thermal diffusion, species and radiation coupling. The complexity 
of these boundary conditions suggest simplifications must be made, 
where possible, to arrive at a tractable set to be used. A simplified 
set of boundary conditions is selected in the next section. 

Typical boundary conditions for the boundary layer equations 
can now be discussed in terms of the ones used for the shock layer 
equations. Outer boundary conditions along a line between the 
shock and the body knov;n as the boundary layer edge are used rather 
than the Rankine-Hugoniot equations. These edge conditions are 
usually obtained using some inviscid layer analysis which is 
bounded by a shock and a streamline. The method of characteristics 
is used for the supersonic portion of the flow and typically a 
Belostserkovskii strip integral technique is used for the near 
stagnation subsonic flow (Ref. 2.2). These methods provide the 
following boundary layer edge conditions 
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u - u. 


V = 


p == p 

e 

h - or g = gg = ^ 

^ C. (assuming che 

j ) (usually not used) 

V ^ v,e 


(2.92) 


;hemical equilibrium) 


The boundary layer wa 
can be written. 


11 boundary conditions that are u 


sually employed 


u = u. 


= (ov).. 


p = p = p 
^ w e 


(2.93) 


h = h^ or g = g„ 


Ci = 


+f ■) = B (usually not used) 
\ V 


If the spectral intensity is eliminated from the previ 

co„.uio „3 ..e, e.u.vax.n. to t.oso ptetentea tn 

, of Ref 2 15. One can observe that the proble,. of 
Chapter lofRel. 

response solution Is eU„lnate,a. Houever. this prohle. Is le t 

unresolved in that the correct edge pressure can be ^ ^ 

aecurately onl, through an Iteration procelure between the rnvrser 
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flow analysis and the boundary layer analysis. It is also signi- 
ficant to point out that, although usually not attempted, it is 
computationally rather difficult to handle boundary layer and 
inviscld flows which ore coupled by radiative transfer. In addition 
to the geometrical integration problems the boundary condition on 
apecific intensity or radiative fltw. is not a single value but a 
frequency dependent function which must be matched at the boundary 


layer edge. 

Stagnation Line Equations 

To this point we have not discussed how initial values for the 
thin shock layer equations may be determined. This problem is of 
near equal importance to the entire shock layer problem and will 
be discussed in the remainder of th^s sect.^cn. T- obt n 
values for the shock layer solution, a reduced set of the thin 
shock layer equations must be solved at x = 0 along y. the stagnation 
line (see Fig. 2.4). The solution of this of equations is of major 
importance because (1) the highest heating rates and pressures on a 
body are experienced at the stagnation point (2) any distributional 
shock layer solution because of its parabolic nature is only as 
valid as its initial values and (3) the thin shock layer equations 
along characteristics x = constant reduce to ordinary differential 
equations like at the stagnation line. Thus by developing a 
stagnation line solution an important problem is solved and a great 
deal of the work is completed which is applicable to the total 
shock layer problem. This is primarily why the stagnation line 
problem has received a great deal of attention in the past decade. 
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The solution to the stagnation line (S L) problem by direct 
methods has been approached in two vays. The »orh o£ Ho and Probstein 
(Ee£. 2.22) typifies the stagnation region solutrons which use cj.pan 
sions of the dependent variables in x to obtain the stagnation and 
near stagnation line equations. The work of Hoshiraki and Wilson 
(Ref. 2.3) typifies the stagnation line solutions which determine 
the stagnation line equations by formally taking the limit of the 
terms in the shock layer equations at x - 0 using symmetry conditions. 
The latter method is used in this development. 

From this point on in the development, attention will be 
restricted to axisymmetric bodies for which the exponent A • 1 (i.e. 

. r). With this restriction noted, let us first examine the 
global continuity equation in expanded dinenslonal form. 


V or' X Qu I ~ .SX = 0 

l^(pu) + + r r Sy 

As X -♦ 0 the following limit is approached 

I5r rx j 1 - u 5r _ ^ 

lim 7 ® r Sx 9x 

x-»0 ^ 

assuming a spherically shaped body at x = 0 

1 ^ _ r sin Q I ^ — 2i = ii 

r By L(1 /k + y) sin 8J I + kV - 


(2.9A) 


(2.95) 


Also, note that 


(2.96) 


Using these conditions the global continuity equation can be 


rewr itten. 
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Global continuity (S L ) 


2 — (pu) + ~ Gtpv) + Hpv - 0 
Sx oy 


(2.97) 


The species continuity equation can be rewritten by subtracting the 
global continuity Eq. 2.52 from the left hand side of Eq. 2.53. 


A '^i A-. ^ 

r pu — + r Kpv 


M V A 

Qr 


3 A 


(2.98) 


3y 3y 

Noting that at . 0. n = 0 and using Eq. 2.96 in Eq. 2.98 yields 
Species continuity: (S L) 


Kpv 


SCi 


dy Sy 




(2.99) 


or by noting 


ey 


5C_. . 2. 


pv 


sy 


i = _ (j ) - -^ j. + uj. 

By '' i.y-* ~ i,y 1- 


(2.100) 


Now cons 


ider the x - mornentum Eq. 2.54 


■£ Ip + Ht puv - ^ ^ 



iy 


By evaluating the above equation at x = 0, relatively little infor 
mation is obtained. Along the stagnation line u = 0 for all y; 
therefore 



0 


( 2 . 101 ) 
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Using this inforiTVition in Eq. 2.54 yields 

/ ^ \ ^ Q (2.102) 

4=0 " 

which agrees identically with the Rankine-Hugoniot equations for a 
synmetrical shock (i.e. cp = 0 at x = 0) . The reduction of Eq. 2.54 
to 2.102 along the stagnation line yields the expected physical 
interpretation that no momentun is transferred in the x-direction at 
the stagnation line. Since this trivial form of the momentum 
equation is not useful, the rate of change of momentum in the x - 
direction is used. Therefore let us differentiate the x - momentum 
equation with respect to x and determine its limiting form along 
the stagnation line* 


A B^u 
r pu — 2 
Sx 



(r^pu) 




+ 


Sx 




5u 

9y 


5 u 


9u 




2 

- " bx bx by Ux^'" ^ oy ^ SxbyV 

3x 


(2.103) 


A. 

Sx 





(r^p) 


After some manipulation and substitution for limit quantities Eq. 


2. 103 reduces to 
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(2.104) 


^ i!5.o 


-k[pv + 1 ^^ + ^|^](sx)“ P " ^2 

H 

Por substitution into E„ 2.104 the stngnetion line global continuity 
Eq. 2.97 may be re'.'jritten. 


or 


2 ^ = - r -i |- (^pv) + xV ] 
9 x L p Sy 

2 ^ _ r - 4- (pv) + 2 kV 1 


(2.105) 


(2.106) 


ax " L p 3y 

Combining Eqs. 2.104 and 2.105 yields 
X - Momentum: (S L) 

% ( “ ¥ ‘ Tp ty + [»" - ^ ( k 


^ ii + j- ^ ]f i- Cpv) + f-) 

+ U ^ IV 9 n ' 2 7 


+ - Kj_pv + P t ^ By JV 2p oy 

H 


(2.107) 


+ p( i-p f ^ ^ U 1=0 ■ ° 

This is a third order Inhctnogcnous ordinary differential equation 
where the rate of change of the pressure gradient in the x - direction 

is an undetermined function of y. 

The y - momentum equation can be evaluated directly by substi- 
tution of the stagnation line limit quantities. The stagnation 

1 “2 . 

normal momentum equation to order p is 
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- Momentum: (S L) 


C)V 

pKV 


Sy 


Sy 


■55 . ^ I ^ (pv) + HVJ 

Sv ^ Sy L 2p Sy J 


^ ^ f ^ (pv) + ^IkV 

3 Sy L 2p Sy 


K 


(2.108) 


4 

3 


^v_,4~S /^N 

+ Sp 3 Sy Sy ^ M-y 3 


rs^ 

K 


{[i 


p. ^ + ^u^ 
K 


p,v 

2 


(pv) 


2 ^ \ 
3^’- Sy J 


—2 

where the terms in the brackets { ) are the terms of order p . By 
dropping these terms only terms of order p remaxn. Since some of 
the terms of order 1 have been expanded in Eq. 2.103 a fcv of the 
terms will combine, 
y - Momentum: (S L) 


~ Sv 


~ .^ ( pv) + kV^ 

H -N„ P" Sy \ 2p Sy ' 


Sy 


(2.109) 


+ II; (t; I? ‘p''’ 2p ly 

Jr ^ 

S , 4 ~ s / i* „ ^ 

+ kV j + 2 hP Qy + 3 5y Sy^ " 3 ^ sy 

It is obvious that either with or without the second order terms the 
y - momentum equation is a second order, inhomogenous . ordinary 
differential equation with variable coefficients. Given a solution 
to the energy equation (i.e. an enthalpy or temperature profile) in 
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principle the x - and y - momentum equations could be solved for the 
normal velocity and the normal pressure gradient if the rate of 
change of the pressure gradient in the x - direction as a function 
of y is given. This pressure term, (S^P/Sx^) . is usually assumed 
since it is a result of the elliptic nature of the problem. 

The energy Eq. 2.57 can be reduced to the stagnation line 

energy equation by inspection. 


xpv-^ = - Cl +^) 


3y 


K 


Sy L 


c^T 


! - Kh ^ + K / h,J, 


V 

L “i“i.y 

i 


( 2 . 110 ) 


_ iiP - 


N 


. -'!i 1_ ('iui - - 

2 L Lm, D, , V p, p,- /J 




Sq 


‘R,V 


H 


3y 


This is a second order, ordinary integrodif ferential equation. It 
is interesting to note that the limiting process has eliminatea the 
viscous dissipation terms. 

The stagnation line conservation equation, obtained from the 
thin shock layer equations, are a set of four ordinary differential 
equations in five unknowns, (i.e. p, v, H, P and C^). In addition 
to the conservation equations, the thermal and caloric equation of 
state is available to provide another independent equation. The 
global continuity equation was used to eliminate the tangential 
velocity gradient in the momentum equations and therefore is not 
needed in a solution of the stagnation line equations. It can be 
used £o^ priori to provide initial conditions for the thin shock 
layer equations. For a shock layer solution the rate of change of 
the pressure gradient in the tangential direction must be specified 
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TABLE 2.5 

STAGNATION LINE BOUNDARY LAYER EQUATIONS 


Global continuity; 


5u _ 2 A- 

Bx p Sy 


(pv) 


Species continuity; 

„„ 1 = . ( j . ) + cu, 

P By oY ^ ’ Y ^ 


X - MoniGntutn: 




- Momentum; 


Sy 


= 0 


Energy; 


BH 


^ k ^ + S h,J, 


By L“ " 3y 

T 


N. D. 

P 11 

N t i-J 


1 r.y 


ihi - ) 1 - 

n. D_. v' J Sy 




( 2 . 111 ) 


( 2 . 112 ) 


(2.113) 


(2.114) 


(2.115) 
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as a function of the normal ditactlon. Comment on how this mi6ht 
be specified is reserved until we have considered the reduction of 
the boundary layer equations to stagnation line equations. 

Dropping normal curvature effects and retaining only first 
order terms Eqs. 2.97. 2.99. 2.107. 2.109 and 2.110 reduce to the 
stagnation line boundary layer equations. These equations are 

presented in Tab. 2.5 


d^P 

^ = 0 and ^ = 0 for all y at x = 0, — 2 a constant 

dx 


Since 


and may be evaluated at any y station. If the boundary layer 
equations are evaluated over the whole shock layer as done at the 
stagnation line by Dirling. kigdon and Thomas Ref. 2.23 we may use 
the Rankine-Ilugoniot relations to determine this constant. 

rrom Eq. 2.61 the dimensional pressure behind the shock can be 

expressed as 


-X 2 ,, 2 

p = (1 - p) cos cp 
s 


(2.116) 


differentiating we get 




Bx 


j = - 2C1 - fJ’J' 


at X = 0, Cp = 0 by symmetry. Therefore 

2 

bk) 


( 2 ) = - 2(1 - p) Vsx) Q « 

bx ^x=0 


U 


(2.117) 


(2.118) 


In order to got the boundary layer momentum equation into a more 

(-Vio T'c|fe of change of the pressure gradient 
common form let us express the rate or cnange 
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in terras of the velocity gradient behind the shock. From Eq. 2.75 
the diraensional tangential velocity behind the shock is 


u =1^ sin cp cos e + p cos cp 


1 

sin e |U^ 


from which we can obtain 


_ r 1 

W L Sx P ox ,Jx=0 




) 


U 


(2.119) 


The rate of change of the pressure gradient, Eq. 2.118, can 
rewritten in terns of the velocity gradient. 

( L ( ^ 

If the shock is assumed to be concentric to the body at x = 0 

then 




= 1 


( 2 . 121 ) 


V=o ' x=0 

This gives a Newtonian velocity gradient used in many boundary 
analysis. Instead of applying this condition behind the shock most 
analyses apply this condition at the edge of the boundary layer 
which is at some intermediate station between the shock and body. 
Using the concentric assumption Eq. 2,113 may be written 
X - Momentum: (B L, S L) 




sy 


_ ,Su s2 

^ .e (1 ^ (pv))^ - 4p (1 - p) ® 

2 > 9y ^ ® x=0 


( 2 . 122 ) 
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It has been den,onstrated that the thin shock layer and boundary 
layer equations can be reduced to ordinary differential equations 
along the stagnation line ulthout resorting to slaiUarity trans- 
forruttions. By doing so one iirportant difference in the resulting 
two sets has becooe apparent. The stagnation line boundary layer 
equations are completely specified by boundary conditions at the 
surface and outer edge. However, an unknown function of y remains 
in the thin shock layer equations which cannot be determined, without 
approximation, by outer and inner boundary conditions. The undeter- 
mined function as stated previously is 

u2. 




( 2 . 123 ) 


This function like the rate of change of the shock angle is, by 
physical interpretation, determined by the flow downstream. The 
downstream flow is to be calculated by specifying these stagnation 
line conditions such that initial conditions may be determined. The 
problem is complicated further by the fact that there is no apparent 
theoretically based means of iterating on this function such that it 
could be assumed and corrected until some satisfactory convergence 
is obtained. The derivation of the stagnation line boundary layer 
equations demonstrates that to a first approximation the function 
F(y) is a constant which can be evaluated at the shock by specifying 
the shock geometry. For usual boundary layer problems the edge 
tangential velocity gradient is specified rather than the rate of 
change of the pressure gradient at the boundary layer edge. The 
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has been correlated as a function of flight 
velocity gradient has oecn ^ 

hr> hp used in blunt body 

conditions and body shspo for n,any oases 

• prHpr to specify this unknown downstream 
boundary layer solutions in order spec 

influence a priorj- . 

in shock layer solutions the shock wave has been considered 

p fs 23 24 2.9, 2.10, 2.17, 2.21, 2.22, 2.23 

concentric by Refs. 2.3, /.4, , 

t- nf these analyses set the function, 
and many others. Furthermore, most 

Eq. 2.123 equal to a constant. The full extent of influe 
assumptions has not been determined for radiation and ablation 
coupled ilous although some radiative coupled results are presente 
in Ref. 2.24. This is the point where engineering iudgement and or 
onporimental results must be used in order to mabe tbe matbematioal 

model useful- 

TO recapitulate the developments made In this section, it Is 
noted that an appropriate order of magnitude assessment of terms rn 
rhe conservation eguations was made. The radiative transfer eguatren 
was developed using a planar slab approximation. Using developed 
and stated expressions for the flux divergence, bulb viscosity and 
conductive flux the thin shock layer equations were stated retarnrng 
second order terms and first order terms with curvature effects. The 
first order shock layer equations without curvature effects were 
found to be the boundary layer equations. Shock and surface boundary 
conditions were developed. Subsequently, the stagnation Une equatrons 
were developed and discussed. This total development provides the 
required information to determining the appropriate equations to 
select for use in the present problem. 
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APT>1 TCATIO^LTO_R,A^il^^ 

• v,hirh provide a reasonable 

The selection of the proper equations 

• t'on of the processes in the shock layer is quite inrportan 
description of the p ^ ^ 

The detailed development of the s o 

, in the previous section provides a proper 
surface balance equations m the pr 

basis from .hich a selection can be made. _ 

^ a set of conditions oi 

X„ora«tca.ealogicalsc™, ase 

„,3ecUvesa.s..e.aana.. Baae. on xle caa.ea. sta. o. a.e a« 

a„a anticipated Icprovenents the foUovlng set o£ condtttons nere 

• to use in selecting the equations to be 
arrived at for a criterion to use 

solved. 

Cr 2 J^a_for_Sel^^ ^ .n^iption of the mass, 

.tioid an accurate eesuripti 
1 The equation j -- 

■ „entu., and energy transfer ulthln the shoch U.r lor 

• .nd flisht conditions of interest, 
the body size and tiig 

nrnst be valid for large and small mass 
2, The equations must 

injection rates fron the surface. 

3 . Xhe equations »ust be coupled property to the surface 

boundary conditions. 

• t detail ™at be maintained in the equations In 

4 . Sufficient detail muse 

.-i rafffsrts of diffusion and 
order to accurately access the effects 

finite rate chemistry. 

5 . The computer solution of these equations should be an 
engineering tool to predict surface heating rates and 
provide the means to determine ahlator behavior. 
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Based on the obiactlvc stated above the £lrst order shook layer 
equations without curvature effects through the shock layer were 
aelected to be solved. The reasoning leading to this choice follous-. 

second order ter»s In the conservation equations were dropped for 
the folowlng reasons. (D They do not significantly alter shock 

( 2 ') There are inherent approximation^, in th 
layer processes. K^) 

1 riof-a iiqpcl to predict radiation, 
theoretical models and experimental data 

thermodynamic, and transport properties. Thus there is no need for 
^re detail in the flou field equations than in the properties used 
in its calculation. (3) At high Reynold’s numbers, typical of 
Earth entry frou, bbtrs. .ost second order teres becoee less 
significant than at the Reynold's nunber used for the order analysts. 

nt-irr l’’n- mathematical rigor required dP/dy 0 
(4) Along the stagnation i^nc s.-i.--. 

t U T .i-inn nf the X - momentum equation from 

in order to uncouple the solution of the x 

itc t-. Stating dP/dy = 0 on the stagnation line 
downstream effects. Stating or/ y 

necessitates dropping second order terms from the y - momentu 
equation. (5) The addition of second order terms would not p 
numerical difficulties. However, they would add numerical computation 
tiu« thus detracting fret, the objective of an engineering tool. 

Curvature effects through the shock layer were not include 

U c fnr fhe large Reynold's numbers and 

the shock layer equations because for the large 

body radii of interest the shock layer is thin thus making the 
curvature effects second order, hikeuise. the ther„al diffusion ret. 
of the conductive heat flur vector was dropped because of its 


second order nature. 
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It has been argued that tor the large Reynold's numbers o£ 
interest viscous etfects could be neglected. Hovever, in order to 
assess tinitc rate chemistry and dlftusion etfects and in order to 
properly couple the flow field and the ablator through surfaces 
balances the controlling viscous terms isust be retained. 

The selected shock layer equations which are applicable for 
large as well as small mass Injection rates are presented in Tab. 

2.6. The.se equations are the boundary layer equations with the 
elception of the invlscld y - momentum equation. The Reynolds 
number limit of validity of these equations is shown in Fig. 2.7 
along with the stagnation line post shock conditions for a typical 
hyperbolic entry trajectory. For the 10-g undershoot trajectory of 
Ref. 2.26 the maximum heating occurs neai r-Lv. - 

trajectory levels out. Thus the equations are valid in the lange of 
interest. Conditions below flight velocities of 36.000 ft/sec are 
not considered in Fig. 2.7 and. In general, throughout the remainder 
of this work since the Apollo flight data is available for these 
lower flight velocities* 

The stagnation line equations which result from Eqs. 2.124 to 
2.128 for an axisymmetr ic body are; 

Global continuity: 


2 ( ) = - - (pv) 

V 5x P sy 


( 2 . 129 ) 
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TABLE 2 . 6 

BLUFF BODY FIRST ORDER SHOCK LAYER EQUATIONS 



Global continuity; 


Sx 



pu) 


+ r (pv) - 0 

w by 


Species Continuity: 


A Bx 
r 

w 

X - Koroentum: 

p" p" 


bx by 


bu 

by 


Momentum; 


2 _ ^ 

K pu - 


Energy: 


bH . M 

P“ f 8y 


i- ik- 2 ^ 

by \ 5 y 


- t ^ >'1 L5 


3 q 


R ^ b 


/ bu \ 


by by 


(2.124) 


(2.125) 


(2.126) 


(2.127) 


(2.128) 



Free Strei’.rr Density (Slug/ft ) 
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Species continuity: 


dC . A . , 

t — (J. ) + 

dT " 



y - monientuni: 


^ = 0 

ay 


Energy: 


dH 

P“dP 


d 


dq^ 


x„ order to .a.e the surface boundary conditions 
„ith the flow- field cuations the boundary conditions 
ain^lified usins the sa„e criterion. The ap.roprrate 
boundary conditions are: 

Species boundary condition: 

(pvc/ + CRx + - V ■ 

Elemental boundary conditions: 

c+ , + 


(pvC/ = (pvV -^^3 


(2.130) 


(2.131) 


(2.132) 

(2.133) 

compatible 

mus t be 
sur face 


(2.13A) 


(2.135) 
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Moroentum boundary condition; 


P+ = p" 


(2.136) 


Energy boundary condition: 


dy 




(2.137) 


, dT 


+ q 


+ 


‘R 


The species and elemental conservation equations are slightly more 
tractable than those presented previously since the internal diffusion 
effects have been dropped. Dropping second order terms significantly 
simplified the momentum and energy equation. The major complications 
remaining in these boundary conditions are the effects of external 

diffusion and surface reactions. 

In general Eqs. 2.129 to 2.132 must be used for proper ablator- 

flow- field coupling if the flexibility of arbitrary ablation rates 

is toaintained in the analyses. However, for the cases in which the 

ablation rate is large (i.c. at least (pv)^ApJJ^) •''5) 
simplifications can be made. If the ablation rate is large the 
convective terms at the surface are much larger than the diffusive 
terms, and consequently the diffusive terms may be dropped. 
Furthermore, in the absence of diffusion or surface erosion, the 
only surface reaction of significance is that of sublimation. In 
this situation the sublimation process can be computed in the 
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ablator 


analyses and thn effect account 


for in the (pvC.)‘ tern.. 


With these qualifications 
equation can be written. 
Species boundary conditi 

(pvC^)" = 

Elemental boundary condif 

Cp^3>" = 

Global mass balance; 

(pv) = Cp''^) 


a simplfff^^^ ae 


t of surface balance 


(2.138) 


(2.139) 


(2.140) 


tlomentum 


boundary condition; 


p" = P 


(2.141) 


+ 


(2.142) 


Energy boundary condition; 

Hfl' '"fl ■ "" 

.hese equations the™, for iar.e ahiation rates, the oniy c^-ate 

the flow- field and the ablator is m the energy 
coupling between the flow 

boundary condition. opction 

The shock boundary conditions derived xn the prev 

.re unaltered by the ,uaUf ications ^de in this section. 
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STAG NATIOlJil^^ 

- , K 1 n f f* body fi^st oirdBir 

I„ this section the stagnation Une , 

ntions arc nondtoensionaliacd and transformed into a 
shock layer equations are 

1 t-•r^n The difference between the 
forn, suitable for numerical solution. The 

. a i-'nit of the first section 

- t.. 1 nzation in this section and tnat ot 

nondimensionalization a. 

this Chapter is in the use of the post shoch density rat er 

J •. As a result a different Reynold s number 
the free stream density. 

appears in the momentum equation than appeared previous y 
llopment begins mith the dimensional .nations selected in 

nrevious section. 

using the foUouins nondimensional .entities . uhich are a 
„aed for the energy and elemental species equations, 


u = 


U 


V = 


CO 


* 


u 

CO 


P = 


JL- /V 9 

,0 c° 


X 

5 


•k 

= 

r'' 


h = 


C = 


P- = * 

^'s.o 

C" T 
_ P ^ l 


T = 


■* * 

k T 


;,o 


s,o 


(2.143) 




-'2 

00 


k = 


s ,o . 


-p u'- R p U 


*3 


J. = 
1 


S.o’co * * * 

J _ J- n U 

* * * ^^R 

Jg _ ^ ^ WhCrB E ~ • * 


Ps ^ 


dy 


U P ^ 

P o ^co CO 

2 129 and 2.130 can be rewritten in 
the continuity and momentum Eqs. . 

nondimensional form as; 

r.ontinuity ; 


2 


(pv) 

p sy 


(2.144) 
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Moment uin: 

u 

Re by V 

S 


oy 


\ b 

/bu\ 

_ 

(H) = ^ 

(2.145) 

)- 

fel) ■ 

^ \b§/ 

\ 2/ 



As developed in the first section of this chapter, the second 
derivative of pressure term in the momentum equation can be evaluated 
at the shock using the Rankine-Hugoniot equations rf — is assumed 


zero. This pressure term is 


2p _ 

= . 2p (1 - p) Va|) 

bfP 


(2.146) 


£or . concentric shock ^ " 1- The concentric shock assumption is 
usually rather good for hypersonic flow, however in this analysis 
S Will be treated as a paramter and thus will be left general. 
sLstltuting this tern, into the momentut, equation yields 


1 ^ f ^ (^)'\ - pv ^ - P (S) 


(2.147) 


+ 2 p(l - P) (H) = 0 


To obtain a more classic 
function is defined 


ical form of the momentum equation a velocity 


f' = ^ / 

J- rr ' 


bu / ,o _ lim 3 function of y 


(2.148) 


b§ ' 

^"s,o ^ bw ^ - be from the Rankine-Hugoniot equations. 


where " b§ b§ 


Substituting into the momentum equation yields 
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1 ^ 

Re By 
s 




2p (1 - = 0 

(Bu /B§) Vb^/ 

^ S ,0 


Thes equations can now be 


transf oi'med using 


transformation 


*n = 



Transforming Eqs. 


p 



6 

(2.144) , 


i- = - i- 

dy 


(2.148) with (2.150) 


Continuity: 


^ _ :i— 1 ^ (pv) 


S1\ 


B^ 6 

Velocity function 


- 1 


f^ = 


3L. 


Bu, 




B(pvl 
27, BTl = M. 

Bu . STl 

hT) 


Momentum: 




~ df , 
V Re^ 6 


Re 6 r 
s ; 

pCdUs^o'^°^^ ^ 



V B5 


(2.149) 


the Dorodnitsyn 


(2.150) 


yields 


(2.151) 


(2.152) 



(2.153) 
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Integrating Eq. (2.152) yields 


pv 


■ - 2 hr ) '> ' 


Equation (2.144) can 


be used to eliminate pv in 


(2.154) 

the momentum equation. 


This yields a third order nonlinear 
Momentum; 

/Su 


ordinary differential equation 


(p P t') J 


(2.155) 


Re r Mlji-Pi (sey - 

L p ^ a. 

^ s , o 




The boundary conditions for this equati 

Tl=l> 


”f\ = 0, • 


= 0, 


f' = 0 


(2.156) 


-(pv), 




in addition, Eq. (2.14) has a boundary condition impos 
the transformed standoff distance, 'K • 


ed on it to 


determine 


-(pv) 


T1 = 1, 


f = f 


2 '6 


The momentum equation can be reduce 
equation by defining 

f' 

C = — 


(2.157) 

first and a second order 


(2.158) 
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and substituting into Eq, (2.1 j>5) 


c" + '5' (-if > 


IC - ra?^i ^ 


2Rc^ "6 p (1 - p) ^ 


The resulting boundary conditions tor Eq. (2.159) are; 


(2.159) 


T| = 0 C = 0 


(2.160) 


= 1 £ = lA 

. ^vantifnrned in 3 tnsnneiT siroilsr to the 
The energy equation is trans 

ation First, however, the equation is written in terms 
momentum equation, tirsi, i 

of temperature as the dependent variable. It is anticipated tu.c 
„ 3 l„g tecpernture rather than enthalpy as the dependent variable will 
aave computer ti„e since both in ther»odynan.ic equilibriun. and 
chemical Unetio subroutines temperature Is used as the Independent 
variable. Thus the use o£ temperature ulll eliminate an iteration 
loop between the energy solution and the property subrouti 

The dimensional stagnation line energy equation is (noting the 

superscript * is omitted lor simplicity). 


dll _ d_ r_ ^ + 2 h 
- dy L dy 


1 125: 

i *^i J " dy 


consider the term on the left hand side of the above equation 


2 2,, 

dH , f h + - — ^ ) 

d;; = d^ 2 ; 
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Noting that 


h = r. C. h. 


dy 


dC. 


dh. 


““i o i- 

‘'i dT" ^ dT dy 


dC . ,T, 

V u i + il 2 C , C . 

^ "i dy dy i 


dC. 

1 

1 dy P dy 


i o <3T 

= £ h. -: — + C — 


(2.162) 


Thus 


■- ^ dT . d(v^/2) ~| 

pv i£h. — + Cp^+ dy J 

by b„bstit„tion of Eq. (2.162) and observing that u 
Substituting Eq. (2.163) into (2.161) yields 

dC . 2 (Jv 


(2.163) 

0 at X = 0. 


dy 


(2.164) 


+ fv ''di’ " 


dy t 1-- *^y 

Additional nianipulation of two terms is necessary. 


dC. 
1 


■^^i\ dT 


pv £ hi - pv £ h^ 1 ^^^ j dy 


where 


dc. 

d y 


dy 


££ = 0 


dy 



dCi\ -q 

.J ' Ty 'J "l(- 

P ®i dy J 


" ■ ^ oT dyJ 
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Notice the tern, (^) can only bn evaluated tor tlovn in local 

chemical equllibriun,. Since come calculationa uslns this assumption 
will be presented, tnrther discussion Is necessary. Substituting 
the above relations into Eq. (2.164) and rearranging yields 


SC, 




pv (^C + > h 


t— > » \ r^T 

V h —^ ] — 

L i ST ) dy 


(2.165) 


2 dv 


dq 


dy ^ dy 


■R 


The terms in the ( ) brackets above are from left to right the 

frozen plus re<acting thermal conductivity, k^, and the frozen plus 

reacting heat capacity, C , respectively. Using th._Et. de.--itic 

-T 


the stagnation line energy equation can be written as 


1- 

dy 


("t f ) - 


del 

dT _ 2 ^ . 

p^ dy “ dy dy 


■R 


(2.166) 


, . i: TE (0 bv the convension stated in 

Nondimensionalization of Eq. (.2. itibj oy me 

(2.143) results in 


dy 


(q- f ) - *--p. ^ 


2 ^ 

dy 


+ tE 


(2.167) 


Eq. (2.167) can now be trans 
(2.150) to give: 


formed using the Dorodnitzn transformation 


d I V 

d-n \P ^ dT^J 


■q pv C 


Pt 


(2.168) 


2 dv . 6 

‘ p'' 


IE 
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2 'T / 2d(pk ). 

AJ _ -A- pvc -:3 )d^ 

2 2pKj- VP Pt 6 dTl ^ 


(2.169) 


2 ^ ^ _iL tE 

‘ ^ p\ 

This later for.. Is second order and linear in temperature. 

The species equation nond Imenslonallzed usin. the convenslon 

stated In (2.143) retains the same form as Eq. (2.130) 


+ cu. 

dy 1 


(2.170) 


Transforming using Eq. (2.150) yields 


dJ. ~ 

3 I o 

— (i). 

dT\ P ^ 


(2.171) 


Por studies of hlnary diffusion and chemical equilihrlum, Eq. (2.171) 
ean be reduced to the elemental conservation equation using the 

following relations 


A .m.C. , 

~ _ V 1 = mass fraction of element j 

Cj = L M. 


" A..m.J. 


= mass flux of element j 


which yield 

c0. 

P'^ dif " dl] 


(2.174) 
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•„ T?n (2 174) since the net 
The geneeatior, tonn <!oe., not appear r„ Eq. (2.17 , 

Eeneratloe of elements Is zero. 

The thin shod laTor equations hevelopea In this section along 
„ith the corresponding boundary conditions are the primary se 

nur^rlcal methods used In the solution ol these stagnation line 
equations and results are presented In Chapter 5. 
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TRANSFOiU?KD AROUNT) T1!E BODY EQHATIOi'iS 

In order to complete the formulation of equations describing the 
hypersonic blunt body problem, the around the body equations are 
manipulated into a form suitable for numerical solution. This is 
accomplished by nondimensionalization and transformation of the 
thin shock equations (2.124) to (2.128). A von Mises transformation 
is used since it yields a set of equations amenable to forward 
integration and finite rate chemistry can be included using local 
one- d imens ional chemistry along streamlines. 

The shock layer equations (2.124) to (2.128) are rewritten here 
with some changes noted. The global continuity and two momentum 
equations are unchanged. The global continuity is removed from the 
species equation. The energy equation is witten with temperature 
as the principle dependent variable. Furthermore, the resulting 
form of the energy equation v;as obtained by subtracting the momentum 
equation, substituting the right hand side of the species equation 
for the left band side and assuming binary diffusion. The resulting 
set of dimensional axisymmetric shock layer equations are: 

Global continuity : 


5x 


(pur^) + r 


w By 


(pv) 


(2.124) 


Species continuity: 

SC 


Sx 


+ pv 


By By 


^ ('LtL 


BC., 


By 




pu 


(2.175) 



Ill 


X - rfiomentum; 


P“ ay 


BP , 9 / 


(2.126) 


- momenturo; 


H P 


By 


(2.127) 


Energy: 


pu 


: ^ 

■"p Bx 


+ P- f “ - ^ "i"! “ 


p oy 


Bx 


(2.176) 


+ ^ €) 

By oy 


Bq. 


‘R 


By 


Plfy) " ^ 


/ V Pi ^ 


L By J By 

i 


It should bo noted that the heat capacity (C^) and thermal 
conductivity (k) in the equation above consist of only the 
parts. This is unlike the final form of the energy equation derived 

in the previous section for the stagnation line. 

The shock layer equations can be nondimensionalised using the 
same set of nond Imensionalislng quantities used for the stagnation 
line equations (Eq. 2.1A3) with the addition of 


R 0)^ 


U). “ 


U 


* 


,0 oo 

With these quantities 
Global continuity: 


the preceding equations may be written as: 


B(pur^) 


+ r 


Blexl = 0 
w By 


(2.177) 
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Species continuity: 


aCi 5C. 

p" if + p'' ay 


J_ iL (_!i ] + oj. 

Le 5y \S„ oy / J- 

^ ^i 


X - iDOinentutn: 


Bu , ^ 

P" By 


y - momentum; 

2 BP 


BE Re By By/ 


K P 


By 


(2.178) 


(2.179) 


(2.180 


Energy: 


r 

BE 


+ 2 


By. 


H- 

r 2L - 
nuG — - 

^ p By 

= - Z h 

i'-*-'i ■' 

/ 

BT'\ 

- 2 K 

. 2P 

7cu'-.' 



Re 

s 

\oyJ 

V 

c 

Pi 





^c. 

By / 

By 



Bs 


(2.181) 


The transformation of the independent variables of the above 
equations from § and y to ? and Y respectively is accomplished as 


follox^js* Let 


(2.182) 


.(1) = -pvP„ 


(2.183) 


Consequently any dependent variable F can be written as 
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F = F(5, 7) 


f M. \ = ( ^ 
\ 



f 7)F \ f ^ \ 
( B7 V 


Substituting Eq. 2.183 into 2.184 yields 


W/ 


y ' 'i' 

Correspondingly 


Y Vby; 


§ 


p 


Uy)^ " Vos^ VoY/^ VSy;*^ 

from which, by using Eq. (2.182), one obtains 


Dur 


/ OFY _ _; '■i ,'iP \ 

KhyJ^ Y Vsy^y 


(2.184) 


(2.185) 


(2.186) 


(2.187) 


The nondimensional shock layer equations after undergoing trans- 
formation with Eqs. 2.182, 2.183, 2.185 and 2.187 are: 

Global continuity: 




_ + ^ f- Y-^') = 0 

hi hyJ Sy V S§/ 


(2.188) 


0 = 0 

Species continuity: 


SCi 


oC 


1 ;■ a_ 

^ “ Re Y BY VS dY / pu 


. U). 

i) + — 

/ pu 


(2.189) 
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X - TnomentuTn: 


" pu UY 


1 ^ fa 5u \ 

Re Y oY '\ SY/ 

q 


(2.190) 


y - monnentuni: 


^ H ^-1 

SY r 
" w 


Energy: 


C 


P SC 


1 ^ 
P SC 


SP , 2 ^ 
Y SY Y 2Y 




E h.o). 

1, I 


-f 2 



+ 


2 a 


Re ^ 
s 



f. ] 

'-^K L.S SY / St J 
Ci 


(2.191) 


(2.192) 


where 


pipur 

y" 


w 


(2.193) 


With the oxeeption of the pressure variation and radiative flux 
divergence tert., the preceding equations are in identical fort, to 
those successfully used to describe shear layer and cor-bustlon 

chamber flows with finite rate chemistry (Ref. 2.27) 

in addition to transformation of the shock layer equations, the 

necessary transformed boundary conditions are. 
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1. 

( 

Sy 

o 

1! 

•= 0 

from y = 0> 

u = 0 


/ 


\ 

p V r 
w w 

from y = 0, 

V = V , p = 

2. 

K 


^y=0 

Y 

w 

w ^ 


i 




from y = 6, 

U = Uj 

3. 

V 

Sy 

) “ 
y=6 


0 

4. 

( 

55 

li 

1! 

y 

6 

from y = 6, 

V = v^, p = 


5 ( — ') =0 


from y = 0, u - 0 


(2.194) 


6. T = T at y = 0 

w 


7. T = at y = 6 

Q 

8. C. = (C ) at y = 0 

1 1 w 

9. C. = (C ) at y = 6 

1 3-0 

The thermal and caloric equations of state are also required to 
complete the set of equations and boundary conditions. 
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CHAPTER 3 


GAS AND ABLATOR PROPERTIES 

In order to iraplemerit a solution to the shock layer equations 
radiative, thermodynamic, transport and ablator properties are 
required. The reliability o.f these properties naturally influences 
the reliability of the flow“ field solution. This chapter presents 
a discussion of the properties used in this work with special 
attention given to radiative properties. A method of determining rad- 
iative ly important species and species mechanisms is presented and 
applied to molecular species. Chemical equilibrium, thermodynamic, 
and transport properties of air and ablator species are discussed. 
Further, quasi- steady ablator response properties for phenolic 
nylon are presented and discussed. 

RADIATIVE PROPERTIES 

Flight conditions, atmospheric compositions and ablator composi- 
tions are the three important constraints which should b.e considered 
in determining what radiative properties are required. The flight 
conditions of interest yield shock layer conditions in the ranges 
(see Fig. 2.7) 

3000 < T < 17000°K 
.001 < 2 atm. 

Gases in these thermodynainic ranges are classified as low temperature 
plasmas. The gas is not highly ionized as in a full plasma state* 
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In a lo« tnlrparature plar,»a, .nolecules, ato»s and Ions s.ay play an 
important role in radiative energy transport and therefore the 
contribution and mechanism of these three types of species must be 
assessed. The atmospheric and ablator compositions define which 
specific species are present and the relative amounts of each 
present. Within this context a selection of the radiative properties 

used in this analysis is made, 

Bound-bonod tr.rnsUlc.ns: Line radiation results from electronic 

transitions between the discrete or hound energy levels in atoms 
molecules. For atoms the concern is with transitions between states 
nf differing principle quantum numbers within a given ionic state 
of the .rtomic species. These transitions are characterized by 
spectral line scries which are limited by the ionization threshold 
of the next ionic state (Ref. 3.1). The importance of a particular 
transition is measured by Us absorption cross section (Ref. 3.2) 


b-b 


X. 


% , tnc L , S.2 2J nn 


if , (1 - e 


-hv /kT 


) 


(3.1) 


which is a function of the £- number, f . and the line half- width 
y in addition to frequency and temperature. The f-number is the 
transition probability and can represent a single line, nmltrplet or 
a collection of line transitions having the same lower state. 

For molecules, the concern is with transitions between different 
electronic states of a given molecular species. The electronic 
transitions are modified by vibrational and rotational effects which 
produce groups of bands called band systems which are composed of 
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discrete line transitions (Ref. 3.1). The importance of a band 
system can be assessed in terms of the spectral cross section for 
the band. The absorption cross section can be written (Ref.' 3.3). 

r f (Acu.)"^ exp [- (arc«^)^/(AcUr)^] (3-2) 

U) o i. 

'where 

9 -13 

r = e /me 2,8 x 10 cm 
o 

f = band system f- number 
AiUj. -■= Aiu^[tanh (0^/2T)]^ 

Am = spectral half width at T = 0°K 
o 

9 = he wj' 

o G 

m " = effective vibrational frequency of the lower electronic 
e 

state 

u) = origin frequency or absorption maximum center of the 
o 

system 

This relation applies to both diatomic and polyatomic molecules as 
reported by Ref. 3.3. The important point here is to note that a 
bands effectiveness in contributing to radiative transport is 
measured in terms of an absorption cross section which is directly 

proportional to the band f- number. 

Bound- free transitions: Transfer between ionic states 

for both atomic and molecular species result in continuum 
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the initial or 

r„r a bound- free process, eithe 

r:::l fncfudes ^ .00 efeerron. Xbe absorpffon of radfefron 

afscrete afo^fc sfafe produced b, a pbofon bavfn. enousb 
from a discrete ,•« the 

threshold results m tne 

process of photoionization. 

(3.3) 

+ 

A + hv -♦ A. + G 

nrocess reccbinatlon. occurs when an ion and an 
The reverse process, 

electron recombine yielding an emission of a phot 

(3.4) 

A+ + e" -♦ A + hv 

u per state for both processes Is continuous, being defrned 
Stnee the upper 

by the kinetic energy absorp- 

.,0 ic continuous, rne ausu h 
A in the respective processes 

or emitted m the v j „„ CRef 3.4) 

,.n cross section for atomic species can be expressed <Kef. 

for hv 


b-f 


a 


= 0 


V 


for hv ^ 


f-b ^ b-f ^ 

V ^ 


G^,,(v)cos^[<pU.f' 

A'=Jt+l 1 1 + ev J 


where 


.-19 


.,2 


C(n) 
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n = effective quantum no 


. = [(z + 


X. 


= ionization energy of i-th state 


f the species and transition 


Z = charge on the residual ion 

G / cp are functions o 

li, »*•*'» 

g / = phase shift 

ij 

= threshold frequency 

Vt 

The hound- free traneUicns o£ »oleeules are not, in general, 
argnlticant enough to consider for this problen,. 

gnee-free transitions: Transitions hetueen tno free energy 

^ ^ 4- T 1 ^ p d 

levels in which free electrons a.c 

final state result In eontlnuu„ radiation. The es.ls=lon process 
IS cos»only hno„n as Bre^estrahlung and taUes place In the Coulo* 

field of an ion (Ref. 3.1) as 

+ _ (3*6) 

+ e" -♦ A + e + hv 


or in the field of a neutral atom as 

The free- free absorption process for an atom 
to the relation 


(3.7) 

takes place according 


A + e + hv -♦ A + e 


(3.8) 
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radiation from these processes may ha aons 

The concrnuum raoiaLio 

, a result of atoms or Ions deflectrng 

resulting In a 

rnciaent electrons and givlns them an acceleratro 

. The free- free cross section can 

photon absorption or emission. 

be expressed as (Ref. 3.4) 


„f-f.l. 26 xlO-'(Z + l)T^e 

V ^ 


(cm^) 


(3.9) 


where 

z = charge on the residual Ion 
X = ionization energy 

Q . partition function of the residual Ion 

^ ^ j ri- correction factor) averged 

— _ (a non-hydroge 1 <- 

g _ oeiC..- 

for a shell 

, with the bound- bound and bound- free the 
Por free- free process ,„,rtiveness 

in terms of the f nuni>er . 

«nny techniques have been used to obtain radiation 
3.31 and the final form of the data may be found in 
from intensity measurements of individual lines to total mt 

•,.Mires at various temperatures 

tneasured from entire gas m dictive 

Accordingly the theories which are used P 

pressures. Accor oj-ub i ..a 

orv in corresponding detail, as 

methods of radiant heating vary m 
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indicated in the previous discussion this »orh uses absorption 
cross sections as a basis for a radiative n.dei. Therefore the 
present discussion is centered around the «asures.nt of the required 
cross sections and the Informtlon needed to use n»dels for 
computing cross sections. 

Quantitative spectroscopic studies uhich yield absolute 
intensities and line half-uidth are required to determine f-nu„bers. 
cross-sections or absorption coefficients (Ref. 3.5 and 3.6). 
Measurements of spectral line profiles of individual lines may be 
„de in the microuave region houever integrated band intensities 
are sufficient for molecules. Additional data such as dissociation 
energies, ionisation energies and line frequency centers along »ith 
a quantum mechanical model for the species under consideration it 
needed to compute the cross section. Intensity measurements are 
usually carried out in isothersul experimental arrangements. For 
example loo temperature measurements can be carried out in constant 
volume cells attached to spectrometers (Ref. 3.5). For high 
temperature studies of continuum radiation shock wave induced 
plasmas and electrical discharges have been used to achieve the 

required thermodynamic state (Ref. 3.6). 

The purpose here is not to describe equipment. The interested 

reader is referred to references 3.5, 3.6 and 3.7 and their 
respective references. 

The selection of the species and their transport processes to 
include in the shock layer radiation calculation is based on the 
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relative contribution to the absorption coefficient. The 

• ed bv the elemental composition of 
species considered are determined by 

the atmosphere (air) and the elemental composition of t e a 

^henolic-nylon). Civen these elemental constraints and the 

te,eratnre and pressure ranges under consideration, thermody^m 
eonstralnts determine the magnitude of species molar composrt o. 

.he selection is divided into too proups (monotomic and polyatomic 

ptee and continui. mechanisms are the mapor contributors to 

3.S). Purthermore. the eiiccts of ionic species H . C . 0 . a 

^ ration- have been shown to be negligible (Ref. 3. ) 
and higher ionization, nav 

n 1 ?n 000°K The radiation model of Wilso 

for temperatures below 20,000 K. She four atomic 

(Kef. 3.2) includes line and continuum mechanisms of the fo 

• reasons and 

species and neglects the ionic species. ^ 

•rse tn his model was fairly simple, 
because adding molecular species 

it was selected as the basis of the present analysis. 

hpreement does not exist on which moiecuiat species and 
cespective mechanisms significantly contribute to the radiation 
rransputt in an ablation coupled shoch layer. The goal here is to 
3 „,,est a means of selecting the species and mechanisms which con 
he significant in the shoch layer radiative transport process. 
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Ablator and air molecular apeciea which exlat In thernK>dyna.ic 
equilibrium in slBnificant quantities ate presented in the folowins 
section. Restrictini; consideration to these species, It Is noted 
that the frequency dependent absorption coefficient is directly 
proportional to the species number density times the species absorp- 
tion cross-section, see Eq. B.14. Table 3.1 presents, for the 
molecular species, a list of estimated maximum molar compositions, 
maximum cross-sections and reported transition probabilities (f- 
number). This information is needed in order to estimate the 
relative effectiveness of a species and mechanism as a radiator. 

It should be noted that the table is not complete. Some Information 
either was not found or does not exist; while other Information 
was reported as being of questionable accuracy. Further, some 
molecular mechanisms were not listed since their cross sections or 
£-„umbers were negligibly small in comparison with the ones listed 

in Tab. 3.1. 

The effectiveness of a particular mechanism of a species is 
principly dependent upon its absorption cross section, ay number 
density and frequency range in which the absorption cross-section 
is the same order of magnitude as the maximum cross-section. Other 
factors which determine the significance of a mechanism of a 
species are the temperature level and the physical distance in the 
shock layer in which the species exist. Most molecules will exist 
in a near constant temperature layer near the body for large mass 
injection. Therefore evaluation of cross-section at and slightly 
above the sublimation temperature of carbon is realistic. Further, 



table 3 . 1 

! SPECIES CONSIDERED 
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TABLE 3.1 (Cont. 2) 
molecular species considered 
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photoionization 
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irrespective of the physical diinensions considered if the relative 
effectiveness of a tnechanisrn of a species is negligible considering 
the absorption cross- section , number density and frequency range 
that mechanism will not participate in the radiative exchange. 
Consequently a measure of the relative effectiveness of mechanisms 
of species in this work is defined as 



. 2 

(Y.) (Ahv) / (ev cm ) 
1 '^max J max J 


( 3 . 10 ) 


where 


ir. - maximum cross section of mechanism i of species j 

(Y ) = maximum possible mole fraction of species j in the 

j'^niax 

shock, layer 


(Ahv).^ = spectral interval over which o is greater than 

(^j^’^max- 

In the above definition mole fraction which is directly proportional 
to number density was used rather than number density. The radiative 
"effectiveness factor" defined above was computed for the species 
and their respective mechanism listed in Tab. 3.1. Figure 3.1 
presents a ranking of the larger effectiveness factors which vary 

over a range of four orders of magnitude. 

Based on the effectiveness factors in Fig. 3.1 a selection was 
tnade of which molecular species and respective mechanisms to include 
in the radiation calculation. The selection is indicated under the 



(Radiative Effectiveness Factor) (ev-ctn 
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10 
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C2(Swan) 
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) 

•CN(Violet) 22 u g 


*- H2(photoionization) 


-♦-CoCFreyma rJ^) q ^ (Fox- Her gber g) 


_^C^(A'rtu~^ 2g ^ 


.C„H(A'£ -X n) 2 + +N 

1 C2H(D S -A Z ) 


-»-C,H(B^A'-X^t) 





^ C2H2(B(?)-X Sg ) 

Other Species and 
Resp6Ctiv6 Processes 




FiE 3.1 A Comparative Ranking of Molecular 
Radiative Effectiveness Factors 
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eo^nt. o£ Ta.. 3.1. Based on thU selection the ^chanls.s - 
„,(Werner) , C^H^c'VX^.) and C^lKE^S^- 

radiation .odel of Wilson's. The hand average cross sections for 

all the species mechanisms along with the details of the r 

medal used in this uorh is presented in AppendU B. The radiation 

it Ice Hincussed and listed in Appendix 
model in co,n.puter program form rs drscus 

nori mAD 3 twelve species are considered 
C. In this program called LRAD 3 tweiie f 


Line and Continuum 


Continuum 


1 -lir. nne-' selected above. Computa- 

vhere the molecules listed are she one,. 

tionally the fre.uencg range 0 < hv s 30. Cev) is hroUen into 
tuelve continuum bands uith nine line groups used located at nine 

line centers. The boundary conditions on intensity entering the 

j ^r*nTTi the body s-Ttg 3.ssurnod 

shock layer from the free stream and from y 

aero. Thus it is assumed that no radiation enters the shock layer 
althrough radiant energy may leave through either surface, 
the assumption m,ade in state of the art analyses (Refs. 3.10, 3.19. 
3.22, 3.27. 3.31 and others). The only questionable part of this 
assulption appears to be in assuming the reradiatlon from the 
surface does not effect the shock layer. This is theoretically 
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correct for an optically thin shock layer. The results presented 
in Chapter 5 show that the magnitude of the reradiation flux is 
small in comparison to the other radiative flux in the shock layer. 
Further, the shock 'layer is observed to be optically thin in the 
lower frequency range where the surface radiation is emitted. 
Consequently the assumption appears valid. The magnitude of the 
energy loss from the surface is significant with respect to the 
energy absorbed and thus is accounted for in the surface balance 
equations presented in the third section of this chapter. 

One additional comment is pertinent regarding the selection of 
radiative properties. It has been assumed in the foregoing discussion 
that the ablator surface temperature is at the equilibrium subli- 
mation temperatuic: of the char and the a0..a..o.. ..espo-.d.^ _i. - q-a., 
steady manner. If however the surface temperature is significantly 
below the char sublimation temperature, which is the case in the 
early portion of a re-entry trajectory, higher molecular weight 
species will be introduced into the shock layer. The radiative 
properties of .such species should be considered. However, by 
and large the radiative characteristics of polyatomic hydrocarbons 
are unknown. In addition, if the surface does not sublime 
completely but rather mechanically erodes, which some experimental 
evidence indicates can and does occur (Ref. 3.16 and 3.17), then 
the radiative properties of solid carbon (soot) should be included in 
the radiative and thermodynamic calculations. The radiative 
importance of this process is indicated by the absorption cross- 
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section of soot for a 250 A particle as coniputed using 

data of Ref. 1.18) which Is a factor of approxtoately 10 larger 
than any gas phase species. Chemical equilibrium thermodynamics 
analyses do not predict the presence of solid carbon at temperatures 
above the char sublimation temperature. Furthermore, mechanisms 
which result in solid carbon in the flow-field are not well under- 
stood and consequently have not been modeled. Thus, for this work 
solid carbon was neglected although if it eaists in significant 
quantities it might well have a pronounced effect on the radiative 

heating. 

In addition to the detailed radiation model and associated 

properties just discussed, a simplified radiation model based on the 

1 rru 4 ^ H 1 OD^d bv 

. o - nr oir V;ay u^cu, 

emissj-un -- 

Ref. (3.19) is a correlation of the radiative flux divergence of 
air presented by Ref. (3.20). The curve fit equations for the 
radiative flux divergence are listed below. 


T^= 1000. logj^Q? + 13800. 


If T < T^, then 


logjo E =.0005 T + 1.15 log^Q? - 3.15 


If T s T^, then 


logpo E - log^(,P + 3.903 


(3.17) 
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In the prececding equations T is temperature in °K, P is pressure 

3 

in atm, and IE is the continuum flux divergence in watts/cm . This 
radiation model will be henceforth called the emiss . ^ mo^. It 
provides a computationally fast means of estimating the flux 
divergence in radiative coupled problems; although, as noted in 
Ref. (3.19), the emission model tends to overpredict the flux 
divergence obtained from more detailed models. 

The species radiative properties discussed along with the 
radiation model presented in Appendix B and the emission model are 
used in a numerical solution of the thin shock layer equations. The 
primary result from such a solution using either of the radiation 
models is the surface heating rate. An additional method was used 

, c — lTclt6S» 

in the present vork to compute ^ 5 

This method is based on a correlation of a radiative cooling para- 
meter presented by Ref. (3.21) which is based on the work of Page 
et. al. (Ref. 3.22). The radiative cooling parameter accounts 
for the radiative energy losses from an inviscid shock layer as 
predicted by the numerical calculations of Ref. 3.22. The cooling 
paramter , F, is defined as 


r = 


2(q.) 

isothermal 


( 3 . 12 ) 


% p U 

^ ’^00 CO 

The surface radiative flux is determined from the relation 


q^ = (0.2 - 0.295 logj^Q ^ 9 ^) 


for .04< r< 1.0 (3.13) 


’R ■ isothermal 

where the isothermal radiative flux must be determined from an 
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independent calculation. Tn the present work the isothermal flux 
is determined using the line and continuum radiation model presented 
in Appendicies B and C . This correlation for radiative heating is 
used in a computer program called RADCOR developed as part of the 
present work. Additional details concerning the use of the corre- 
lation along with the RADCOR Program are presented in Appendix E. 
Results and comparisons with more complete analysis are given in 

Chapter 5. 

To recopntolate, the molecular and atomic mechanisms which 
produce line and continuum radiation have been discussed in relation 
to the current problem. A method tor selecting the mechanisms 
which signitlcantly contribute to the radiative energy transfer has 
been developed. This method based on a radiative effectiveness 
factor was applied in selecting molecular mechanisms. The equations 
for an emission model of air was stated. Equations for the radiative 
cooling parameter were stated. These equations were used in the 
RADCOR computer program for simplified heating rate calculations. 

THERMODYNAIIIC PROPERTIES 

Two sets of thermodynamic properties and computational methods 
were used in the present work. One method is limited to 
species of air whereas the second is applicable to arbitrary gas 
mixtures. The method which is limited to air species was used, 
because of the computational speed, in solving the Rankine-Hugoniot 
equations .ind in computing thermodynamic properties of shock layer 
gases where no ablation products exist. The general method was 
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used to compute thermodynamic properties of air-ablation and ablation 
product gcis mixtures. 

Equilibrium thermodynamic properties of air which includes its 
major components (0^. N^, 0, N, 0+, N+. e") were obtained in closed 
form from a set of approximate partition functions using the method 
of Hansen's (Ref. 3.23). The partition functions and thermodynamic 
relations for an ideal gas were programmed (Ref. 3.24) such that the 


total enthalpy, entropy, speed of sound, average molecular weight, 
heat capacity and species concentrations can be found for a given 
temperature and pressure. The following section discusses the 
transport properties which were also computed using Hansen's method. 
The computational scheme for both the thermodynamic and transport 
properties uf air is embodied in subroutine GAS which is listed in 
Appendix D. Figure 3.2 presents a comparison of species number 
densities computed using this method with number densities from 
two different free energy minimization methods. This agreement is 
typical and lends validity to the use of the species number densities 
as well as the overall thermodynamic properties obtained from the 
approximate partition functions of Reference 3.23. 

Thermodynamic properties for arbitrary gas mixtures were 
computed using a free energy minimization program, CHEMEQ, developed 
by Ref. (3.26) which is a modification of a program reported by 
Ref. (3.25). Curve-fits of thermodynamic data for individual species 
are required as the basic information for a calculation. The forms 
of the curve- fit for the standard state properties are listed in 
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Tab. 3 . 2 . From the pure component properties vhtch arc represented 
by the curve-fit equations the mixture thermodynamic properties ore 
computed by minlmlring the free energy of the system. 

The selection of the species to include in the shock layer 
calculations was based on the ablator composition, the shock layer 
temperture, and the pressure ranges of interest. Air and the 
phenolic ablator selected for study contains the four elements, 

H. C. 0, and N. The temperature range considered was from the sub- 
limation temperature of the ablator up to approximately 17000“K and 
the pressure range considered was from .001 to 2.0 atmospheres. 

Twenty species were found to have a significant concentration in 
temperature range of interest at 1.0 atmosphere pressure. This was 
taken as represee t.at ive in the pressure range. These species are listed 
in Tab. 3.3. from Ref. 3.32 along with their respective heats of forma- 
tions where the reference elements are Hj. Ny, Oy. solid carbon and e'. 

The curve-fit constants obtained from Ref. 3.32 for these 
twenty species are listed in the block data package of V1SR.AD 3 
presented in Appendix D. Two sets of curve-fits are used for the two 
temperture ranges 1000“K to 6000"K to 17000"K. Details of the 
Fortran nomenclature for the curve-fit equations are given in 
Appendix D. 
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TABLE 3.2 

A SUMMARY POLYNOMIAL EQUATIONS FOR STANDARD 
THERMODYNAMIC PROPERTIES 


Specific Heat 


= A^ + A^T + A^T^ + A^T^ + A^T^ 


Enthalpy 


^2 ^3 2 '^4 3 ^5 4 ^6 

i = A^ + 2^T+3^T +/t +-T +- 


Entropy 


■^3 2 ^4 3 ■^5 4 

. A^lnT + A^T + f T + + A^ 


Free Energy 


F. 

RT 


T ^2 ^3 2 \ 3 ^5 4 "^6 

= A.(l-lnT) - — T - — T - - 2 qT + .j. 


- A. 



1A2 


TABLE 3.3 

COMPONENTS AND HEATS OF FORMATION 



0 

(AH f)298.16 

(Ah/)o 

Conrponent 

(k/cal/emole) 

kcal/ j^mole 

N 

112.951 

112.507 

0 

59.544 

58.972 

C(gas) 

171.301 

169.990 

N+ 

449.709 

447.564 

0+ 

374.867 

372.942 


432.357 

429.537 

H 

52.098 

51.620 

C- 

197.0259 

195.000 

z 

c_ 

189.6115 

188.000 

3 

C 2 H 

117.6448 

116.700 


53.8670 

54.300 

CO 

-26.4179 

-27.202 

CN 

109.7865 

109.000 

HCN 

31.1895 

31.281 

154.000 

C.H 

4 



127.100 

C^H 


0.0 

0.0 

No 

0.0 

0.0 
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trarspqrt properties 

The transport properties used in the present work 
Hansen (Ref. 3.23) for air. As pointed out by Hansen the knowledge 
of transport properties of air at high ten^eratures is in an 
ele^ntary state as cohered to the knowledge of the ther.odynas.ic 
properties. The basic Information needed for calculation of 
transport properties from kinetic theory xs the species 
potential such that the collision integral may he computed. At 
high temperatures the atom-atom, atom- ion and atostmolecule potentials 
are required but are not known well enough for accurate transport 
ealculations. The problem is further compounded when ablation 
products are introduced into the gas mixture. Thus as an engineering 
approximation, the closed form expressions for transport coefficients 
(r.e. thermal conductivity and viscosity) for air were used over the 

entire shock layer. 

The calculation of transport properties by the Hansen method is 
based on simple kinetic theory of hard spheres. The viscosity is 
calculated by a simple sunn.atlon formula (Eq. 66 of Ret. 3.23) 
which is a function of species density, mean velocity and mean free 
path. These variables are then related to the temperature, pressure 

• 4 -’ tnf mixture The thermal con<3uctivity 

and air species composition of the mixture. 

is calculated as the sum of two effects. The first effect accounts 
for the energy transfer by molecular collision which is processed 
for ordinary thermal conductivity. The second effect accounts for 
energy transfer by diffusion of species which are involved in 



equilibrium reactions at each point in the flow-field. The first 
effect is expressed as a function of species density, mean velocity, 
mean freepath, molecular weight and heat capacity. The second is 
written in as a function of the temperature, species diffusion 
coefficients, concentration and stoichiometric coefficient as well 
as the equilibrium constant. These two effects are the frozen and 
reacting parts respectively of the thermal conductivity discussed 
in Chapter 2 (Eq. 2.169) when developing the temperature form of the 
energy equation. 

Using air transport properties is justified not only because 
significant uncertainty is present in the basic data, but, more 
importantly, the resulting heating rate from a shock layer solution 
is essentially unelfected by the difference in air and air-ablation 
product transport properties. Rigdon et. al. (Ref. 3.27) showed for 
a typical flight condition a 1.5% change in radiating heating rate 
as a result of the difference in air and air-ablation product 
transport properties. More recently Esch (Ref. 3.26) has further 
substantuated the negligible effect of different transport properties 
on radiative heating. Specific calculations substantiating this 
assumption will be given in Chapter 5. 

ABLATOR RESPONSE PROPERTIES 

The present work uses a quasi-steady ablator analysis of a 
phenolic- nylon composite ablator. Furthermore, the surface tempera- 
ture is selected to be the sublimation temperature of the char. These 
restrictions imply appreciable ablation rates (at least 5% of the 
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free-stream mass flux); the surface balance equations (2.138 to 
2.142) are appropriate to describe such flow- field ablator inter- 
action. Results presented in Chapter 5 will show post prioiii that 
this is a valid approximation for most flight conditions of interest. 

Since the elemental species, total mass flow rate and pressure 
are identical across the interface for quasi-steady response with 
appreciable ablation rates, only the energy boundary conditions is 

considered here. 


, dT\ 

\ dy dy 


(2.142) 


The left hand side of the above equation represents the heat flux 
into the body. The quasi-steady model of the ablator can be used 
to relate this influx of energy to the heat absorbed and converted 
away. This relation is 


\ fy - 


(3.14) 


CoiLbinins Equations (3.14) and (2.142) yields 


dX"^ + 

(pv) AH - dy “ ‘^R 


(3.15) 


or 

„,4 (3.16) 

(pv) AH^^ “ '^T “ 

where q.j. is the total heat flux (convective plus radiative to the 
surface from the flow- field and AH^^^ is the heat of ablation. 

The phenolic-nylon composite considered in this work is 



146 


1 ,-n Rpf 3 28, This composite has a mass 

described by Pike et. al. m Ret. i.^b- 

,, .,.„bon. ;.?W hydrogen. 4.96*4 nitroECn 

elemental coroposition of 
and 14.724 oxygen. 

The heat o£ ablation o£ phenolic-nylon »as computed hy balho££ 

(Ref 3.29) using a checlcal equilibrium program developed by 
(Ref. 3.26, uhtch accounts for solids as mell as gas co„.onents. The 
heat of ablation uhich includes the heat of pyrolysis (300. BTU/lb„) 
is presented in Pig. 3.3 along .1th the computed sublimation 
temperature. The computed sublimation temperature and heat of ablation 
data .ere least squares fit as a function of pressure for compu- 
tation work. The equations from the curve fittrng process are 


O/. nr, /, 4- 1R7.0 log,„P + 9.^.3 (logj^QT) 

sub 


gH -9179.7- 114.81 log^„P + 329.64 (log^„P)' 

ab 


(3.17) 


(3.18) 


bslng E,s. 3.17 and 3.18, the surface halance equation (3.16) 

„as solved for various pressure levels providing the results sho.n 
in Fig. 3.4. The results indicate that the ablation rate Is a .eak 
function of pressure and a .eaker function of the sublimation 
tes,.erature for a given flow-field heating rate and constant 
surface emlsslvity. A conservative value of emissivlty is 
.6; such as that reported by Pope (Ref. 3.30) for phenoUc-nylon ebar . 
The average value reported In Ref. 3.30 was .67. this is in contrast 
to an average value of .95 for the data reported earlier by 



AH , (rTU/lbm) 



Phenolic Hylon Heat of Ablation 
anrsublimatlon Ten,peratare as 
a Function of Pressure. 


SubliiT'.ation Temperatur 
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.. .K. 

;ing two value 

calculations were made for tvo v.iuu. ■ - - ' 


ing two values 


■ ■ -tv namely: 1.0 and 0.4. The results are shown an 

of enu-ssivity, nameiy. , , n 

f nsf> lb /ft^ sec for 1.0 
FIB 3.5 and Indicate an uncertainty o£ .0. 

„ 031 lb /£t^ sec at 0.1 atmospheres. These values 

atmospheres and . „ „^aentage uncertainty 

are independent o£ heating rate; however the percen 

increases with decreasing heating rate. 

The results presented in Fig. 3.4 are used in Chapter 5 to 

obtain coupled flow- field and ablator solutions. 



Total Heating 


Phenolic Nyloi^ 

j. tr 6 


0-1 atm 




Ablation 


Rate (Ibm/it 


-- region 

for Various Pressure 



Ablation Rate (Ibm/ft^ - sec) 


3,5 The Effect of 

Surface hmissrvi-uy 



151 


CHAPTER 3 

Ref£T?ncP'-. 


Kltievieir Pub, Co.) 1 >' 68 , 
, reoffrey. Plasma SpecW sco£y. El 

j.l Marr, Geoffrey. Coupled Viscous- 

V H "Stagnation Point Ana ^^ 5 ;^ CR- 1548 , 

3.2 Wilson^^K^^_^ Flo/with Massive Blowing, 

June 1970 . 

"Equilibrium Opacities 1 

V ndW E Hicolet, "Spectral Absorption^^^^^ 

3 ., Wilso^K.^a^^^ 1967 . 

Radiat. Transfer, Vol. /. PP 


1.5 

3.6 

3.7 

3.8 


Radiat. , 

, S S 

. ,-o„ 

. „ "StudieB of 

cairns, R.d- and J. a- ^ Sitrogcr.," MASil r.n- .9 . 

by Atomic Hydrogen, Ox>oen, 

March 1968. Jomir 

, S s and D. Clfe, Radl«irfMX. 

Penner , ^ ^ • 

Press ) 1968 . 

Hicolet, W. E.. Sntfimrnlfe'Lrndr^^^ 

Transport m Ablation 

Layers," NASA CR- 1656. Sept. 

F Lasher, "Convective and Radiative 

”-i:ffVa;sffr «ee. 

,.„e.ed-U,.oc.dndd« ssr 

1 1 H and K, H. Wilson, t m^c Report No, 

3.10 Boshizaki, H. and Entry.' ^EC Kepo 

Transfer During P a,- _-^ies and Space Co., June 
4-43-65-5, Lockheed Missiles 

H G J T. Ohrenberger and T. R. Thompson^ Em^ 

•Ab;;rption of Badlant - a ^ . 3,03. 

, , n atAA Journal, voi. ->> 

Atmosphere, 

A T Wilson, "Measurements o - 

•n A A Textoris and d* tti troszen and 

3-u and f---rnt‘:“s;ir^^* 



152 


3.13 


3.14 


3.13 


X, "Fxoeri mental Measurements of 

ThoMS. G. M. and M. Ldintion from Planetary 

Nonequilibrrum and EquUr^ ^ ^ ^ ^ 1966. 

Atmospheres," AT'^^ 

^ V, L woouwata, 'Studies of 
Arnold, J. 0., V. H. Reis Entering Planetary 

Shock Layer Radiation o ^1 1965. 

Atmospheres," AlAA_Jour n^l . Vol. J. 

Golden, S. A., Stolecu^^^ 

,,;etronic ^ 7, PP. 223-230, 1967. 

Snectrosc^_MHi^i=-:— i 

^ . 1 AKI 


3.16 

3.17 

3.18 

3.19 

3.20 
3.21 


n tt n K "An Experimental Study of a Carbon-PhenoUc a- 

'"“tiin »te^ial," NASA ™ D-5930, Sept. 1970. 

• nf Carbon Particles Produced 

:\?eL?Mrtor:"'rinri''Ke^ort GD/C-I.BE66-006 . General 
Dynamics Corp., May 1966. 

W Snradley, "Radiation Absorption Effects 
Engel, C. D. ^nmuncered During Hyperbolic ^^n_^ry , 

on neaiiug Lb.ere... --- ^ ^ 1969, pp. /oa-<t>o. 

j Spacecraft Rue ., \ ol . , 

Avco Corp. 

,-n- "Planetary Entry Body Heating 

Livingston, P. and J . i ra ^ Atmospheric Gas up to 

r-^5ToO»K:"“^. 


3.22 

3.23 

3.24 


3.25 


’X = iJUUW , 

ij T Borucki, and D. L. Clifforie, 
Page, W. G., D. b. Compton, ‘ Monadiabatic Stagnation- 

"Radiative transport in 68-784, June 1968. 

Region Shock Layers," AlAA Paper 

• i-'nnt; for the Thermodynamic Properties 

Hansen, C. F. ’’APP^-^ir " TR R-50. 1959. 

of High Temperature Air, 

J p n Pnpel "A Com.puter Program for 
Snradley, L. W. and C. D. Enge , Radiative Heat Transfer 

predicting Coupled Reentry," LMiSC/HREC 

to a Blunt Body rsp^ca Co. . Huntoville. Ala., 

A791350, Lockheed Hissiies f 

Hay 1968. 

Stroud, C. W. and 

Ablation Materials Including 

D-5391, August 1969. 



153 


3.26 

3.27 

3.28 

3.29 

3.30 

3.31 

3.32 


Escb D D. . staCTation Iteatini; of n rhe..ollc-Ky .loa 

La .. , >vag ^ ^ • ' 

. n c. R B Dirling, Jr. and M. Thomas. "Stagnation 
PpInt'He;^ l.ypervelccity Atmospheric Entry, 

NASA CR-1462, Feb. 19/0. 

r-. A A T-i 1 F n del Valle and S. Hacker, On 
Pike, R. W. 5 G, A. Ap , • • ... pvrolvsis Products 

Methods of Determining Co^posrtron f (1970) 

from Ablative Composites, Vol. , 

Balholf, J. F., Research Associate, Louisiana State University, 
Private Communication. 

r Total Surface Emittance of 

Esch. D. D., A. Siripong. R H. PlRe. ^“I^rogon and 
Properties in T, KFL TR-7C-3, 

S. Ocpt., Eoniaiana State University. Nov. 1970. 



QIAPTER 4 

'Jill T,''^ T r A T PP oc edures 

and energy equations 

... stagnation Una continnity. .o»antn„ en .y 

1 a hv finite differences using a 
j -rt rhaoter 2 are solved by fm 

selected in Chap ^j ,^3 

tini called VISRAD 3 deveoped as par 
eocputar program 

ahaptai desctihas the techniques used in VISKAD 

numerical solutions. herein has been 

developed nerexu 

. ,• ii finite difference method aev 

The implicit fmi 

1 fRef 4.1), Blottner (Ref. 

^ a in oart by Adams et. al. ( 

presented m p . 

n • fRef 44), Edelroan and Hoffman ( 

4.3), Bavis CRcf. 4. ,P,,erential 

. .. /n,.f /, T). The stagnation tme 

Pay ana nay „arlvatlves are replaced by three 

equations are quaslUnearia . ^ 

point unite diUerences, An l.pUcit set ol 

^ 11 from the previous steps for each 

of the tridiagonal for.> result from P 

he rapidly solved psing an 
Thiq set can be rapa-ux^ 

second order equatio . ^ 

algorithm for tridiagonal matrice . 
implicit in the shock layer coordinate. 

Xn Chapter 2 the continuity and momentum equations were 

.ransformed and split into the following first and second or er 
ordinary differential equations. 

f' 

c = r’ 


(2.158) 


and 
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Cpv^) c" + L" 


C 


+ (piO' ;C 


^3 

Re 6 
s 


(^■)c' 

V c)^ 2 


(2.159) 


- i' ' 


Si 


for Eq. (2.159) are: 


The boundary coi^ditrons 


T\ = 0 


- 1 


c - « 

r -- 1^6 


(2.160) 


Fa. 2.159, is second order and 

The tnoroentuir, equatxon, . • • • it is desirable to 

„„„Une„. X„ .o .0... - 

obtain a lineal seconl ordei equation 


a\i dw , _ 

5-i. + a 't:; ^ “o” 
2 1 II 

dT^ 


(4.1) 


"3 


tfa f 4.2 and 

,a n is the lineal foin. used by Blottner (Ref. 

(Lf qM.H-l.an and H„ffe.n (Refs. ..WU-.O, 

„ £ 4 7) and others to solve see or all of the 

,nd pay and Kaye ( e . completeness. The 

conservation equations »i ^ manner 

. (2 159) can be quasilineatized in 

nonlinear term m Eq. (2.1^ 
of Lee (Ref. 4.8) . 


(0 = (c') " ^ ^ ^ 


(4.2) 


number Substituting Eq. (^-2) 

where k is the iteratron number. 

. . ^ -Linear equation of the form of Eq. ( • 

(2.159) yields a linear h 
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c' + L pu t 


S5 


PM- 


pp. 




(4.3) 


.V pd - P) I [- j 

2 \“2^ ^ + 2p^x U| J ^ J 


s ,o 

In Eq. (4.3) the superscript k denotes, computed from the k-th 
iteration, and the superscript k+1 for the current iteration has 

been dropped. 

The three point, variable step size finite differences stated by 
Davis (Ref. 4.4) and given below were used in Eq. (4.3), 

^ Vi 


/dw\ 


' A Vi \ + 4 VP 


^ 'Vi Vi 


(4.4a) 


\ 


Vi'* \ Vi> 


= r w , , 

V" *V‘\ + ‘Vi> 

. 2 „ (4.4b) 

i T)„ 4 Vl " 

4 Vl^‘ \ + * VP 

Using the above relations for the difference equations, Eq. 2.23 


can be written as 
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. 2- ^ ^ \ • 


1 •" 


|- a^(A \ - & Vi^ 

I- 6 \ ^ Vl ^ 


+ a 


i]Cn 


(4.5) 


^ 2 + a^ A T\j^_ 

+ li TV,(4 H„ + A V 

Equation (4.5) is of the form 

A C„.l + » C„ + «n+l ■ ® 


(4.6) 


„hlch yields a tridiagonal matrix in the vector equation to he 
solved. Equation (4.6) can be »rltten in the following matrix form: 







Cl' 



“ 

'^1 Co 

^2 





C 2 


^2 




.“3 <=3 



C 3 

— 






A B 
n n 


Cn 



<=„ Cnil 


(4.7) 


Equation (4.7) was solved using the tridiagonal Inversion algorithm 
program given by Conte (Ref. 4,9). This algorithm gives a rapid 

and accurate solution of the tridiagonal system. Ho significant 
Improvement In accuracy has been noted in the use of double precision 

in calculating variables in this subroutine. 

The first order equation (2.158) which accompanies Eq. (2.159) 

was integrated 


f='6j 

‘*0 


(4.8) 
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where 


' ( ov*) 

2 '6 


(4.9) 


using a simple trapezoidal scheme. 

The transforaed standoff distance,'?. «as computed by using the 

boundary condition 

-(pv)^ 


T^=l, 


( 4 . 10 ) 


2 % 


and Eq. (4.8) in the following relation 


f - f 

s w 

6 


(4.11) 


,1 


0 


j co^-oared to the assumed 6 i-or con- 

This computed value wat, c.U«-n co...paiut. 

t-»QTT ^ 1 c: a<^sumed and the entire set of 
vergence. If necessary a new 6 is assumea au 

equations is solved again. 

Finally the actual standoff distance is computed using 

(4.12) 

6 ='K [ P d T1 

‘‘o 

with a simple trapezoidal scheme. 

The transformed energy equation developed in Chapter 2 


_ -I— fpvC 

2 2pltj Pn 


KpKpK dT 


dvi 


dT\ ) dT] 


(2.169) 


~ 2 . 

6 V 

k„ dT| 2. 


lE 


"T p ^ 

is linear in temperature and thus quasilinearization of terms is 
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not neodea. In order to nee the above equation in the fort, of Eq. 
o- .1 , -r pviv, IF must be a known function of 11. 

Numerically this is accomplished by determining the flux 
profile fro» the previous iteration and this profile is then used 
in the solution of the energ, equation. Additional information 
concerning the iteration technique is given later in this chapter in 

the discussion related to Fig. 4.8. 

The boundary conditions for the proceeding equation are: 


11 = 0 , 


X = T 


w 


(4.13) 




T = T 


. - K firs pihlAtcr -c-nonse (for ablation coupled 

where T is deterrainea by the ablator - 

, ” ^ tp is chosen) and T is determined by the 

analysis sublimation 

solution of the Rankine-Hugoneot equations. 

Following the development of the momentum equation, finite 
differences (4.4) can be substituted into Eq. (2.169) for all 
derivatives. This yields 
2-6 


'l-‘Vh‘’V ‘Vh -’'"' '' 


+ L 6H„6Vl 


-fn 


(4.14) 


b>\<ATln + ‘Vd-’ ^ 


where 
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*1 " 2pKr 


d(pKj,)-, 

' J 


(4.15) 


^ 2 

„ 1_ ^ ^ — E 

°3 " ’S' p\ 

Equations (4.W) and (4.15) yield a tridlagonal Mtrix like Eq. 

(4.7) »hich is solved using Che algorithm described for the momentum 

equation, 

in the present analysis a simplified species equation solution 
was used. It was assumed that the elemental composition was constant 

and equal to the ablator elemental composition from the ablator 

• ^-1 ^ \7=0^ Likewise from the stag- 

surface to the stagnation point (i.e. v u;. 

nation point to the shock the elemental composition of air was 
assum.ed. Results using this approach are in very close agreement 
with results containing solutions of the elemental species equation 
as demonstrated in Chapter 5. If the elemental species equation 
was to be solved numerically, the same procedure discussed for the 
Other equations can be used. 

In order to take advantage of the variable step size finite 

£ •---a-m Tiac apvelooed to determine the step 

differences a set of criteria was deveiopeu 

aite used in various regions of the shock layer. The nondimensional 
temperature was used as the variable to specify the step srze 
pattern since in general it exhibits regions of more rapid change 
as a function of 1, than the velocity function. Numerical experi- 
mentation resulted in the following sufficient conditions : 
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,00. < 0T< .0. (TKC noo.i^nslo„al tc.eratv.re change 

^ t; than 005 unless condition 
a sren -an not be less than .UUD 

applies ot greater U«au . U5 . 

.V preater than .OA; no lower 
2. -0^ (The step-size rs never g 

limit is set.) 

t.n is used to maintain accuracy xn the 
The upper limit on step- size x 

The H step-size is updated between 
-intejiration. ine 
„entur, equatvoa integ 

. energy solutions as the temperature 
momentum and energy 

convergence. t-taa in the VISRAD 

.he numericai methods discussed s,ere implemented 

„hich is documented in Appendix B. Ihie 

3 computer program 

,lso contains additional information on the 

. , , of the stagnation line equations, 

used in the 

SOLUTION liEHAVm 

-—Tin this section are presented to demonstrate the 
validity of the momentum and energy equation selected to illustr 

f aonlicability of the numerical methods used 
the wide range of applica y 

TO flchicvGci under 

to exhibit the manner in which convergence 
different situations. 

fs the results of four different methods used 
Figure 4.1 presents the tesu 

to solve the shoch layer momentum and continuity equations. 

thods agree reasonably well for this case in which there rs 
four methods agree 

no mass infection and no radiation coupling. The present metho ^ 
and the Adams-MouUon predictor corrector method used by Howe an 



(du/d^) /(dUg/dj;) (velocity function) 


1.0 



y/ 5 , Nondirnor.i Lonsl Shock Layer Coordinate 


Fig. 4.1 Comparison of Stagnation Line Momentum Equation 
Solutions for the Case of No Radiation Coupling. 
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Vegas (Ref. 4.10) agree quite well. Both of these methods as well 
as the integral method of Spradley and Engel (Ref. 4.11) include 

, ;e method reported by Wilson 

(Ref. 4.12) does not include this effect and it is probable that this 
effect accounts for the lower values of the velocity predicted by 
Wilson’s method. The deviation of Spradley and Engel’s integral 
solution from the present finite difference method may be accounted 
for by two effects. First, and foremost is that integral solution 
is limited by the degree of the polynomial selected in its 
implementation. Secondly, second order effects are included in the 
momentum equation solved using the integral approximation. A 
comparison of computed stand-off distances between the present 
implicit method and the integral method for the case shown in Fig. 

4.1 indicates agreement within 67c. 

Figure 4.2 presents a comparison of velocity profiles with 
variable and zero d(p|j,)dT| for zero and 20 percent mass injection 
(i.e. (pv) /(pv) - 0 and .2) . The results indicate that for no 

mass injection, as anticipated by Wilson (Ref. 4.12), neglecting this 
effect significantly changes the velocity profile near the wall. 
Further, the results indicate that for cases when the boundary layer 
is blown from the wall by large ablation rates the effect of this 
term is not as significant in determining the velocity profile. 


* 


Note the negative sign of (pv)^/(pv)g has been dropped. 



7], Transformed Shock Layer Coordinate 


Fig. 4.2 Conrparison of the Effects of (pp.)^ on the 

Momentum Equation for 0 and 20 Percent Mass 
Injection 
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4.3 presents e comparison o£ temperature profiles co„u- 

rea „itl tte Present -I'™* 

i , I I ef. 4.11)* agreement is 

tiOli ■•■ '■'•^ ,.)^, I. c.t. i c I 

..nPte ,oo.. Both solutions use Hansen's thermoHynamic and transport 
properties. Thus the only factors contributing to the dtfferences 
p„ the temperature profiles are the numerical teehnipues used, the 
effects o£ second order terms in the equations solved by the 
relaxation technique and the difference in the velocity profiles sho„„ 

in Fig. 4.1. 

Figure 4.4 presents results of the present method for the 
momentum equation coupled to the energy equation vith emission 
radiation for various mass iniection rates. The velocity profiles 

In. ^ ion. coupli— o th3t 

exhibit ctianges due to muss injecti-ort cu u .ctu. 

are expected and reasonable. Furthermore, no numerical problems 
mere experienced In obtaining the momentum solutions. This is In 
contrast „ith unstable characteristics reported for the momentum 
equation by Home and Vegas (Ref. 4.10) for their - -7. 

need by Orison (Ref. 4.12) to go to an alternate method for large 

blouing rates. The present method has been used with no 

to solve the momentum equation for a mass injection rate of (p v )„/ 

= (pv) = .50 which corresponds to Howe and Vegas 

^ . r A r^Tt been found to converge 

numerically, the velocity function f has been 

quadradically in a quasilinearisation sense as shown in Fig. 4.5. 
However, the Iteration on's, see Fig. 4.6. required a damping factor 
to insure convergence, bihewise, for some cases, the temperature 


66 
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profile is weighed with the previous guess for a new guessed 
profile. The use of a maximum of 60 points in the flow field have 
ul:; 5 t , Lhe finite difference solutions 

of both the momentum and energy equations. Radiative uncoupled 
solutions converge when each point of the input and output profiles 

rsj 

IS within one percent. Convergence on 6 was also set at one percent. 
However, the addition of radiative coupling requires more stringent 
tolerances to insure convergence to a unique solution. Figures 4.7 
and 4.8 present some results of the effect of emission radiative 
coupling on the temperature profile for different convergence 
criteria. The temperature profile was selected to study emission 
radiative coupling effects on convergence since temperature is the 
most sensitive variable to both the radiative flux term and to the 
convergence tolerances . 

Figure 4.7 presents temperature profiles for no mass injection 
that result from different convergence criteria. It was found that 
this case requires more stringent convergence tolerances than most 
mass injection cases. This case was studied in detail to determine 
the least stringent tolerances needed to insure convergence to the 
correct answer. There are three convergent tolerances of direct 
importance - the tolerances on each point in the temperature and 
velocity function profile and the tolerance on the ^ change due to a 
change in the temperature profile. The tolerances are denoted as 
the energy tolerance, E; momentum tolerance, M; and the energy- 
momentum coupling tolerance on ^6 change in Fig. 4,7 and 4.8 Several 



Temperature (®K) 
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Fig. A. 7 Temperature Profiles for No Mass Injection 

Resulting from Different Convergence Criterion 


Temperature (°K) 



T|, Transformed Shock Layer Coordinate 

Fig. 4.8 Temperature Profiles for 10% Mass Injection 

Resulting from Different Convergence Criterion 
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conclusions can be obtained by close observation of Fig. 4.7. First, 
i^ 0 itiber l.OVr- or 0.3.4 tiolerances on E, M and 6 arc sufficxent to 
jLiisure _u n\ j.i. ; 1 - ^ ~ i.lit: w.^iqui ;,u ivi i. j-ou obtained using 0.14 for all 
three tolerances and also obtained using 0.1% on E and M and 0.057o 
on 'b. Secondly, it is noted that only 3 iterations between the 
momentum and energy equation were required using the 1 . 0 % or .5% 
tolerances on E, M and 6 - Thus the radiative flux divergence term 
had been calculate.d only 3 times. It v/as natural to question whether 
the integral nature of the energy equation was satisfied with this 
number of Iterations. This leads to the third conclusion. Six 
or seven energy- momentum iterations appear necessary to insure that 
the coupling between the radiative flux divergence and the remaining 
portion of uhe energy equation is correctly computed. This is 
demonstrated in Fig. 4.7 by the case where .5% tolerances on E, M 
and % vias used again but in addition it was required that the 
program make 7 energy- momentum iterations. The agreement of the 
temperature profile for this case and the unique solution is much 
better than when only 3 iterations were used. Fourth, the energy- 
momentum coupling was found to be adequately satisfied when the 
computed temperature profile did not change 6 more than 0.17o. This 
is demonstrated by two cases. One case was run with 0.54 on E and 
Hand 0 . 2 % on ' 5 . These tolerances yield a solution significantly 
different from the unique solution. The solution satisfied the 
tolerances in 8 energy-momentum iterations. Thus the radiative 
flux divergence coupling in the energy equation was apparently 
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satisfied. However, the momentum- energy coupling was not satisfied. 

Another case was run with 0.5% tolerances on E and M and with 0.1% 

profile is in quite good agreement with 

the unique solution Indicating that these tolerances are sufficient 
to insure progcr convergence. Thus. fro» these numerical experiments 
one concludes that 0.57. tolerances on E and M. 0.17. tolerance on 
'6 is required in order to insure proper convergence. Further, a 
minimum of six or seven energy- momentum iterations appear required. 

The results presented in Fig. 4.8 indicate that even though the 
above convergence tolerances are required for some shock layer 
solutions, others may not require such stringent tolerances. For 

1 • un-fT A R all the tolerances employed yielded 

the problem stated in Fig. 4.8 aii tne 

essentially the same tusuUs. It should be noted however th.it all 
the cases were required to iterate betueen the energy and a»mentut. 
equation.s at least 5 times thus satisfying the radiative flux 
divergence and energy equation coupling. Other conditions not 
studied here or in Chapter 5 might require different criteria. 

The establishment of the necessary and sufficient conditions 
to assure convergence to a unique solution discussed above was 
carried out using the emission radiation model. These criteria 
were found to be quite satisfactory for calculation which included 
line and continuum radiation coupling. However for some flight 
conditions difficulties were experienced In converging to the 
required criterion* 
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The principle cause of convergence difficulties is the very 
strong non-linear coupling between the temperature and radiative 

i ... i icultics were experienced for 
cases of no mass injection but convergence became more difficult 
with increasing pressure for mass injection cases. The extent of 
radiation coupling is, of course, increased rapidly by increased 
shock layer pressure. At pressures near one atmosphere the temperature 
profile exhibits large oscillations over the entire flow-field and 
instabilities principally near the stagnation point if the profile 
is not consti'ained in some manner. This behavior is similar to the 
behavior discussed by Anfimov and Shari (Ref. 4.13) for the same 
flight conditions. Several methods were tried to assure and to 
speed convergcncp - A copibination of methods was necessary to 
achieve satisfactory convergen.ee- per formance . This combination of 
methods will be discussed first followed by the reasoning leading to 
the use of the individual parts. 

The convergence logic of the VISRAD 3 computer program is 
presented in Fig. 4.9 in block diagram form. The oval ended blocks 
in this figure denote operations associated with convergence. The 
first three of four parts of the convergence schem.e are simple 
weighting factors. The "6, flux divergence profile at each point and 
temperature profile at each point is weighted with their respective 
computed values of the previous iteration to provide a guess for the 
current iteration. These three variables generally exhibit oscillatory 
behavior if unwciglied, thus the vjeighing procedure tends to 
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dsrnpcn tli0 oscillcitions • Th6 fouirth psirt of the convoirj^oncG scIigitig 
is a constraint on the guess of temperature profile. The guess is 
kept within 8 percent of the temperatu.. e profile from the previous 
iteration. Furthermore, the guessed temperature profile is used to 
compute updated density and viscosity profiles to be used in the next 
inomentuin solution. This procedure vjorks satisfactorily if the initial 
guess provided to start the program is close to the final result. 

This can be achieved by running cases at constant starting with 

a low pressure (i.e., 0.1 atm.) and using the converged results as 
a guess for the next pressure level. 

The first part of the convergence scheme weighs the new and old 
values of "6 in an attempt to reduce the momentum energy coupling 
betvjeen iterations. The effect of this weighing was found to be 
only slightly influential but beneficial in convergence. The 
second part of the convergence scheme vjeighs the new and old 
radiative flux divergence at each point in the flow- field. This has 
the result of dampening the oscillations in the flux divergence 
profile with the main effect being in the stagnation point region 
where the flux divergence changes sign. The third part of the 
convergence scheme weighs the new and old temperature at each point 
in the flow-field. This v;as done in order to dampen the oscillation 
in the temperature profile, improve the properties for the next 
momentum solution and to improve the flux divergence for the next 
energy solution. This part of the convergence scheme is the 
principle contributer in obtaining a converged solution. The specific 
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wsighing values for the preceding variables were determined by 
numerical experimentation and appear to be satisfactory; however, 
due to the nonlinear coupling involves uiJ the Urunge in coupling 
for different flight conditions the weighing factors are a compromise. 

A converged solution was obtained using the first three parts of 
the convergence scheme. The behavior from iteration to iteration of 
stagnation point temperature and radiative flux divergence are 
shown in Figs. A. 10b, 4.10a and 4.11 respectively. The stagnation 
point values of temperature and flux divergence were chosen for 
presentation since they are generally the last value to converge 
in the profile. Both ~ and temperature exhibit large oscillations 
during the first five iterations. Corresponding to the large 
temperature oscillations, Fig. 4.11 show.- a rrach larger change in 
tlie flux divergence. Even though convergence was achieved for this 
low pressure the results indicated that the temperature had to be 
constrained to achieve convergence for higher values of pressure. 

This realization lead to the constraining procedure of part four of 
the convergence scheme. By not allowing the guessed temperature to 
deviate from the previous profile at each T] location by more than a 
fixed percentage the large oscillations were reduced. Fig. 4.12 
shows a typical convergence procedure for a pressure of 0.5 atm. 
using a maximum percentage change of 4.04. The temperature profiles 
in Fig. 4.12 are for intermediate iterations and the converged 
profile is not shown (it is reported in Chapter 5), however, the 
progress toward convergence is shown. The eleventh and twelvth 



179 



Fig. A . 10a Temperature Convergence Behavior 



Fig. 4.10b Energy Coupled 6 Convergence Behavior 
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iterations are near convergence in the near wall and near shock 
region. The dip in the temper profiles near the stagnation 
point is characteristically the last irregularity to disappear 
before convergence. The results shown in Fig. 4.12 indicate that 
perhaps the 4.0 percent change constraint was somev/hat smaller than 
required. Subsequent calculations indicated a value of 8.0% is 
sufficient and provides a good compromise between stability and 
time required for computation. 

From the results presented in Figs. 4.1 and 4.3 it can be 
concluded that apparently the second order effects do not signifi- 
cantly contribute to the description of the flov7-field characteristics 
and therefore are justifiably neglected. Furthermore, the present 
method agrees with ether described methods for the case w^here these 
methods are applicable. The numerical solution of the momentum 
equation is quite rapid and exhibits no numerical difficulties- 
Numerical difficulties were experienced in the solution of the 
energy equation for mass injection cases where the radiative flux 
divergence profile changes from a large negative to large positive 
number in the stagnation point region. The convergence scheme 
presented in Fig. 4.9 is shown to satisfactorily overcome these 


difficulties . 
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CllAFl^ER 5 


STAGNATTO'^ TJNF RESHT.TS 

Two of the basic objectives of this research, as stated in 
Chapter 1, were to develop a computer program solution for the 
stagnation line flovj field equations and to perform param.etric studies 
for hyperbolic entry conditions. The results presented herein 
illustrates that these two objectives have been realized. 

Two computer programs have been developed which provide two 
levels of detail in defining stagnation line characteristics. A 
coupled set of first order stagnation line equations are solved 
numerically by implicit finite differences in the program called 
VISRAD 3 which provides the most complete analysis. The equations 
and nun^erical procedures used are given in Chapter 4. A second 
program called RADCOR was used to compute radiative heating rates 
for no ablation using the radiative cooling parameter correlation 
discussed in Chapter 3 (Eq. 3.12 and 3.13). 

The two computer programs have been used to perform parametric 
studies and the results of these studies are presented herein. The 
parametric studies were conducted for continuum flight conditions 
corresponding to hyperbolic earth entry velocities. These conditions 
produce radiative heating rates which are much larger than convective 
heating (Fig. 1.7). Flight velocities below 36000 ft/sec were not 
considered since the Apollo flight data and other re-entry data are 
available and define the heating and material response behavior for 
these less severe conditions. The conditions considered are: 
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Free stream velocity: 36000 to 58000 ft/sec 

Post shock pressure: .01 to 2.0 atm. 

Fig. 2.7 illustrates that these coad it j on;,, ■■ o. in the range of 

applicability of the shock layer equations and are appropriate for 
hyperbolic entry trajectories which lie in the domain of interest. 

In addition to flight conditions, additional parameters must be 
specified to obtain a shock layer solution at one point in a trajectory 

The basic specifying variables are: 


Variables 


U or T \ Free- stream 
“ ^ ■ or 


} or P 

00 O 

R 


Post shock conditions 
Body Radius 


1 


(de/d5)^_Q- Initial shock curvature | 


(pv) 
'•H 'v? 

T 

w 

C. 

iw 


- Mass injection rate 

- Surface Temperature | 

- Elemental mass fraction 
at the wall 


Specified by 


Trajectory 


Vehicle Shape 

Assumed (zero for a 
concentric shock) 


Ablator Response Parameters 


Throughout the results presented the »ell temperature and elemental 

species composition at the wall have been specified as: 

T = 3450°K (sublimation temperature 

of phenolic nylon at ^ 

P = 1.0 atm) 

C. = .7303 carbon j 

.0729 hydrogen 
.0496 nitrogen 1 

.1472 oxygen I 


407o nylon 

607o phenolic resin 


quasi- steady 
state ablation 
assumed 
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The primary emphasis in the parametric studies was to quantitatively 
determine the effects of ablation products and radiative energy 
transport on surface heating races, iruiii this xniorination coupled 
ablator- shock layer solutions were obtained. In addition, results 
of the calculations v;ere studied to provide additional understanding 
of shock layer processes such that the importance of the various 
transport and coupling processes could be assessed. 

In describing radiation, ablation and other effects the word 
"coupled*' has been extensively used in this chapter. It might be 
well to reflect on its connotation as used here before proceeding. 
Mathematically, the stagnation line equations are solved as a set 
of ordinary differential equations in one dependent variable 
each. Each of the equations contains variable coefficients and 
terms which are functions of and/or include the dependent variables 
from the other equations. Thus the equations are coupled in that a 
solution to the set of equations can not be obtained unless itera- 
tions between equations on the dependent variables are made to satisfy 
the set simultaneously. The extent of coupling may be described by 
a brief discussion of weak coupling. If one equation is weakly coupled 
to another equation large variations in the dependent variable of 
one equation produces only small variations in the dependent variable 
of the weakly coupled equation. The integrodif ferent ial nature 
of the energy equation produces another type of coupling. The 
integro term in the energy equation has been treated as a separate 
equation with respect to the remainder of the energy equation. 
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Iterations are performed to satisfy both parts simultaneously. The 
integro term is the radiative flux divergence which is quite 
obviously dependent on temperature and xs thus coupled to the differ- 
ential portion of the energy equation. The phrase "ablator- shock 
layer coupling" has a slightly different connotation. It is used 
to describe the ablator and shock layer conditions which yield a 
coropatable set of surface boundary conditions. 

The first four sections of this chapter present parametric 
studies of the stagnation line flow- field obtained using the VISRAD 3 
program. Some results from the RADCOR program are also presented 
for comparative purposes. In addition to the parametric studies the 
fifth section of this chapter presents heating rates obtained from 
a radiative cooling parameter correlation used in the RADCOR computer 
program. The main results in this section are a set of graphs for 
heating rates which were obtained using the RADCOR computer program. 
These graphs permit hand calculations of hyperbolic entry radiative 
heating rates for no ablation, A method is suggested for estimating 
the effect of ablation on heating rates. The ablation adjusted rates 
could then be used to obtain quasi- steady state ablation rates. 
EFFECTS OF R\DIATION COUPLING ON THE SHOCK STAIHI-OFF DISTANCE 

The shock wave location is a boundary condition for the thin 
shock layer equations. Mathematically this boundary condition is 
known as a free boundary and is determined by the solution of the 
equations in the bounded region. For the problem under investigation, 
the shock stand-off distance is determined by the radiation- gas 
dynamic coupling which occurs in the shock layer. 
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The purpose of this section is to present results which detnon- 
3 trate the radiation- ga. dynamo coupling offsets on the stand-off 
di 3 ta„cs. Those results were omnlned using the V.«« 3 and RABCO. 
a„.putor program uhleh are discussed In Appends D and . respectively. 

A aystemtlc study Is presented for the stand-off distance for Ulght 

velocities between 36000 and 58000 ft./sec. assuming T„ = 3450 K, ^ 

I oV ThP influence of the concentric 
(pv) = 0.0 and a concentric shock. The intiue 

shod assumption is assessed at a typical flight condition. 

addition, the location of the stagnation point for mass Injection 

cases is eKamlned for a number of flight conditions. 

Stagnation line solutions were obtained for free-strean, velocities 

between 36000 and 58000 ft./sec. at post shock pressure levels of 

„ , line and continuum 

p = 1 . 0 , .50, .10 and .01 acmospheieio us^nt, --- 

radiation model. The nondimens lonal stand-off drstanoes resulting 

from these calculations are shown in Fig. 5.1 as a function of 

cVinrV Alonr with the computed results are 
density ratio across the shock. Along wic 

plots of two commonly used correlation equations which were 
developed to predict nonradlatlvely coupled stand-off distances. It 
ia noted that the correlation of Inouye (Ref. 5.1) predicts l.arger 
atand-off distances than any obtained by the present method; however, 
the correlation presented in Ref. 5.2 agrees to within 3 % with the 
four points for 36000 ft./sec. computed with the present flov^field 

aolutlon. It will be shown that the effects of radiation coupling 

1 S shows the maximum effect), 

is small for this flight velocity (Fig. 5.19 shows 

c Rr.f S 2 nrovidcs a quite reasonable estimate 

Thus the correlation from Ref. 5.2 pro 4 

of the nonradlatively coupled stand-off distance. 



Lonal Stand-Off Distance) 


.06 



Fig. 5.1 


Radiative 

Distances 

Pressure 


Coupled Shock Stand-off 
for Various Post Shock 
Levels 
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The present results shown In Fig. 5.1 indicate some additional 

interesting processes. Most obviously, the stand-off distance is a 

^ . .1 u - -iistant pressure 

strong function of the post pi v 

lines show a double value for 6 at constant p. This effect, which 
occurs at the higher velocities, is a natural result of air therco- 
dynan.ics being used to solve the Ranhine-Hugonlot equations. As 
discussed in Chapter 3 the Rankine-Hugonlot equations are solved 
using the air tnodel of Hansen (Ref. 5.17) which assur.es an ideal 
dissociating and ionizing gas. It is also noted that the turn in the 
constant pressure line gets sharper with decreasing pressure until 
the stand-off distance becomes double valued for the pressure 
p = .01 atm. 

A measure of the radiative- gas dynamic coupling effects 
atand-off distance can be shown by observing the actual stand-off 
distance. 6, divided by the radiationless stand-off distance. A. The 
nondimensional radiationless stand-off distance. A, selected for the 
present work is represented by the equation 


A = p/(l +V^) 


(5.1) 


obtained from Ref. 5.2 This equation was used since it appears to 
be compatible with the flow- field predictions for small radiation 
coupling as shown in Fig. 5.1 and the use of Eq. 5.1 makes the infor- 
mation to be presented more accessible to other investigators than 
radiationless flow-field results. Accordingly the deviation from the 
radiationless shock layer distance is expressed as 



o 


o - 

a - 
A - 



r (Radiative Cooling Parameter) 


Fig. 5.2 Radiative Cooling Effects on the 
Stand-off Distance for Constant 
Post Shock Pressure Flight 


6/A = 6(1 +yIsp/3)/p 


(5.2) 


Figure 5.2 shows the effects : ? i I on the stand-off 

distance as a function of the radiative cooling parameter, F. The 
cooling parameter was computed using the same line and continuum 
radiation model used in the flow-field solution and is based on the 
radiationless stand-off distance, A. The results indicates that the 
ratio 6/A varies linearly with F for small values of F . It is noted 
that there is an apparent translation at 6/A = 1.0 for different 
pressure levels. This is attributed to the representation of A with 
Eq. 5.2 rather than using the actual radiationless stand-off distance 
A comparison of the present results with previous estimates of 
Goulard (Ref. 5.3) shown in Fig. 5.2 provides additional insight 
into the radiation coupling effects. Goulard formulated the stag- 
nation line problem with a simple inviscid flow field model and a 
gray gas radiation model. The effects of radiation - gas dynamic 
coupling on the stand-off distance were studied using a perturbation 
scheme for small values of F. Results were obtained for optically 
thin and thick limits. The optically thick result, Eq. 63 of Ref. 


5.1, 


6/A = 1 - % F E2(: 


■T1. 


(5.3) 


is in terms of a parameter, 1)^ = ^6^6^ * parameter is the produ 

of the post shock density, frequency averaged mass absorption 
coefficient and A. The optically thin results were identical to 
the preceding expression with the exception that the exponential 
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integral, E^, »as equal to one. For the optically thick Italt, 
Goulard found a mini»u>» value of iM for horlrontal flight and this 


minimum corresponded to a vai 


uu "il* 


,a)d iiii. optically thin 


results are shown in Fig. 5.2 The optically thin li„it predicts a 
larger reduction of the stand-off distance than that observed for 
P = .01 atm; and this result is in agreement with Goulard's assess- 
m!„t that his optically thin relationship would overestimate the 
cooling effect. It should be emphasised also that Goulard's analysis 
is valid only for F « 1- T>>'= optically thick line tends to under- 

predict the radiative cooling effects for small values of F. It is 
interesting to note, however, that the trends shown, by the simplified 
model are in agreement with the present results. That is the 
increase in optical thickness corresponding to increasing post shock 
pressure tends to decrease the radiative cooling effect of reducing 
the shock stand-off distance. 

The functional relationship of the radiative heat transfer 
coefficient and 6/A is shown in Fig. 5.3 for the same conditions as 
those in Fig. 5.2. These results indicate the nonlinear relationship 
between the stand-off distance and the radiative heating. This 
contrast to adiabatic inviscid analyses which show a direct 
proportionality relation. 

The results presented in Fig. 5.2 und 5.3 are for constant post 
shock pressures which is approximately horizontal flight (see lig. 2.7). 
The results presented in Figs. 5.4 and 5.5 are for vertical flight fer 
a fixed sot of velocities. The optically thin and thick estimates 



Vertical Entry 
R = 9 ft 




r (Radiative Cooling Parameter) 


g. 5.4 


Radiative Cooling Effects on 
Stand-off Distance for Vertical 
Entry Flight 
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„£ Goulard appear to be la agrocaant whore applicable with the 
pPesont reauUs for the U„ - 36000 ft./see. curve. The aet of curves 
appears to torn a surface which hu„ au .a .- otic Unit on the optically 

I- t-ci r> the shaded line in Fig. 5.4) 
thick side. Further, the asymtote (r.c. the 

1 r Ipcic than .5 and corresponds 

appears nearly linear for values of F less 

to the P > 1.0 and 0.5 at., line in Fig. 5.2. For larger values of 
r a n.inl!u» value in the asyctote for 6/6 appears to he approached. 

It is anticipated that for flight velocities higher than 58000 ft./ 
sec. precursor radiation will become significant resulting in an 
increase in 6/6. since a portion of the energy lost by radiation 
through the shock wave which reduces the stand-off dxstance will b 
absorbed by the on coming gas and returned to the shock layer. It 
„ight be pointed out that this minimum was understandably not 
predicted by the analysis of Goulard. This minimum appears to occur 
for large values of F while Goulard's analysis was limited by the 

assumption that F « 

Figure 5.5 points out features not apparent tn Fig. 5.6. The 
actual nondimensional stand-off distance and 6/6 ate plotted against 
post shock pressure for an intermediate velocity of 50000 ft./sec. 
lu Fig. 5.5. The stand-off distance for this vertical flight case 
exhibits a maximum whereas 6/6 monotomlcally decreases with 

increasing post shock pressure. 

u. • f- in t-his section of the results presented are for 

To this point in tnis seccxu 

a single body radius, k - 9 ft. Figure 5.6 presents the functional 
dependence of the stand-off distance on F for various values of body 
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radius. Even though 6 verses V is linear, plots of 6 verses the 
corresponding R or would not show a linear dependence. Also, 

since Fig. 5.6 is for one flight cond.: ; : : 5 ^/A exhibit the 

same dependencies. It is interesting to observe that the smaller 
nose radius body will produce the largest nondimens ional stand-off 
distance. This, of course, is a result of the snialler radiative 
loss experienced by the smaller shock layer. It should be pointed 
out that although the results given in Fig. 5.6 show a linear 
relationship, the flight conditions lie along the asjnntote shown in 
Fig. 5.4. Consequently, for flight conditions other than those 
along the asymtote one might well expect a nonlinear relation between 


6 and F. 

Duririg the development of the stagnation equations in Chapter 2, 
the shock bluntness at the stagnation line was demonstrated to be 
unknown. The bluntness parameter, (de/d§)^^Q, has been assumed 
zero corresponding to a concentric shock in most of the present 
work. To evaluate the effect of this assumption on the stand-off 
distance and radiative heating rate, (de/d§)^^Q was left as a para- 
meter in the momentum equation (Eq. 2.153) where dcp/d§ = 1. - de/d§. 
Parametric variation of (de/d§)^^Q resulted in the effects shown in 
Fig. 5.7. The results indicate that the shock stand-off distance 
and radiative heat transfer coefficient both increase in a near 
proportionate manner with (de/d^)„_g. Further, the effect of 
(de/d§)^_g is approximately the same for the two radiation models 
used. The work of Burns and Oliver (Ref. 5.4) indicated that 
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«, /dO " for a hemspharlcal shaped body at the sa« flight 

.oldltl!:! used to obtain Fig. 5.7. Consequently one uould expect a 
c to to,, one, cate In stand-off distance and heating rate coefficient 

cr *- TKt c? tdGX C cut 3.CG ITlciy y of COUlTSG y 

because of non- concentric effects. This perce 

chat with flight conditions and body site. Additional 
change somewhat witn ri-xgui. 

♦-•esm linp shock bluntness are 
observations related to the stagnation line 

given in Chnpier 6- 

Xhe shoch stand-off distance Is. In addition to the previously 
etated variables, detercined by the amount of ablation products 
being injected Into the shoch layer. Kass injection essentially 

existing on the «all side of the stagnation point and air species 

TO- ^ R and 5 9 show the stagnation point 

*.u -r'- c-r'-io Fieures 5.o ana j.j 

on tne sUk>cu c 

location as a function of rmtss injection rate and free stream 
velocity respectively. The nondimensional stagnation point location, 

_ . is shown as a function of (pv)„ for three post shoch 

° pi„ 5 8. The results indicate that the stagnation 

pressure levels rn Fig. 

point location is essentially independent of shock layer pressure 
level. Also shovm in Fig. 5.8 are results from Ref. 5.5 at (pv)„ - 
.20 for a slightly different body radius and carbon phenolic ablation 

products . 

The stagnation point location, for a given post shock pressu 
level, as a function of free stream velocity is given in Fig. 5.9. 
Present results are compared to those of Ref. 5.5 which are for 
different pressure levels shown. The present results and those of 


(Stagnation Point Location) 
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Ref. 5.5 indicate that for the given mass injection rate the ablation 
product layer thickness increases on a percentage basis as the free- 
str.aa. .etoci.y rs increased. Moreover, the ablation layer thickness 
reported by Ref. 5.5 is approximately 6.0 to 7.07« of the shock 
layer smaller than the present results. This agreement appears quite 
reasonable in the light that Ref. 5.5 results are for a different 
ablator and slightly larger body radius. The results presented in 
Figs. 5.8 and 5.9 may be of benefit for investigators using inviscid 
analysis approximating the inner and outer regions of the shock layer. 

heating rates and ablatoe. coupling 

The main concern in designing an ablative heat protection system 
is the prediction of the rate and total amount of energy which 
would be transferred to the ablator during atmospheric entry. For 
hyperbolic entry conditions the main mode of energy transfer to the 
ablator surface is by radiation as illustrated in Fig. 1.7 for no 
ablation. The effect of mass injection due to ablation is to reduce 
the convective heating to a negligible level compared to radiative 
heating. Accordingly, the results presented in this section concern 
the quantitative definition of radiative heating levels for typical 
hyperbolic entry flight conditions, A systematic variation of problem 
defining variables was made to establish cause and effect relation- 
ships. Coupled ablator and shock layer solutions were obtained for 

flight conditions typical of vertical entry. 

Many of the radiative heating rates results presented in this 
section and others are in terms of the nondimensional radiative heat 
transfer coefficient which has been defined as 
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(5.4) 


C = ^ 3 

^ % Pc 


.His soefiicleni represents t.o traction ot the flux o£ £lo„ energy 
entering tHe sKocU iayer per nnit area nhic. is transierre. to the 
aurface by radiation, the energy not transferred to the surface y 
nadiation or convection is radiativeXy lost through the shoch or 
stained by the gases surrounding the vehicle". These gases flon 
over the vehicle and equilibrate in its »ahe vith the atmosphere. 

Radiative heating rate results, for no mass injection, are 
presented in Tig. 3.10 for conditions nhich through traditional use 
have become a standard for comparing radiation calculation methods. 

c ,t. rms^nt-d .11 contain line and continuum r.odiataon and 

The results 

agree to oithin ± 107. of the average. The comparison made is 
essentially between computational methods since as pointed out by 
Wilson (Ref. 5.6) the basic data for the radiation calculation rs 
common in the worR compared. It is noted that tends to an 

asymtotic limit at the higher free stream velocities. 

This asyantotic behavior is exhibited to a lesser eKtent as the 
post shocR pressure is lowered for the .same flight velocities as 
shown in Pig. 3.11. Figure 3.11 presents results from the radratrve 
cooling parameter calculations and shock layer calcu 

Shock la.er calculation results are in general slightly larger than 
the correlation results except at the lowest pressure. 


* Neglecting precursor heating effects. 



(Radiative Heat Transfer Coefficient) 


1.0 atm 


O 



U , Free Stream Velocity (ft/sec) x 10 

CO 


Fig. 5.10 Methods Comparison of Computed 

Radiative Heat Transfer Coefficients 
for No Mass Injection 



(Radiative Heat Transfer Coefficient) 
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demonstrate the nonnegligible effect of the shock layer pressure 
level on radiative heating. 

ihe effect of body size on radiative heating without ablation 
is demonstrated in Fig. 5.12 for one set of flight conditions. All 
of the results presented contain line and continuum radiation 
coupling. The results are markedly different than the linear body 
radius dependency predicted by adiabatic shock layer calculations. 

In fact the present results show a near linear dependence on (-Cn R) . 
The results of Ref. 5.7 also appear linearly dependent on {in R) 
for values of R greater than 1.5. Calculations at additional flight 
conditions are necessary to determine if the Xn(R) relation is 
characteristic of coupled shock layer solutions which include both 
line and continuum radiation. 

The role of line and continuum radiation coupling is important 
throughout the peak heating period of hyperbolic entry. To illustrate 
this point, stagnation line heating rates were computed for a 10-g 
undershoot boundary trajectory with an initial entry velocity of 
55,000 fps. The undershoot trajectory, presented in Fig. 11 of Ref. 
5.8, was used by Engel and Spradley (Ref. 5.9) to assess the role of 
radiative absorption effects on heating loads. This earlier work did 
not include the effects of line radiation. Fig. 5.13 presents a 
comparison between the present results and those of Ref. 5.9. All 
of the results presented account for radiative losses and thus point 
out the importance of using a radiation model which includes line 
and continuum radiation without optical depth restrictions. From 



u = 50000 ft/sec 
00 

P. = .583 atm 
0 

(pv) = 0 

w 

A Rigdon et. al. (Ref. 5.5) 

O Page et. al. (Ref. 5.7) 

• Coinputed Points 



R ADC OR 
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Fig. 5.13 it it observed that both the emission model and the optically 
thin continuum model of Ref. 5.9 overprcdict the radiative heating 
rates co,;ii..teo with the RADCOR and VISRAD 3 programs. Further, it 
is noted tint the results of Rot. 5.9, using continuum radiation 
only, significantly underpredicts the line and continuum heating rates. 
The present results obtained using the cidiative cooling parameter agree 
quite well with the results from the shock layer solution. It is 
interesting to note that the present results and those of Ref. 5.9 
show a net absorption in approxinntely 97. of the shock layer near 
the wall. Thus the line and continuum and continuum only results 
Indicate that including line radiation does not significantly change 

this observation. 

- .. nhipctives of the current research is to 

One or tuti ~ ^ 

evaluate the the effects of ablation product mass injection and to 
obtain ablator shock layer coupled solutions. It was found that 
previous work. Ref. 5.5, 5.6, 5.10 and others, obtained heating 
rates for specified ablation rates from shock layer solutions which 
exhibited no apparent correlation as demonstrated in Ref. 5.5. Most 
of the results of the previous work was for various flight velocities, 
altitudes, body radii and ablation rates of carbon phenolic. 
Unfortunately either too few results were obtained in one analysis 
and/or more than one specifying parameter was changed from case to 
case. This has lead to an Incomplete understanding of the effective- 
ness of ablation products in reducing the radiative heating to the 
body. Specification of the altitude as an independent parameter 
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exe^puries the dltUceUics encountered. I£. £or exa^U. - were 
to hotd the altitude consfint as uell as the body site, wall te^era- 
ture ana ..ass inlectlon rate and vary the llisht velocity the follow- 
ing would be observed. As shown in Fig. 2.7 the post shock temperature 
end pressure levels would both change. This is in addition to a 
Change in post shock velocity as computed using the Kankine-Hugoniot 
eguations. Thus b, holding altitude (or free stream density) 
constant and changing flight velocity, a change in three shock layer 

parameters occurs rather than one as desired. Since the shock 

•ciQvaKIp pffect on radiative heating 
layer pressure level has a considerable ettect 

as shown in Fig. 5.11. it is desirable to specify this parameter 

rather than altitude. Accordingly, for a constant a change rn 

_ only one thermodynamic parameter, T^, 

free stream ‘ c. ^ j 

T it is not possible to change 

and the post shoch velocity. In general it rs n 

night conditions without altering at least two shock layer parameters. 
Consequently the present parametric studies were conducted at 
specified pressure levels for various free stream velocities . 

Radiative heating rate results for - 50000 ft/sec, R = 9 ft. 

= 3450“k are shown In Figs. 5.14 and 5.15 for different levels 
of polt shock pressnre and ablation rate. The results in Fig. 5.14 
illustrate the post shock pressure effects on heating at constant 
(pv)„ are somewhat similar In shape but are not similar to the no 
„ass”lnJectlon line. The results show chat an ablation rate (pv)„ 

= .05 reduces the radiative heating rate substantially and addrtronal 
ablation has a much smaller percentage effect. This effect is also 
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shown in Ft.. 5.15 which presents another projection of the C^, 

P . (pw) enrface. The constant corves e.hihit a characteristic 
“ but do not overlap when nondta,ensionallzed by 

the heatin. rate Efficient for no «.ss injection. (C,^)„. 

Host investigators present heating rate results for nass 
injection in a plot of C /(C )„ verses (pv)„ irrespective of the 

free streas, velocity, post shock pressure or body radius. As pornted 

*- , 1 s' eaCT-rppment between results, 
out above this had led to apparent disagreemen 

the (ov) parameter is not sufficient 
The obvious conclusion is that the P 

to correlate the results. Fig. 5.16 was prepared to. in part, point 

out why (pv)„ is insufficient. The results show a ..uite obvious 

1 - cthnck ores sure and exhibit a minimum m the 
dependence on post shock pressu 

constant (pv)^ curves near = .3 atm. Two values of nondime 
heating for carbon phenolic are also shown in Fig. 5.16. These 
values, although not at the precise flight conditions and body srte 

of the present results, were thought to be near enough for comparative 

c TT-;ior»ri 5.6) shows ^ sUghtly 

Thp = .10 case of Wilson o. } 

purposes. Inc 

higher rate of heating than extrapolation of the present results 
would yield and the (pv)„ = case of Rigdon et. al. (Ref. 5.7) 
is slightly lower than present results. The agreement is quite 

reasonable since the results are for different ablato 

■ t-c: indicate that the same kind of post shock 
importantly, the two points indicate tna 

t From these observations it is realized 

pressure dependence is shown. From these 

that plotting C /(C )„ versus (pv)„ irrespective of pressure level 

„y load to errof.eous\oncluslons regarding the effectiveness of 
ablation in reducing radiative heating. 
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Fig. 5.16 Ablation Rate Effectiveness in 
Reducing Radiant Heating as a 
Function of Post Shock Pressure 
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To examine the effect of free stream velocity on the radiative 

heating at a constant ablation rate Fig. 5.17 was prepared. As shown 

lalf of the figure the effectiveness of the ablation 

products for the conditions examined exhibits a minimum (i.e., 

C /(C ) exhibits a maximum) near 50000 ft/sec. The reason for the 

\ \ 0 

maximum in 0 upper half of 

Fig. 5.17. Even though the actual heating rate increases with velocity 
for (pv) = .20 its rate of increase is smaller than for no ablation. 

In the lower half of Fig. 5.17 values of given for 

different post shock pressure levels at = 50000 ft/sec to 
indicate how the pressure level would shift similar curves verses U^. 

In addition, three computed points from Rigdon, et. al. (Ref. 5.5) 
for carbon phorolic are shown. These results for a constant altitude 
are in reasonable agreement if one considers the post shock pressure 
level shifts. 

Of prime importance in the design of an ablative heat shield xs 

the accurate determination of the amount of material which will be 

lost during entry. To compute the amount of material lost by ablation 

processes requires a coupled solution of the flow field and ablator. 

For quasi- steady ablation the coupling is adequately described by the 

surface balance equations presented in Chapter 3. Fig. 5.18 presents 

a graphical solution matching the flow field and ablator response for 

five different post shock pressure levels. The point of crossing of 

the constant P shock layer lines and the linear segments for the 
6 

given Pg of the ablator response as read from Fig. 3. A yields the 
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coupled solution. Several observations can be »ade fro„ n,. 5.18. 
- „ , , u ,.t ,,yer coupling becomes important at approalfatoly 


P . .10 a». cor the Clight velocity and body sire speclCled. 
p!rther the ,uasl- steady ablator response assumption would not be 
valid below Pj = .10 atm. since the heating is insuCClcient 
the surface at the sublimation temperature unless the char surface 
la reteved by other than the sublimation mechanism such as oridation. 
Also for coupled ablation rates at - . 10 and lower the species 

boundary condition of the third kind at the surface may become 

j PcpVi fRef 5 11). However to obtain coupled 
important as indicated by Esch (Ref. b.ix;. 

ablator- shock layer solutions for = .10 and larger the surfa 
balance equations used appear quite adequate. In addition, it is 
observed that as the post shock pressure is increased the intersection 
engle of the lines from the ablation analysis and the shock Uyer 
analysis becomes increasingly acute. Thus errors which are inherent 
in the heating analysis are reflected in a greater uncertainty in 

the ablation rate at the higher pressures. 

in sus»ary, the results presented provides the most complete 

description of ahlator- shock layer coupling which has been reported. 

radiative r.OlIPLED SHOCK LA YER CHARACTERISTIC S 

Examining characteristic cf radiative coupled shock layers not 
only leads to a better understanding of the processes which occur but 
also permits the assessment of the relative importance of the various 
processes. Accordingly, this section is devoted to observations of 
shock layer characteristics such as temperature and velocity profiles 
under a variety of different flight conditions. 
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The ofCects o£ Increased radiative coupling is shoun in ligs. 

. no Nondi^enaional shocK layer te^erature profiles are 


• Fie 5 19 for a constant post shock pressure and six free 
shown xn ixg. 

i-nri fh-it as the free stream velocity is 
cairar-it-ips It IS noted that as 
stream velociuieb. 

lacreased. the „ondls,enslonal temperature profile is louered as a 

result Of increased radiative loss (i.e. cooling). In general the 

• of-oiv linear between (y/6) 
are approximately 
temperature profiles are niip 

_ f- T,ifh the constant temperature 
and (y/6) - .8. This is in contrast with 

profiles hypothesised hy Conlard (Ref. 9.3) for optically t re , 
shock layers. Xhe regions of rapid temperature change near the «a 

aud the shock, anticipated hy (Ref. 5.3).are exhihited hy the 

-TVs. efforts of radiation coupling has 
current shock layer solutions. The effects 

„ r--- oa the coraentun. transfer In the shook l.ayei as 

a much smaller elreou on aic 

shoun in ng. 5.Z0. Tuo ohservations ...ay he made houever. Prrst. 
the slope of the f' profile at the wall is decreased as the free 

stream velocity and corresponding radiative cooling is rncreased. 

*- • Aoi ^n one concludes that 
Since the shear at the wall is proportrcal 

effect of radiative coupling is to reduce the wall shear stress, 
is in agreen.nt with previous ohservations discussed hy Anderson 
(Ref 3.12). secondly, radiative cooling effects produce corresponding 

temperature and thus density changes in the shook layer. The density 

. * from the wall* 

in a nonlinear variation m f aw y 
changes result m a nono-i-utt 

stated another way the radiative cooling effect makes the whole 
shock layer exhibit viscous hehavior. Again this concurs with 
puhUshed work of Hoshirakii and Wilson (Ref. 3.13) as well as others. 
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Temperature profiles from shock layer solutions using the 

tion model are given in Figs. 5.21, 5.22 and 5.23. 

A comparison of temperature profiles obtained using the emission 
model and LRAD 3 model is given in Fig. 5.21. The results show that 
the emission model yields tlie proper temperature solution behavior 
though it tends to overpredict the radiative coolxng effect. 


even 


T!,3 use of the emission model to isolate cause and effects is 
desirable since the model is only dependent on shock layer temperature 
and pressure\ With the use of the emission radiation model, the 
effects of including the density variations due to ablation products 
was examined. The results shown in Fig. 5.22 are for axr C^, k and 
p. properties. The effect of including ablation products in the 
density calculation, rather than using an air value, is to shift the 
temperature profile and stagnation point. Further, the stand 
distance is substantially increased. The effect was considered 
significant enough that it has been included in all of the ablation 
coupled results of this work. Rigdon, et. al. (Ref. 5.5) did not 
isolate this effect and used air values in most of their calculations. 

A comparison of the temperature profiles computed usxng a bxnary 
solution of the species equation (Ref. 5.11) and usxng the constant 
elemental two zone model of the present work is given in Fig. 

The effects on the temperature profile by including a binary speexes 
solution appear negligible. It should be noted that the radiation 


* It is noted that the model was developed as a correlation of air 
radiation characteristics only* 
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coupling duo to opocics variation, ia not accounted for since the 

K 1 . V ^ ! ^'his comparison. The results 

indicate, however, that the densit, variation resulting fro„ diffusion 

is adequately approximated. 

The effects of radiation coupling changes on the teeperature 
profile Which arises fro., the two zone species .odel approximation 
la Illustrated in Fig. 5.2A. It is noted that only a slight change 
in the temperature profile is observed. Further, the percentage 
change in the radiative heating rate for the two solutions was 
observed to be 4.0%. In addition Each (Ref. 5.11) showed that multi- 
component diffusion effects tended to reduce this difference. A 
eomparlson rs made in Fig. 5.25 of elemental expositions for typical 
shock layer conditions showing the approximation made in using the 
two zone constant elemental model, the multicomponent solution, 
which was an uncoupled analysis, of Ref. 5.11 agrees more closely 
with the present model than the binary solution". From the results 

presented, it Is concluded that the two zone constant elemental 

.4 nt-o for use in an equilibrium solution of 
model appears quite adequate for use in a q 

the shock layer equations. It is realized, hewever. that if finite rate 
chemistry effects of the ablation species are found to be pronounced, 
this species approximation would need revision. 

An additional discussion regarding properties is in order at 
this point. The present results are for air mixture values of C^, k 

" It is existed that a selection_of ^^ 

Rrf!“!ri"ould'r'educrthe‘i?fferenc^ in the binary and r„lti- 

component concentrations. 
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I 1 1 o'.ror The results and 

and ^ which were used throughout the shock lay . 

hv .t. «l. (Ref. 5.5) for carbon 

.henolrc ahlatloo products and those ftesented by Bsch (Ref. 3.U) 

for Phenolic nylon ablation products indicate the use of air stature 

values for k and p throughout the shock layer is sufficiently a^ua 

pc For the cases examined by Ref. 5. an 
for calculation purposes. For tne 

Kef 3.11 the xnclusion of ablation product chan.es in the mrxture 

thermal conductivity and viscosity chan.ed the heatin. no more than 

h-Pd affect the shock layer profiles 
^r^t- These properties artecc 
1 0 to 2.0 percent. ihese p p 

nl), in a s.all rosion near the stasnation point enplainins therr 

1 t-- -n Results presented by Esch 

ra t-he overall solution. Resuicb p 
small influence on the over a i 

f,,. 3.11) Which include ablation species in the heat capacity s o 
a Significant decrease in the tenper.ure Profile near the wall, fro. 
those shown in Fig. 2.23- The reacting heat capacity f 

nylon ablation products was reported to be as Urge as 17 tit«s 

('Rpf 5 11) The radiative energy a sor e 
greater than that of air (Ref. 5.U). 

• .In c rransferred into internal modes 
by the ablation products is thus 

rather than translational energy .odes as predicted using air C 

Even though the teeperature was decreased significantly 

a,PP C the radiative heat transfer was 
wall using the ablation pro p 

Ef„ 3.2.. Other cases reported showed even smaller percentage 
ehanges in heating. Consequently the use of air heat capaertres rn 
me present work does not appear to have introduced significant 

j f. raffprts on the heat capacity, 
ablation product effects on 
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,-Ur. of the ablation rate 

It is interesting to examine the effects ot tne 

and post shock pressure on the ' o,-k . tespernture. norm! 

velocity and tangential velocity. Fids. 3.26. 3.27. 5.28 present 

tes,>eraturc profiles for different ablation rates at three post 

. , - Increasing the ablation rate decreases the tempera- 

shock pressures, increasinb 

cure near the »all as expected. Thus the convective heating Is 
reduced to a negligible level. The temperature profiles exhibit 
plateaus near the stagnation point and in the ablation layer. The 
plateaus appear characteristic of the ablation injection cases and 
were aslo observed by Rigdon ct. al. (Ref. 5.5) and Esch (5.11). 

The plateau near the stagnation point occurs near the peak carbon 
atom concentration which has a primary role in absorption of radiation. 
The plateau closer to the wall is a result of molecular absorption 
of radiant energy which gives rise to the increase in temperature 
over nonabsorbing results such as those exhibited by the emission 
model in Fig. 5.23. These plateaus appear more accentuated by 
increasing shock layer pressure levels. Increasing pressure 
naturally increases the radiative coupling effects and thus gives 
rise to the accentuation. 

Shock layer pressure levels have a much smaller effect on the 
velocity field and thus only one pressure level is presented herein.^ 
The normal veltcity. v/ and the tangential velocity function, f ' . 
are given in Fig. 5.29 for four ablation rates as a function of y/6. 
The effect of mass injection is to change the character of the f 
profile while the same character is maintained by the normal velocity. 
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Fig. 5.29 Effects of Ablation Rates on 
f' and v/(U ) Profiles at 


.50 atm 


v/ ( U ) (Normal Velocity) 
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^ ..f the f' profile is changed from one like a boundary 

The character of the t proLxi.o 

t * 4-u rrtocc: iniection exhibiting 

= 0.0 to those vjith mass 

layer profile for vp 

; j. h. stagnation point, 

linear f' variations in regions uu ca. . 

is noted near and on both 

The region of significant vrscous effects 

,f.es of the stagnation point; this point is indicated by vertical 

arrows on the f' profiles. The results in Fig. 5.29 also show a 

in V/ profiles between the wall and the stagnation point. 

u- nf the density decrease in this region. 

This behavior is a result of the density 

BADIATIVE transport CII^ ACTERIST^ 

one of the main processes studied in this work is that of 
radiative energy transfer within the shock layer. The results in 
this section illustrates the Important radiative characteristres 
observed which determine the radiative heating rates to the body's 

surface. 

The effects of phonolic nylon ablation products on the line and 

oontinuum radiative flun to the surface are shown in Tabs b.l and 

f mo flux for the twelve spectral intervals 
5.2. The continuum surface flux t 

• -r n 1 Tn each spectral interval the 
considered is shown in Tab. 5.1. In eacn si 

„«,„it„de of the flur is given for tbree different ablation rate 
cases. By comparing tin magnitude of the fluK for the three ablation 
rates in each spectral interval, one observes that the ablation 
pcodnets block most of continuum flua above hv - 8. ev. further, the 
continuum flur below 5 ev is not absorbed but slightly enhanced by 
increased ablation rates. Tab. 5.2 presents the line flur counterpart 
of Tab. 5.1 where the line flur is located in the spectrum at nine 
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• o* At the two lowest frequency line centers 
line center frequencies . At the two 

the nux i. changed only slightly by ablation product injection 

„hereas at larger ireguencies the eiioe e »..c pronounced. The 

j ^ ^ n c nb*^ervGcl th6 linG 
greatest flua reduction by ablation products 

renter hu = 10.4 ev. Thus the nylon phenolic ablation products are 
the least effective in absorbing energy in the frequency levels 
below 5 ev. Unfortunately a significant percentage of the radiant 
energy Is shown to be transferred in these lower frequency levels 

(i.e. 797o for (pv)^ = . 10) . 

The solution of the energy equation is directly coupled to the 
radiative transfer through the radiative flux divergence 
this equation. The radiative flua divergence profiles for aero and 

, • TTn- rr S ^0 It Is notcd that the 

207o ablation rates are shown in Fig. 5.30. 

effects of ablation products is to reduce both the wall and shock 
values of the flux divergence. However, the negative area, denoting 
net absorption, shown in Fig. 5.30 is substantially increased by 

injection of ablation products. The sharp dip in the profile is 

ma ?n7 ablation. This sharp change in the 
evident for both zero and 204 ablation. 

• H i f f iculties , if a small step size is not 
profile presents numerical difficuitie , 

used locally, as discussed in Chapter 4. The swall peah near (y/6) ■ 
.gg occurs near the ^wi™» carbon atos. concentration and represents 

s„all net emission primarily attributed to this species. 

The influence of ablation rate on the total line and total 


1 for line flux for some line groups is a result 
rrunfa^rorrtion of continuum flux a, noted by Ref. 5.6. 
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contlnuu. nux at two locations in the shocU layer for two post 
shock presaore levels are shown in Figs. 5.31 and 5.32. These 
£i,„res present the line and contrnuu.., itua to the surface at (y/5) = 

0.0 and the fluxes toward and away from the surface at stagnati 

point, (y/6) V (y/6)„=o- "" 

tlux toward the surface at v = 0 are Increased by increasing ablation 
rates. This is in part due to the flux from the ablation layer to 
the air layer shown as the lower curves in Figs. 5.31 and 5.32. 

The upper and lower set of curves gives an indication of the radiative 
coupling between the ablation and air layer. In a previous section 
(see Fig. 5.15) it was shown that the total radiative heating to 
the surface was decreased by increasing the ablation rate. The middle 
curves of Figs. 5.31 and 5.32 show that both line and continuum 
contributions to the total flux decrease as a function of increased 
ablation. However, for the cases shown. Increased ablation beyond 
(CV) . .05 has very little effect on the continuum part of the total 
flux! Thus regardless of the ablation rate above (pv)„= .05 the effec- 
tiveness of the ablation products in reducing the total surface 
heating appears limited by the continuum radiative processes. 

It has been demonstrated in previous sections that the post 
shock pressure is quite important in determining radiative heating 
rates. To illustrate the pressure dependence of the line and 
continuum parts of the radiative flux, Fig. 5.33 was prepared. The 
line and continuum flux toward the body at V = 0 and the resulting 
two parts which arrive at the surface are shown in the lower half of 
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Fig. 5.35- results demonstrate several Important effects. 

I„ the pressure range considered, the line and continuum flua tovard 
the body at V . 0 vary approximately Unear ly with pressure. Also 
the line flux at the wall shows a near linear dependency. The 
continuum surface flux, however, is lower than the surface line 
flux at lower pressures and increases nonlinearly to values larger 
than the surface line flux at higher pressures. This cross over 
effect with pressure is seen ,uite clearly from the results in the 
upper half of Fig. 5.33. The two curves shown are for the line 
and continuum ratios of the respective flux at the wall. (q^)„, 
divided by the respective flux at the stagnation point, 

These curves show, that for a constant ablation rate, line radiation 

through the ablation layer is attenuated less at low pressures than 

- t-viio fnr continuum radiation, 
at higher pressures and the converse x 

This effect is the reason for the minimum in nondimens ional total 
radiative heating rate curves of Fig. 5.16 as xllustrated for one 
ablation rate in Fig. 5.34. Furthermore, since the continuum flux 
is essentially unaffected by ablation rates above (pv)^ - .05 as 

indicated in Figs. 5.31 and 5.32 the continuum mode for radiant 

rtf for ablator coupled 

energy transfer is the predominate mode of heattng 

shock layers at the higher pressure levels considered. 

in Chapter 3 an effort was made to select the molecules which 
would have an effect on the radiative transport. A shock layer calcu- 
lation was made to assess the Importance of including molecules in 
the radiation calculation. The results shown in Fig. 5.35 indicate 
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•c anite irnportant which is in 
rhe role of molecular radiation is q 
that the r t, f 5 6 Excluding 

•tv. the assess,.cnt t,ade In Ret. 5.6. 
contrast with -eased tl heating by 

„Veentes t„ and stsnttteantt, 

^ : Ta^et te.etatnte .otUe. Xt ta aVao noted that 

the stagnation pom assessment of 

,„£Ue near the wall IS observed. 

tea^erature pr vade with an uncoupled analysts 

molecular radiation tn ^ 

(t.e. using a ter.peratu 

-, ^ This resulted m a sman 

tv, the ettects of .oleeular radiation were 
heating and conse,u» ,,, 

.onsidered negUgt ^ 

radiant heating may ^niroling to the tempera- 

rofile the molecular radiation couplmg 
tetnpcrature p > coupling is to 

f-ia is significant. The final effec 
ture profile i o .-le and corresponding shock 

.^reciably^nse^^eh^;;;” ,,,£.ee heating signiiieantly. 

layer prop illustrate the 

T*. r^T-<^c:pnted in Fig. 

The results pres 

^ av.^ral carbon radiation. A 5Z pe 
iv^ortance oi .to.s 

,er£ace heattng resu ^peulation for the 

but including all other specie ^ ^^^.^^Von was not included in 

ease studied. Carbon and hydrogen 

. „f Ref 5.10. Smith et. al. (Re£. 

the calculations of Ref. ablution to results 

,„.„lts for phenolic nylon ablation 
Vbeir heating rate ^ iucludes line 

for carbon phenolic of Rigdon et. a . 
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, a hyaro^en radiation for identical conditions and found 

earbon and hy e ^ ,„,atantially larger. The 

the heating rates of ( 

. , ,y i ..ii M. the present 

iHiportance of the carbon radxaixou a reported 

investisation susBosts t^ t the .alor dtf 

.esults «y uell he attrihuted to the differences .e 

an their res.ectiye ^„_d in 

Radiative transport results for 

, 36 am 3 37 . A »aai™ in the line contrihutron to 

P . 35 a». ia ehoun in the upper 

aurface heating at approximate y • 

ya- 76 This behavior rs ditterem. 

plot of Frg. 3 . 36 . 

^ .1- • the surface heating seen m Hg 

attenuated line contri utio larger at the 

fVasat- the line contribution is 1 g 
The lower plot indicates that 

, „„ Pcee stream velocities than at higher ones. It 

lower tree sui oi dered both 

-’nd velocity ranges consider 
ramtire pressure and veiucj-t-i 

that over the enti P ^ y-ficantly to 

p,ae and continnum radiation processes contribute signrfica 

rue total surface heating 

Total radiative heating 

,lysis are compared In nondl^nsional form with resu 
the present analys ,„„elates 

Of other investigators rn Tig. • • 

,nite well with the radiative cooling — ’ ^ 

•d check layer analyses (.Rets. 

the results of two invisci pl„ 5 16 lie below the 

Ml of the shoch layer results presented in g. • 

t c p^f 5 7. As pointed out by Fage 
transparent gas theory resu ts o • • ^ 

1 in Ref. 5.7 the transparent theory is not a eq 

• ce the true optical properties of the shock 
predicting heating since 
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layer are not uniformly optically thin. The Important point Illus- 
trated in Fig. 5.37 is that the radiative heating for aero ablation 
can be correlated with a single parameter. C. and mat both viscous 
and invlscld analyses yield similar results. The success of this 
correlative par.nmeter . which is the basis of the KADCOR program 
calculations, leads to speculation that perhaps a correlation which 
includes ablation product coupling might be developed. 

An attempt was made to determine if the heating rate results of 
the present work both with and without ablation could be correlated 
with shock layer parameters. Some of the difficulties encountered 
are shown in Fig. 5.38. First, the results show that the heat 
transfer coefficient for the aero ablation cases are correlated 
rather well by the cooling parameter a, empected from the results in 
Fig. 5.37. The constant ablation rate lines for 11^ - 50000 ft/.ec 
exhibit an analogous shape and relationship to the rero ablation 
line as that shown in Fig. 5.14. It was illustrated in the second 
section of this chapter that 0^ was not a simple function of ablation 
rate that could be represented by one curve. Thus an appropriate 
means of reducing the zero ablation curve and the thiee curves for 
different ablation rates at = 50000 ft/sec has not been found. 
Moreover, results for a constant pressure, P^, and ab.ation rate 
various free stream velocities shown in Fig. 5.38 quite obviously 
Indicates that varying at constant 0^ has a different effect 
heating than varying at constant P^ when both are done at the 
same ablation rate. This is different from the zero ablation cases. 
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. c: QQ t-iiG body radius may also be 

Beyond the effects shown In FiE. 5.38 F 

Thnq it is concluded that 
Important In ablation coupled problems. Thus 

. iM ..t I s in addition to 

at least two parameters, excluding rauiu. , 

1 a hr- reouired to correlate t e 
the radiative cooling parameter would bo requ 

radiative heating results presented. 

HBATIhG wa TKS FKOM C 

, nature of the shock layer at the stagnation line has 
The complex nature oj . 

been amply demonstrated in previous sections. It Is highly desirah e 

to being able to perform detailed shock layer coa.utatlons. According 
this section presents radiative heating rate results base 
radiative cooling parameter correlation which permits hand calcu- 
lations Of hyperooUc entry heating. The results presented were 

computed using the kABCOR computer program documented 

as also used to computc heating rates for a 
The RADCOR program was also c 

c, on. 1 -rate the effects which may be realized 
100% co^ atmosphere to demo 

by entry into atmospheres typical of Mars or Venus. 

Xr was shown In the previous section that the radiative heating 

rates for no ablation were correlated quite well by the radiative 

cooling parameter. Fig. 5.39 presents a comparison of results from 

Che present shock layer calculations with a correlation equation from 

blvlngston and «lllard (kef. 5.15). The shock layer results show 

only a small pressure dependence and are slightly above the 

, e- Heating rate results from using the correlation are, 

correlation. Heating 

1 -fK-n the present uncertainty limits, i.e. ± 10%, of 
certainly within the presem. 
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... ... .1 .h. ...u« 

, ,,,, , : . justified. 

heating rates for prelin^inary des.gu 

. set of fisures .hich can be used for hand calculations o 

i.npd using the RADCOR program and 
radiative heating rates were developed 

are designated as Figs. 5.40 to 5.45. Ihcse figures p 

radiative heat transfer coeffioient as a inaction oi free strea. 

£ t-n 2.0 atmospheres, 

velocity for seven post shock pressures ro 

• fnr n snecific body radius ranging from 
Each of the six figures xs for a spe 

l.o to 11.0 ieet. 

Civen a speciiic traoeotorv, defined hy p.. and tine, and 

Ki-j developed from Figs. 5.^0 
f, heating rate history can be devcjppe 
body radius, a heating 

■int- in a trajectory is determined 
. ra'-c ^or one point in a ci-cij 

to 5.A5. Tne heat^ne ra.e _o i 

r from the figure for the appropriate 
■ rr n and P and reading C from tne r, 

^ ^ ^ te can then be computed from the 

body radius. The actual heating r 

. • of C using p and U . The required values for the 

definrtron of C using co . 

, P am usually confuted from and m the 

post shock pressure, Pg^> 

u ,-iMPr if the normal shock information for 
trajectory analysis. Howt- , 

M. Fig 2 7 can be used to obtain iron, the 
is not readily available, Fig. 

This process may be repeated for selects 
free stream conditions. This p 

citing in a radiative heating rate history 
points in a trajectory resulting 

3 ,.fi.r to that Shown in Fig. 3.13. Graphical integration of the 
heating rate curve yields the total heating load at the stagnai-ion 

point neglecting ablation effects. 
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Fig. 5.40 RADCOR Radiative Heat Rate Results for a Body 
Radius = 1.0 Feet 
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Fig. 5.45 RADCOR Radiative Heating Rate Results 
for a Body Radius = 11.0 Feet 
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An estimte of ablation adjusted heating rates may be made 
using Fig. 5.16 and 5.17. If quasi-steady ablation appears appro- 
priate for the conditions under consideration, sufircient information 
is presented herein to estimate ablator coupled heating rates and 

ablation rates. From the results presented in Figs. 5.16, 5.17, 

. r Kc 1 = 7 for all values of P. and would 

and 5.18, setting ^ 6 

yield a conservative estivate of the ablator- shock layer coupling. 
Ablation adjusted heating rates could then be computed using 
C /(C - .7. An estivate of the ablation rate can then be vade 

„!"ing 3.A. Finally, the total n.ss loss can then be estivated 

by graphical integration of an ablation rate verses time plot. This 
procedure nuy be useful in defining a body site and r.anges 
flight condition lor vote detailed analysis. The radiative cooling 
correlation has been used by Livingston and Wllllard (Ref. 5.15) 
and Stickford (Ref. 5.16) to calculate stagnation point heating 
rates for atvospheres containing different percentages of CO^ and 
The radiative cooling paraveter correlation «as shown to predict 
experin«ntally veasured radiative heating rates within the scatter 
„f the data for both air and a 90% CO, - 107. N, atvospheres. The 
experivental heating rate data was obtained for a 5 inch diavater 
hevisphere and a 1.25 in. diaveter truncated cylinder in a shock tube. 
Data was taken over a free streav velocity range of 18.000 to 30.000 
£t/sec and a post shock temperature of 7.500 to 15.000 K and 
pressure range 1.0 to 7.0 atmospheres. Thus, although the free 
atream velocities are lower than those considered in the present work. 
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the post shoch conditions otc sln,ilar. The agreement between theory 
and data shown by Ref. 5.15 lends validity to the use of the cooling 
parameter correlation in the present work. Fur u.csmore. entension 
of the present work using the RADCOR program to include CO^ - 
atmospheres is justified by the agreement of these results. 

To Illustrate the difference In radiative heating rates resulting 
from different atmospheres Fig. 5.46 was prepared. The results 
shown are for air and 1007. CO^ at the same post shock conditions 
which correspond to a representative flight condition for Venus entry. 
The radiative heat transfer coefficient is given as a function of 
stand-off distance which, of course, is a function of body radius. 

The results show a much larger radiative heating for the 1007. CO^ 

luu . Further, the 

atmosphere than for air at the specn-ueo t-o - 

* '.ii.T rt-iT- T 1 V fir; a function of stand~ 

heating rate increases more rapidly for 

off distance than for 100% 00^. 

CHAPTER CLOSURE 

To recapitulate, stagnation line shock layer solutions were 
presented for various specified ablation races and for coupled 
shock layer conditions. The characteristics of radiation and ablation 
coupling in the shock layer were quantitatively shown and discussed. 
These characteristics Included the radiation and ablation effects 
on the stand-off distance, temperature profile, velocity profile 
and radiant heating. A simplified heating rate calculation method, 
based on the radiative cooling parameter, was used to develop a set 
of graphs which can be used to make hand calculation estimates of 
hyperbolic entry heating rates. 
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CHAPTER 6 


AROUND THE BODY RESULTS 

The purpose o£ this chapter is tuo-iold. First, procedures used 
to con.pute the shoch shape around the body are discussed and results 
con^aring different methods are presented. The effects of post shoch 
Ptessure distributions on the shoch shape are studied paraa.tricaUy 
demonstrating shoch location sensitivity to such variations, 
secondly, radiative heating rate distribution results from radiative 
cooling parameter calculations are presented. These results are 
compared uith more detailed analyses which require flow- field 
aolutions and at most are found to differ by about ± 61. 

SHCX:K StL\PE C.AL ClUuATIONS 

A„'y discussion of calculations of the location of the bow shoch 
wave produced by a blunt body naturally Involves a discussion of the 
post shoch pressure and surface pressure distributions. This is 

true since the location of the bow shoch is determined by the post 

. . a- o t-n a first approximation directly 

shock pressure which in turn rs to a iirst app 

determined by the surface pressure. 

There are three methods which may be used to determine the 

shock shape and pressure distribution around the body. First, we 
will consider the technique used in Refs. 6.1 and 6.2 and others. 

The shock shape Is specified a from which the wall pressure 

distribution is calculated as the solution proceeds around the body. 
An output shock shape is calculated from the geometrically relation 

(Ref. 6.3) 
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e = tan"^ [^(d6/d|)/(l + r6)j 


( 6 . 1 ) 


where 6 is calculated as a result or tiie x - momeutran solution. This 
output angle in cohered „ith the input angle. II the input and 
output are nearly the ea»e the solution is said to be converged. 

A second technique Involves specifying a wall pressure distribution 
a priori. Preferably this distribution is Vnown iron experimental 
data for hypersonic mch numbers. The change in pressure due to 
radiation coupling is Justifiably neglected (Ref. 6.1). A shock 
shape is also assumed. The shock layer equations are solved around 
the body and the calculated and input pressure distributions are 
compared. The shock shape is numerically adjusted according to the 
pressure difference. The solution is repeated until sansfactory 

IQ nbtained The third technique involved is 
pressure convergence is obtainea. 

a simultaneous solution of the geometrical relation 


6 = (1 + h6) tan ed§ + 5 q 

^ A 


( 6 . 2 ) 


with the shock layer equations. The post shock and surface pressure 
distributions are automatically calculated as part of the 
layer equations solution. Only ono around the body iteration is 

needed for this technique. 

The first two techniques have been implemented in a modified 
version of the computer program described in Ref. 6.1 and 6.2. 
Unfortunately both of these techniques have limitations. The first 
technique consusQS computer time because of the many around the 
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•r each estimate of 

• c which are necessary. In addxtro , 
body iterations which are ne 

..e SHOO. s.ape .e .V 

second technique also consumes a great uea I oi ae.ae.ter tin 

the iterations around the body that are necessary for 

^ an advantage over tne 

t'Hp second technique does hav 
convergence. The secon 

„rst Since the input and output pressures are smooth; thus updatr 
of the shoch shape can be done automatically. However, this 
technique does not necessarily satisty the geo.trical differential 
equation stated above. Since the first two techniques are su 

to undesirable Umitations.it is the purpose here to explore 

t the third technique. 

feasibility of implementing the 

^ h pressure P differs, in general, from the 

The post shock pressure, 

u ho-v in'>-'e location. The pressure 

v-rr p for the same bouy ant,**- 
surface pressure, P„, 

ond temperature variations across the shoch layer are shown r . 

, f These results, obtained from unpublished worh of Sprac. ey 

;::i .ef. Show that the post shooh pressure is smaller than 

neater than the wan pressure for larger l.ch numbers, .he resu 

a . Pig 61 were obtained using a computer solutron of 
presented in Fig. • . a 6 2 The referenced 

shock layer equations documented in Ref. 

.ethod used a modified Karman-Pohlhausen integral method to solv 

the X - momentum equation and a finite difference relaxa 

• An inviscid y - momentum 

procedure to evaluate the energy equatron. An 
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equation solution for the given x - moinentum solution is used to 
obtain the y pressure variation, Results from this program yield 
the required post shock pressure variation for the present investi- 
gation of shock shape calculations. 

The geometrical integral equation (6,2) was integrated with a 
simultaneous Rankine-Hugoniot solution to obtain the bow shock shape 
for a given post shock pressure distribution. Specifying the free- 
stream velocity and density and the post shock pressure is sufficient 
to solve for the shock angle, e, (from Eqs. 2.79 to 2,82) using a 
thermal and caloric equations of state. This shock angle is then 
used in Eq. 6.2 to solve for the local stand-off distance, 6* The 
integration of Eq. 6,2 was carried out using a simple predictor 
corrector method. Numerical experiments shewed that a step size of 
“ .5 deg. was sufficiently small to insure convergence to a 
unique solution. 

Figure 6.2 presents results of the integration and Rankine- 
Hugoniot solution. Two aspects of shock shape calculations are 
demonstrated in this figure. First, a comparison of the present 
integration method and the shock shape computed using the flow- 
field numerical solution of Ref. 6.2 is presented (i.e. PR = 0.0). 

The results are in quite good agreement as expected for the following 
reason. The post shock pressure distribution used in the shock 
shape integration was obtained from the flow- field solution. The 
flov;- field solution was obtained by specifying the wall pressure 



6(Nondimensior;al Stand-off Distance) 
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distribution and iterating on the shock shape until the wall pressure 
converged. The wall pressure was specified using 

~ “ 1.25 sin^ 0 + 0.284 sin^ 0) (6.3) 

o 

which is an inverse solution correlation of Inouye (Ref. 6,5). Thus 
one aspect demonstrated in Fig. 6.2 is that the numerical integration 
of Eq. 6.2 does yield a solution comparable to the flow-field solution. 

The second aspect of shock shape calculations shown in Fig, 6.2 
is the sensitivity of the stand-off distance to the post shock pressure 
distribution which was studied by parametrically varying this 
distribution with the equation: 


= (V ) , (1 -f PR 

6 t u-iruii-iJ 


( 0 °/ 34 °)) 


(6.4) 


In the above expression PR represents an error that might be 
experienced by the forward integration procedure in any flow- field 
solution. The resulting stand-off distances for PR = + .05 and + .01 
are shown in Fig. 6.2 These results indicate that the normal 
direction variations in pressure as those shown in Fig. 6.1 are 
quite important in determining the shock shape. Moreover, a high 
degree of accuracy must be maintained in a flow- field calculation to 
prevent computational inaccuracies from being amplified in the 
resulting shock location. 

Assuming that numerical accuracy can be maintained, results 
presented in Fig. 6.3 indicate that the present integration method 
is to be preferred over the two other methods previously used with 



Shock Angle (Deg.) 
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flow-field solutions. This figure presents the shock 

angle computed by the present method and the output shock angle from 


the flow- field solution. 


computed from 


differencing the output stand-off distance of the flow-field 
solution and using a two point difference derivative in Eq. 6.1. 
Small fluxuations in the output stand-off distance result in the 
irregular changes in e and acculumutive error. Using the present 
method a smooth shock angle is computed; this is necessary if the 

shock shape is to be updated by an iteration procedure in a computer 
program. 

The set of thin shock layer equations for around the body flow 
are parabolic. Thus, one expects the initial conditions, i.e., 
stagnation line conditions, to influence the down stream solution. 

The initial conditions are determined, in part, by the shock 
curvature at the stagnation line. As pointed out in Chapter 2 this 
stagnation boundary condition is unkno;^, and it is usually assumed. 
The true elliptic nature of the problem indicates thar this boundary 
condition is determined by downstream effects. As pointed out in 
some early work of Hoshizaki (Ref, 6.6) a downstream boundary 
condition could theoretically be substituted for the initial shock 
curvature. This downstream condition is the irtviscid shock angle 
far from the body where all disturbances from the body are negligible. 
Since it is not practical to attempt to satisfy this downstream 
boundary condition, the initial shock curvature, (de/d§)^_^, must 
be specified. 



In order to examine uncertainty in the post shock pressure 
distribution resulting from uncertainty in the shock curvature at 
the stagnation line, Fig. f ; d. The flow-field solution 

of Ref. 6.2 was used to compute the post shock pressure distribution 
for a given initial shock curvature, (de/d§)^^Q. The three distri- 
butions shoxvn were obtained using Eq. 6.3 and the surface pressure 
convergence metfiod. The maximum percent difference in pressure 
noted at 9 “ 34 is approximately -5%. As was demonstrated in Fig, 
6.2 a change of pressure of this magnitude causes a very major 
change in the shock shape. Accordingly, we may conclude that 
uncertainties in the initial shock curvature, which result from the 
elliptic nature of the problem, may preclude an accurate estimate of 
the shock shape. 

As a practical matter the shock shape must be approximated to 
compute the surface heating rates. Therefore, it is suggested that 
the initial curvature be assumed such that calculations may be made. 
The accuracy of the results should be regarded with an awareness 
of the error which may be introduced by the assumed curvature. 

To conclude, it is observed that the integration of the 
geometrical equation (Eq. 6.2) results in a smooth shock shape. This 
is important if the shock location is to be iterated upon numerically, 
Further the shock shape location is strongly coupled to both the 
initial shock curvature and pressure changes across the shock layer. 
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RADIATIVE HEATING ?IATE DISTRIBUTIO^:S 

The success of using a radiative cooling parameter correlation 
to compute stagnation lin; . rompted an 

investigation into its applicability to non- stagnation line calcu- 
lations. The goal here is to obtain a computationally rapid means 
of estimating distributed heating rates for no mass injection. The 
effects of mass injection may be accounted for by using stagnation 
line results presented in Chapter 5 or correlations v;hen available 
to determine the ablation - shock layer coupled heating rates. 

One of the basic assumptions usually made in calculations of 
shock layer radiative transfer is that the shock layer can be 
treated locally as a planar infinite slab for calculations along 
the stagnation line or around the body. This assumption is inherent 
in the radiative cooling parameter correlation presented in Chapter 
3. Thus this assumption is consistent with radiative transfer 
calculations made in flow- field analysis. In order to use the 
radiative cooling parameter correlation, developed from stagnation 
line calculations, for nonstagnation line calculations it is necessary 
to assume that the radiative transfer process and the local tempera- 
ture profile are similar to those of the stagnation line. Fig. 6.1b 
provides a qualitative basis for judging such an assumption. The 
constant temperature lines shown in Fig. 6.1b were computed using the 
flow- field analysis of Ref. 6.2 with radiative coupling of the 
emission model. The results indicate that although there is a 
change in character of the temperature along constant lines, the 
temperature level changes rather slowly as a function of distance 
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along the body. This is in contrast to the more rapid change in 
pressure as shown in Fig. 6.1a. Consequently one expects the 
similarity assumption to be ^ i , l oi i i ct. 

Given the expected limitations, the radiative cooling parameter 
was used in the computer program RADCOR (see Appendix E) to compute 
radiative heating distributions. The isothermal flux was computed 
locally using the post shock temperature and pressure across a slab 
of thickness equal to the chosen stand-off distance. In the event 
of the lack of a better estimate a concentric shock is assumed. No 
provisions have been made to account for radiative transfer blockage 
by ablation products. Radiative heating distributions computed 
with this method are presented in Figs. 6.5, 6.6 and 6.7. 

Figure 6.5 pre.scnts a comparison of heating rate distributions 
from the present method and an inviscid flow-field solution of 
Barnwell (Ref. 6.7). Barnwell used a time-dependent finite- 
difference technique to obtain numerical solutions for the problem of 
inviscid flow of radiating equilibrium air past spheres at hyperbolic 
speeds. The results of Ref. 6.7 were computed using a two step 
absorption coefficient model which included the effects of line and 
continuum air radiation. The results of the present method were 
obtained using the shock shape computed by Baruwell.- 

Fig. 6.6 presents a comparison of the radiative heating 
distributions about a sphere from the present method and the 
viscous shock layer solution reported by Chou and Blake (Ref. 6.8). 
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Chou and Blake solved the thin shock layer equations using a Blasius 
type series expansion technique. Solutions were obtained using 
three terms in this series. Radiative transport v/as computed using 
a three band continuum model for the absorption coefficient of air. 

The agreement shown in Fig. 6.6 is not as good as in Fig. 6.5; 
however, it is noted that line radiation effects were not accounted 
for in the results of Ref. 6.7. As demonstrated by Chou and Blake, 
viscous effects do not significantly alter the radiative heating 
distribution for these no ablation cases. This provides additional 
credence to simularity assumptions inherent in the present method. 

Fig. 6.7 presents results for three flight conditions using the 
present method. Each heating rate distribution case presented took 
less than 2.0 minutes of IBM 360-65 computer time. The heating rate 
distributions were computed assuming a concentric shock for a 
spherical body. The results show the same trends reported using more 
detailed models. The main trend demonstrated is chat the heating 
rate decreases, as a function of body angle, more rapidly for lower 
flight velocities than higher ones. Additional cases not presented 
indicated a weak dependence of the heating rate distribution on body 
radius which agrees with the results of Ref. 6.7, 6.8 and others. 

In order to estimate the effects of initial shock curvature on 
the heating rate distribution, equilibrium around- the- body solutions 
were obtained using the computer program documented in Ref. 6.2 
Figure 6.8 presents the results of this parametric study. The initial 
shock curvature, (de/d§)^^Q , was varied from zero, the concentric 
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Fig. 6.8 Comparison of Heating Rate Distribution 
for Different Initial Shock Curvatures 
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.hoc. assu..ption, to a valoo o£ 0.1. The heating rates, for the 
1 toot sphere considered, was significantly increased. ~ 257. higher 
.r ,o“. by assuming an inlti.nl curvaU.rc of 0.1 either than aero. 

The initial shock curvature was reported by Burns and Oliver (Ref. 

6.9) to be approximately 0.0705 for a hemispherical body at similar 
flight conditions. It should be noted that the emission radiation 
model was used in these calculations and therefore only trends, not 
a quantitative measures, are established. Nevertheless, the large 
changes shown in Rig. 0.3 indicate that one amy safely conclude that 
fhe shock curvature at the stagnation line can have a significant 

effect on the heating rate distribution. 

From the results presented in this section one can conclude 

^ r parar,eter to compute heating rate 

that the use of a ox 

distributions appears to yield satisfactory results for preliminary 
design work. This conclusion is indicated by the agreement with 
other methods shown, the small computation time required and the 
remaining uncertainty in downstream heating rates resulting from 
initial shock curvature uncertainties. Although it has not been 
studied in this section, the present method of using the cooling 
parameter might well be applicable to non-sero angle of attack 
problems near the stagnation line if a shock location estimate is 


availfjb le . 
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CHAPTER 7 


CONCLUE I 


The results In the previous chapters give a more extensive and 
detailed quantitative description of hyperbolic entry heating than 
any previously reported studies. The shock layer and ablator 
analysis Includes all of the known significant processes. These 
processes Include ablation and radiative coupling effects within the 
viscous shock layer, radiative line and continuum processes of both 
air and ablation species, local chemical equilibrium throughout the 
shock layer and quasi-steady ablator behavior. The overall analysis 
includes tl>e best available models of each process a.,d no single 
other analysis has included all of these processes. Specific 
processes were studied with respect to their contribution to the 
shock layer heating. Where several models of the same process were 
found to yield comparable results the most simple one was Incorporated 
into the overall analysis. The detailed discussion of the mathe- 
matical model used and the results obtained provide a sound basis for 
understanding many of the characteristic processes of hyperbolic entry 
heating. In addition to the stagnation line work, the radiative 
cooling parameter, previously used only for the stagnation line, was 
shown to be applicable in computing heat rate variations around the 
body. The computer programs developed are engineering tools which 
can be used to quantitatively define .aerothermal environments not 

already considered in this study. 
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CONCLUSIONS 


Based on the results presented in this work the following 
conclusions are drawn. The 'he stagnation line 

unless stated otherwise. 

Mathematical Model 

1. As demonstrated in Chapter 4, the thin viscous shock 
layer flow processes which occur on the front face of a 
bluff body during hyperbolic planetary entry are 
accurately described by the bluff body first order shock 
layer equations stated in Tab. 2.6. 

2. Numerical solution of the stagnation line momentum and 
energy equations using quasilinearization and implicit 
finite differences was found to be both satisfactory and 
reliable. In contrast with other reported methods, no 
numerical difficulties were encountered in using this 
method for the momentum equation. Thus quasilinearization 
used with implicit finite differences to obtain a numerical 
solution to the momentum equation is highly recommended. 

Shock Location 

1. The shock stand-oft distance as a function of the radiative 
cooling parameter approaches an asymtote at high free 
stream velocities. The asymtote appears to approach a 
minimum as the cooling param,eter nears a value of one. 
radiation perturbation results of Goulard show the 
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which accounts for line and continuum transport to obtain 
estimates of the radiative heating rate variations away 
from the stagnatioii . . .k. uy ■ .uiu.Kis . i i^ed in Chapter 

4. The stagnation point location as a function of mass 

injection rate is essentially independent of the shock 
layer pressure level at = 50000 ft/sec. The stagnation 
point location was observed to move slightly away from 
the body as the flight velocity is increased at a specified 

ablation rate. 

5, The nondimensional ablation rate parameter (pv)^, is 

insufficient to correlate the nondimensional heating rate, 

P //(-. ') The nondimensional heating rate changes 

Wiuh pus.- shock pressure and free stream velocity in additi 
to (dv) . The results available for comparison indicated 
that computations for carbon phenolic ablators exhibited 
the same pressure and free stream velocity dependencies 
as the present result for phenolic nylon ablators. 

6. The ablator- shock layer coupled results indicates that 
below = .10 atm at = 50000 ft/sec the surface 

heating is insufficient to maintain the surface at the 
sublimation temperature. Consequently the quasi-steady 
approximation will not be valid unless other mechanisms, 
such as oxidation, remove the surface rapidly enough to 


on 
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maintain a constant char thickness. In addition, the 

ablator coupled results showed that at larger pressures 

(i.e* P > .5 atiuj - ai a i. i'i j- att. Increasingly 

6 

sensitive to small changes in surface heating rate. 

7. Ablation products of phenolic nylon are least effective 

in absorbing radiant energy in frequency levels below 

hv = 5.0 ev. For a typical case considered (P^ = .5 atm, 

(qv) = .10, U = 50000 ft/sec) approximately 79 percent 
^00 

of the radiative flux arriving at the surface was in the 


frequency range below hv “ 5.0 ev. 

As the shock layer pressure is increased the continuum 
contribution to the surface flux is Increased and the line 
coiLtiibutior. is decreased at a constant and (pv)^. 

This change in the relative contr ibutions from the two 
radiative mechanisms is responsible for the presssure 
dependence of the nondimensional heating, 

The continuum contribution to the surface radiative 
heating is essentially unchanged by increasing the 
ablation rate above (pv) = .05. Thus any reduction in 
radiative heating rate below that for (pv)^ = *05 is 
primarily due to blockage of line radiation. 

Sufficient information is given graphically to permit hand 
computations of hyperbolic earth entry heating rates for 
no mass injection. A method of obtaining an adjusted 
estimate of the ablator coupled radiative heating is 


suggested . 
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Gas Properties 

1. Comparison with multicomponent and binary diffusion 

species equation . O!:.; : , ; es the two zone 

constant elemental approxim,ation of this work is suffi- 
ciently accurate to predict ablative coupled heating 
rates (i.e. within 4.0% of calculations using binary 
diffusion) . 

2. The use of air values for viscosity and thermal 
conductivity rather than including ablator species effects 
is justified for engineering analysis of the shock layer 
when the main concern is prediction of surface heating 
rates. The maximum heating rate percentage change 
observed v.’as 2.0 percent. 

3. The influence of ablation product species referenced to 
that of air on the reacting heat capacity and thus the 
temperature profile and radiative heating is more 
significant than the influence produced by the transport 
properties. For the case studies the maxim.um change in 
heating rate due to difference between air and ablation 
product heat capacities was 4.8 percent. 

4. Molecular absorption of radiant energy in the ablation 
layer reduces the radiative heating rate significantly 
(i.e. 52% for the case studied). 
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rf.comi-iendations 

Recommendations for improving current ablation coupled heating 

rate calculation capability are as ioiiuvn., 

1. An analysis to determine the magnitude of the initial 
shock curvature for different ablation rates is 
recommended. The stagnation line and distributed heating 
rates are influenced significantly enough by this parameter 
to warrant a better definition than presently available. 

2. Additional experimental verification of frequency dependent 
radiation data is needed to improve heating rate calculation 
reliability. For some species, C^H and C^H, a complete 
lack of data was found. In addition, carbon soot has been 
experimentally observed in tne shock layers of ablative 
models. Since carbon soot is a strong absorber and radia- 
tively active below 5 ev, mechanisms for carbon soot 
injection into the shock layer from the ablator need 
mathematical definition. 

3. Additional cases could be run with the VISRAD 3 computer 

program for flight velocities and body radii not considered 
herein to provide a larger range of calculated results for 
hand calculations. The effects of ablation products on 
heat capacity should be included in these additional 
calculations . 

Recommendations for future analyses of ablative heat protection 

sys 


terns are as follows: 
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1. Establishment by NASA of a set of flight conditions and 
body sizes for analysis by investigators would permit 
more direct compar isuui> of rcauiLs and nuuiericai methods 
than can now be achieved. This type of coordinated work 
has been achieved for orbital entry conditions by AGARD, 

2. Studies for entry into the planetary atmospheres of Mars 
and Venus could be conducted using the VISRAD 3 program. 
Only minor changes to the program to account for arbitrary 
free stream gases would be necessary to obtain this 
additional analytical capability. 

3. The VISRAD 3 program has the flexibility to be used for 
different ablator composition. A comparative analysis vjith 
this program to examine the effectiveness of ablation 
products of different ablators in reducing radiative 
heating rates would yield information complementing test 
results . 



APPENDIX A 


GENERAL CONSERVATION EQUATIONS 

CONSERVATION EQUATIONS OF A IIirLTICOMl’ONENT , RADIATING, CHEMICALLY 
REACTING n.UID 

The conservation laws for mass, ntomentuni, and energy will be 
presented for a continuum, multicomponent fluid whose internal degrees 
of freedom are in thermodynamic equilibrium. The assumption of thermo- 
dynamic equilibrium implies that no matter how small a volume of 
fluid we are interested in there are enough molecules within the 
volume to give meaningful average properties and that regardless of 
the flow velocities of interest a temperature may be ascribed to the 
fluid. This is roughly equivalent to assuming the first postulate of 
nonequi 1 5 b''" 1 uri thermodynamics, see Fitts P.ef. A.l. 

A general property balance can be made on an element of volume V 
moving with an arbitrary velocity similar to that given in Ref. A. 2. 

The property (mass, momentum, or energy) per unit volume is designated 
by 3 . The flux of a property through a control surface is denoted 
by B (property x length) / (volume x time), and the generation of a 
property within the control volume is denoted by IB (property)/ 

(volume X time). The differential form of the general property 

balance can be written in terms of the above definitions (p 31, Ref. A. 2). 

•|^ + V • P V + V • B - IB = 0 (A.l) 

ot 

if the control volume is subsequently assumed fixed in space. Thus 
for a control volume stationary in space there is a convective flov; 
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through the control volume which is identifical to the motion term 
associated with the moving control volume, This means that if the 
general balance is derived for a moving ouuiiui volume it may be 
used for a fixed control volume, with B maintaining exactly its 
same definition. This allows B to be interpreted as a diffusive 
flow. This is obvious for the moving control volume, but the common 
practice of lumping all kinds of effects into this flux term for a 
fixed control volume effectively redefines B. Therefore, the general 
balance equation is derived in the form for a moving control volume, 
but it is fully intended to be used to describe a stationary volume 
in space. 

The general property balance Eq, A.l can also be written: 

■p + V-V3 + 3V-V + V-B - IB = 0 (A. 2) 

O t 

(1) (2) (3) (4) (5) 

The meaning of these terms is, for a control volume: 

(1) the accumulation of 3, 

(2) the convective flow of p, 

(3) the dilation of the flow, i.e. the change of B when the 
fluid is compressed or expanded, 

(A) the diffusional flux, 

(5) the generation of 3- 

Using Eq. A, 2 and specifying 3> ^ write the 

conservation equations. Consider first the conservation of mass by 
specifying 3 = P (mass/volume), B = B ~ 0. 



300 


Substitution into Eq. A. 2 yields 
Continuity: 


-|£ + V-Vp + pV-V = 0 (A. 3) 

Before proceeding to the other conservation equations let us rewrite 
the general property balance equation in another form by substituting 
3 = bp into Eq. A,2. By using this substitution and noting the 
continuity equation appears as a product of b, the general property 
balance relation can be expressed as; 

p ^ + V-B - IB = 0 (A.A)"^ 

This equation will be used to evaluate the rerr.aindcr of the conservation 
equations . 

Consider now species conservation by specifying 


b = C. 
1 

B = J. 

1 

IB = cL, 
1 

where 



Substitution of the above relations into Eq, A, 4 yields 
Spec ies Continuity: 

DC. 

P ^ + v-j. - J). = 0 (A. 5 ) 


Db 

Dt 


Bb 

is the substantial derivative of b which equals — + V*Vb, 

ot 
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Let us accept the second postulate of thernKxiynanilcs of irreversible 

processes which states that If the fluid Is not too tar fros, 

O.OUS functions of 

equilibrium, fluxes and cuiii. 

the driving force. Using this postulate the mass flux vector of 
specie i. Eq. A. 5, can be written as the sum of contributing vectors. 


(Si Vi-^ (^2 Vi+ • 


(A. 6) 


where 


L = Transport Coefficient 

Im 

X = Driving Force 

and vhle subscript "D" indicates dlffuslonal transport coefficients. 

The number of necessary terms to consider can only be discussed 
in reference to a particular application. Four terms are stated below 
from Bird ^ , Ref. A. 3, for consideration. 

n ^ r • \ 

<hiVi='^i'“' I WiiUi 


k=l 


(LpjXPi = 

(hdVi = "i'” 


- D. V Xn T 
1 




p RT 
2 


(h4\>i = '"i 


-r (g) t 


— V M.M.D. . [y M. VPj 

RT L. 3- J 3-J L J J 

‘'k - .1 




— — — / . 1 1 . . . I * . * * ♦ \ o » / ri 

p RT 1 J rju J J J A- p 


(A. 8) 
(A. 9) 

(A. 10) 


where 


= Concentration in 
Y = Mole fraction (X in Ref. A. 3) 

j J 


in total no. of moles/volume (C in Ref. A. 3) 
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Gj = Gibbs' free energy 

®ij “ ^Multicomponent diffusion coefficient 

T 

- Thermal diffusion coetiicieut 

E,. A. 7 expressed the nttss diffusion vector. Since the driving 
force is of the su,,k= measure as the flux, they are called ''conjugate". 
The conjugate transport coefficients, are the largest, l.e. mass 

is diffused primarily by mass concentration gradients. Eq's. A.8. 

A.9, and A.IO represent the mass flux vector contribution from thermal 
diffusion, pressure diffusion, and forced diffusion respectively 
There are also fluxes due to inertia and viscous terms, but they are 

very small, see appendix in Pitts Ref. A.l. Electrical and magnetic 
effects can also create fluxes. 

The definition of flux as a linear function of coefficients and 
potentials and the realisation that fluxes are tensors of various 
ranhs leads one to speculate on „hat type of cross effects can exist. 
Curie's theorem states that "fluxes whose tensorlal characters differ 
by an odd Integer cannot interact In isotropic systems," Ref A 1 
This means that the mass flux tensor and the heat flux tensor which are 
both vectors are not coupled to the reaction rate tensor (a scalar), 
or the momentum flux tensor (a second order tensor) but may be coupled 
to each other. Also, It should be observed that momentum flux tensor 
either as a second order tensor or in contracted form as a scalar may 
be coupled to the reaction rate tensor. 

With the foregoing information In mind consider the conservation 
of momentum. For substitution into the general balance equation 
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b = p/p = V 

B = - ? + IP - Pj^ 

“ ■ I ’’A 

i 

Using Eq. A,4 for momentuni conservation yields 
Momentum; 

p - V • (T - IP + Pj^) - V p.i. = 0 (A. 11) 

i 

Note that in the above equation the radiative pressure tensor, is 

included for completeness. This term is negligible for practically 
all non-nuclear problems (Ref. A. 5). 

Let us now apply the general balance equation to conservation of 

energy by specifying 

2 _ _ 

b = Q + ^ ^ p^V • = E (energy/mass) 

i 

B = q^^ (energy/volume) (length/time) 

— B = V • V • Ct - IP + Pj^) • V + ^ I 

- Sp (energy/volume - time) 

= generation by radiation + pressure tensors + external forces 
+ heat sources internal to the C.V.; i.e. induction heating, 
resistance heating, etc. 

Substitution of the above into Eq. A. 4 yields the total internal 
energy form of the energy equation 

T\-p = «=> 

P ^ - IP + Pr)*V 

+ i. • - Sp . 0 

1 


(A. 12) 



where 


a = diffusional heat -,-t;or 

4d 

0 = radiative heat flux vector 

g.*J. = heat generated in ^he system by a gravitational field 

1 


Y 

L 


Let us investigate further the diffusional and radiative heat 
flux vectors. Again imposing restrictions from thermodynamics of 
irreversible processes, the diffusional heat flux vector may be 
written as a sum of vectors 


Id - + 


(A. 13) 


where 

4i^i = I Vi - 

i 

= energy transport due to the Dufour effect 

The term is the conjugate term for this flux vector. It should 

be noted however that the right hand side definition is an arbitrary 
one. The Dufour effect arises due as the conjugate of the Soret 
effect in mass diffusion. Additional cross effects from other 
coefficients and potentials will not be considered. 

Radiative transfer of heat is propogated in an entirely different 
manner than diffusional heat transfer. Duffusional heat transfer 
mechanism depends on gradients in the gas, such as temperature, species, 
pressure or external forces as pointed out by Planck, Ref. A, 4. 
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Radiative transfer of heat is in itself entirely independent of these 
gradients in the medium through which it passes. In general, radiation 
i.b a LcU tuuie complicated phenomenon than diffusional heat transfer. 

The reason for this is that the state of the radiation at a given 
instant and at a given point of the gas can not be represented by a 
single vector as the diffusional mechanisms can. All radiative 
energy rays v;hich at a given time pass through the same point in a gas 
are independent of each other. Therefore, to specify completely the 
state of the radiation at a point the radiation intensity must be 
known in all directions which pass through the point under consideration 
Special attention will now be given to the development of the 
radiative flux and flux divergence terms which are needed in the 


Lluatiou of energy conservation. Starting with the basic concepts 


of radiative transfer in an absorbing and emitting medium. Ref. A. 4 


and A. 5, a definition of the spectral radient energy density is 


developed . 

Let f(v,r, Q^,t)dv dQ be the number of photons in the frequency 
interval v to v + dv, contained at time t in the volume element d¥ 
located about the point r, and having a direction of motion within an 
element of solid angle dLQ about the unit vector The function f 

is called the distribution function. For this definition to be 
meaningful the linear dimensions of the volume element must be 
larger than the largest wavelength C/v . 

Each photon possesses an energy hv . Therefore, the , spec t£a_l 
radiant encr gy density may be defined as the radiant energy of frequency 
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V included in a unit spectral interval and contained in a unit volume 
at the point r and at the time t by: 


lU^ (r, t) -- hv 


f dii 


4tt 


In a like mnner, the spectra 
First recall each photon possesses 
c. Therefore the quantity 


(A. 14) 

1 radiation intensity can be defined, 
a velocity equal to that of light 


hvc f(v,r, dv dO, 


(A. 15) 


represents the radiant energy in the spectral interval dv passing 
through a unit area in a unrt time in the direction within the solid 
angle dH about . The area is located at r and is normal to 
This statement is not necessarily obvious. In order to clearly rndrcate 
how and what area is located at point 7 let us follow the derivation 


of the spectral radiant energy density given by Planck Ref. A. 4. 

Consider an infinitely small clement of volume dif, at the point 
which has an arbitrary shape Fig. A.l. In order to allow for all 
rays to pass through the volume cW, we can construct about any point 
7 of dV a sphere of radius a . The radius a is selected to be large 
compared with the linear dimensions of dV but still so small that no 
appreciable absorption or scattoring of radiation occurs in the sphere 
Each ray which reaches d^ must then originate from some point on the 
surface of the sphere. 

Let us determine the amount of energy contained in d¥ which 

element of surface area da. The surface area is 


originated from an 
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chosen such that its linear dimensions are small compared to those of 

<JV . Consider the cone of rays which start at a particular point on 

a. ... . -Uiic ^.onsists of an infinite number 

of conical elements with a common vertex at a point on da each cutting 

out of the volume d¥ a certain element of length s. The solid angle 

2 

of such a conical element is .AA/a where AA denotes the area of cross 
section normal to the axis of the cone at a distance a from the 
vertex Fig. A, 1. 

In order to find the energy radiated through an element of area 
let us first define hvcf 

(r,Q^,t) - hvcf (A. 16) 

which is caJ led the spectral radiatio n intensity . Using Eq. A. 15 
and A. 16 the monochromat ic energy v;hich has passed through da and is 
in d ¥ is: • 

I dQ(s/c)da = hvcf dCl (s/c) da . (A. 17) 

V 

2 . 

where dCl - — and s is the path length in d ¥. The energy which 
a 

enters the conical element in d ¥ spreads out into a volume AAs. 

Summing up over all conical elements which originate in da and enter 
d ¥ yields 

— ^ / AAs = — ~d¥ =— dQd¥ (A. 18) 

c 2 u. c 2 c 

a a 

The symbol, dQ, used in Eqs. A. 17 and A. 18 has two different meanings; this 
difference is seldom noted in the literature. This represents the entire 
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monochromatic radiant energy contained in volume d^ resulting from 

radiation through the element of area da. To determine the total 

uioiiOi-hi uiiici I i' ..a-iLiiuL eii i^y ,..;aiiained in cfV we must integrate over 

all elements of area da contained in the surface of the sphere. For 

the procedure of this integration observe Fig. A. 2. In this case the 

incretnent in solid angle dfi = — 2 which corresponds to a cone with a 
_ (J 

vertex at r . Integrating the right hand side of Eq.*s A, 18 yields 
the total energy: 


The monochroinatic radiant energy density is obtained by dividing by 
d ^ . 

tU =- Ft di'2 (A. 19) 

V C J V 


Since the radius a does not appear in Eq^ A. 19 we can think of as 
the intensity of radiation at the point r itself or the intensity of 
radiation passing thru a unit area at r in the direction . This 
clarifies a difficult concept which is avoided in many derivations. 

From the definition of I it follows that the radiation heat 

V 

flux is a vector of magnitude 

qR(r,t) = J c lU^dv = J J (A. 20) 

in the direction of photon propagation. Let the normal to any 
surface thru point r be called n . Therefore the magnitude of the 
heat flux passing thru a unit surface area normal to n from photon 
propagation in the direction is: 
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Nonabsorbing 

Nonscattering 

region* 


^Radiation in minus 0]^ direction _to 
from the unit area in the minus 


the unit area equals the radiation 
direction. 


Fig. A. 2 Geometric Relations 


for Calculation o 


f Radiation to d¥ 
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(n-q^) - qj^(r,n,t) = J j (n.Qj^) (r,Q ,t) dO dv (A. 21) 

4rr 

ho • , ’iol ' 'f - ^ " ir can be written 

_ _ <4n _ 

"J j I (A. 22) 

v=0 Q=0 ^ ^ 

Therefore is defined at any point r^^ and time t in space. 

For the use of the radiative heat flux vector in the energy 
equation, it is desirable to be able to calculate a component of q^ 
in any coordinate direction of an orthogonal coordinate system and 
to calculate V-q^^ . These calculations may be accomplished in a more 
expeditious fashion by first vjriting the equation of radiative 
transfer . 

The radiative transfer equation states that the rate of radiative 
energy accumulated in a volurae element plus the rate that it flows 
thru the element equals the rate of generation within the element. 

The generation of radiative energy is accomplished by emission and 
absorption. The general property balance can be used by defining 


3 = 1. 


B = 0 


IB = c 1 J (1 + 

L v 


1 


2hv 


o I ) - K I 
J V vJ 


where 


J1 ^ ~ Radiative emission coefficient 
K = Radiative absorption coefficient 
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Substituting into Eq* A. 2 


'Vr 

dt ' 'l' "‘v ' v" 


=[. 


(1 + 


2hv 


3 


and noting that 


V 



- 0 


we can OTite 





J 

V 


(1 + 


2 

c 

2hv 


I ) « K I 

V 


V 


(A. 23) 


(A. 24) 


which is identical to the expression given by Zel*dovich and Raizer 
Ref. A. 5, In order to simplify Eq. A. 24 the following observations are 
made. The emission term J[ ^ can be expressed 


J1 


V 


a B 

V V 


(A. 25) 


by using Kirchoff's law and assuming local thermodynamic equilibrium* 
Note that the effective volumetric absorption coefficient 


= Kj^l - exp(-hv/k^T) J 


(A. 26) 


is the product of the absorption coefficient and the induced emission 
term. Therefore the emission terra J) ^ has both spontaneous and 
induced emission taken into account. The spontaneous emission term is 


the Planck function. 
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(A. 27) 


2hv ■— — i — -r r 

B„ = — r exp(+hv/k^l) - 1 


Uslns deItnitio„!i Eq. A. 24 can 

1 + n, • vi^ - 


be rcvnritten as; 


(A. 28) 


, A t8 la -.ItipUeO by « and intCEratad 

I£ the radiative transtcr Eq. . obtained 

the conservation of radiation equations 
over all directions the con-er 


\ 

9t 


where 






(A. 29) 


U 


vp 

Let us assume 


!J!lv , i f ^ dn = 0 

c J 


St 


St 


1 n Fn A 26 for the radiative 
alnce c is very large. Then we .ay solve Eq. A.2 

flux divergence. 

to ,4 tt _ \ (A. 30) 

- 1 _ f „ ('‘aB - J ^v^'> 

’'■V S " -^o ^ '' 

.he contribution ol the radiative ilu. divergence ter. in the energy 

oquatlon has Lportant .athcatical ra.illcations. It should be 

noticed that the ilua divergence ter. is evaluated by ” t 

• thP energy equation are differ 
n cinace The other terms m the energy q 
over all space. therefore 

1 , ..llv The radiative flua divergence term 

calculated locall> • 

' ' vrial-dlfferential equation. 

the energy equation an integro-parti 


makes 
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eojgKRV^ 

section « hcve a vector formulation of the 
, 3 Sic conservatton e^uatio.n. to, a teacttns, raiiating. conducting 
fluid. Moat flou problems are represented by the conservation 
e,uations in orthogonal coordinates. In this section the basic 
conservation laws will be written in general curvilinear orthogonal 
coordinates. This permits one to select a useful coordinate system for 
a particular problem and thus determine the appropriate coordinate 
atretching functions. Substitution of the stretching functions into 
the conservation equations in curvilinear orthogonal coordrnate 
„iu yield the appropriate governing equations for the problem of 

interest . 

a 1 --t physical tensor operations for 

Table A.l presents a oot c- ^ > 

orthogonal coordinate systems. By using the information in this 
fable we are able to write the conservation egnations in curvilinear 
orthogonal coordinates. The statement of these conations have been 
made in part by Bach. Tslen, Brodhey. Refs: A.6, A.7, and A.2 

respectively, and others. 

The steady state conservation equations can be stated as follows: 
Global Continuity: 


(A. 31) 


VpV = 0 


Species Continuity: 


y(p.V) + VJ. = a:. 


(A. 32) 
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(Table A.l continued, 
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r 






(Table A.l 
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Home n turn: 


= 0 

i. X 


(A. 33) 


Energy: 


p(VV)E + V*(qp + V I ‘ 


+ P„)-V - Sp = 0 
K 


The energy equation may also be written as follows. 

p(V-V)H + V-(qj^ + - V-C? + Pj^)-V - Sp - 0 (A. 35) 

This form is presented by Ref. A. 3 without the radiation or internal 
heat source tori:’. 

For tho purpose of writing the eonsorvation equations in curvi- 
linear orthogonal coordinates, the coordinates are chosen to be 5 ^. 

5j. and 53 corresponding to and of Tab. A. 1 respectively. 

The differential elements of length in the respective coordinate 
directions are h^dj^. h^dC^. and 1.3053 such that a differential arc 
length can be expressed as 

(dC)^ = + h^^dSj)^ + hj^dSj)^ 

where h, , h„, and h„ are called the "stretching functions" in the 
respective coordinate directions. In the following equations u, v, 
and w are the velocity components of V in the direction of increasing 

§1> and §3. 

Applying the V operator from Tab. A.l the global continuity 


equation becomes 
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. 3(h„h pu) oCh^h^pv) ^ ^ 

1 . . ^ + -17? ^ -J 

IT L ^^2 3 


(A. 37) 


h> -.ho 

In a 


..e .CC..S c„n..„uu, can 


be 


written 










[ 






S§. 


5§, 


(A. 38) 


u). = 0 
1 


In order to 


and J. F 

X ^ b 2 

flux vector 


evalnacc the a.oce e, cation the co»ponenta 

, j test be specified. The ™ss 

f the case flea vector ^ 


for a wide range 


,f fluid problen.s is well repree-.. 


by two terms 


J. 




^(D) 

The expressions 


for there two mass 


flux vectors are 


n SG . 

^ “iVij ^ 
t-i 


n^ 


- 0>) _ _L_ 
Ji - piRT 


V Y, 


] 


k= 


(A. 39) 


(A. 7) 


1 (T) ^ V Jin T 

J. 0- 


The trans 


formed components 


of the above equations are 


(A. 8) 
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Tot L L''j ><3Y 


,0» . 


1 


k=l ■^X??j,k 
k?^j 


(D) 


J. 


^ > 5 9 P 


1,§. 


2 n f^Y -. 

.j;i_ y „,„,o,,[y, yc^)„, (A.«) 

p RT A 1 j ij L .] dY, T,-V,Yj h„ 0 ^ 2 '^ 


j=l 


2 n 


1-1 2 ^2' 
k— 1 j »k 

k?fj 


dG, 


1 SY, -, 
1 k ! 


p RT L ^i^j^ij L *^3 ^^3 

j=L k=l ^ UiX 

kT^j 


,(T) ^ " ^i 5(-en T) 

^i,§l 5 


^(T) _ ~ T) 

^ > §2 ^2 ^ ^2 

T 

(T) ^ dUn T) 

i-h ' S ® 53 


(A. 41) 


For substitution into the species continuity equation 


(D) (T) 

J. _ = J. c + F 

(D) (T) 

J. - = J. r + J, F 

i>§2 ^>32 ^>§2 


J 


i,§ 


3 


(D) (T) 

. + J , p 


(A.42) 


This completes the necessary operations to explicitly write the 
species continuity equation in general orthogonal coordinates* 
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Before writing the momentum and energy equations in general 
orthogonal coordimtes the radiation pressure tensor and external 
torce rieia ue.u. .re ex. pp.d/ The resulting vector form of the two 

equations are 
Momentum: 


p(V • V)V - V .(t - IP) = 0 


(A. 43 ) 


Energy: 


p(V . V)H + V • (q^ + qj^) - V . (t) . V 0 


(A. 44 ) 


If the need to account for the additional effects should arise, the 
appropriate terms could be added to the governing equations in an 
analogous manner to the terms which will be considered. 

Using the definitions in Tab. A.l, the momentum equation can be 

written in the three orthogonal directions. 

E, - momentum: 

^1 


u 5u ^v + ^ — 

^ ^ ^^2 ‘>3 ^^3 '’ l ‘’2 ^^3 

1 1 SP 


dh 


1 ^ uw 


dh. 


2 dh„ 

V ^ 


2 dh- 
w 3 


V _ W j ^ 

dC^ ' h^h3 p d§^ 


(A. 45 ) 


Lh 


1*^2^3 


(■ 


,d(h2h3Tii) ^^^lVl3^ ^ 


■^12 ^ '^13 

+ + 


T22 5h2 


T,o Sh- 
33 J = 0 


hih2 d§2 ^ ^ 1^3 ^^3 ' ^1^2 ^^1 ^1^3 

These terms are not usually significant for gas dynamic problems. 
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F - momentum: 
^2 


Sh„ 

vw I 


'hi "B§i ' h2 a§2 ^"3 ""^3 

2 9h 


Bh„ 

\'.V 2 _l_ 

i»l'h 2 S§i h,h, o§ 


2 3 


2 Bh_ 
w 3 


_u 1 _w 3 i 1_ ||_ 

hih2 a§2 P ^2 


, r 1 ,3(t'2l'3''l2> , ^ 3(hi>>2V') 

■ 7 Lh^ 852 


_^’jj_^^-'23 ®”2 3hi 


“^33 1 _ 


^ 1^2 ^ 2‘"3 ^^3 ' ^ 1^2 ^^2 ^ 2'"3 °^2 


! = 0 


(A. 46) 


F - momentum; 

'=’3 


3L. ^ + 

hi o§i 


^ 3h. 

V 'dw Sw ^ ^ _ ■ _3 4 . 


ah_ 

wv 3 




h h ^ 

"13 -n 


SI 2 


2 

u 

ahi 

2 ah„ 
w £ + 

1 _1 

_ IL- 

^1^3 


*^2^3 °^3 

P h 

3 °^3 

1 r 

1 

^'^ 2 ^ 3 ’’^ 13^ 

a(hih 3 'i 
— 

" p Lh 

1 ^ 2 ^: 

} ^ °-l 


o §2 

'^31 

ah 3 

^ ^23 ^*^3 

■-U 

^ 1^3 

S51 

+ h3h3 a?2 

*^1^3 ^^3 


ah. 


(A. 47) 


in the above equations, the subscripts 1. 2, and 3 in the symmetric 
stress tensor denote the coordinate directions and §3 

respectively. In order to evaluate the three momentum equations the 
components of the viscous stress tensor must be defined. For a Stokes- 
fluid the stress tensor is defined by. Ref. A, 2, in terms of 
rate of strain tensor e . 


T = f(e) 
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t form for this equation in three dimensions is (Ref. A. 2) 

(A. 48) 


The simples 

j - Ai + Be + Ce ■ e 


For a Newtonian fluid 

A = - (-|i - p)V • V, B = + C = 0 

The stress tensor may now be written as 
^ = XIV • V + lie 

The components of the stress tensor are 
Tfi = XV • V + pe^j^ 

^22 ~ ^ ^^22 
T33 = XV • V + pe 33 

”^12 " '^21 " ^®12 

■^13 " ”^31 " ^^13 

■^23 " '^32 " ^^23 
Which may be written 


pSCh^h^u) aCh^h^v) a(h^h^w)-, 
1 r-r- d >r ■* 1 


A 

■^11 “ hj^h^h 

3 

B§1 

’ S^2 

3^3 J 




Sh, 

5h- 

+ 2\i 

“1 

Su , 

V 1 ^ 

« 

Lh^ 

3§1 

*^1^2 °^2 

^' 3*^1 °^3 


(A. 49) 


(A. 50) 


(A. 51a) 


L (A.51b) 


(A. 52) 
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r-9(hJuu) a(h.h v) 9(h h w)-, 
_ X ' 2 J j l_ L_± — I 

’’^22 ^ °^2 ^^3 ^ 







9hrt -j 


S§2 ^ 

\vf 


u 

hj^h, 


X r 

a(h2h2'^ 

0 a(h^h3’ 

4- 

a(h^h2V7) 





S§2 






dh 


ah„-, 

ri 

Lh^ 

5w 

S§3 ■ 

u 

>’l'’3 


V 

h^h 

3 1 

3 ^“2'^ 

■ "21 ' “[■ 

*^2 a 
h, a§, 


^2 

a 

S§2 

<^>] 


5 ^ _L. 

"^13 " ^31 ■' o§3 d§^ 


"■^3 _d_ 

■^23 = '^32 ■ d§2 ‘b S 


>3 ^2 


(A. 53) 


(A.54) 


(A. 55) 


(A. 56) 


(A. 57) 


With the preceding definition of the stress tensor, the tnomentuni 
equations become three equations expressed in the three components of 
the velocity vector. 

The total enthalpy form of the energy equation Eq. A. 44 may be 
written in general orthogonal coordinates by noting the form of the 
three operators expressed in Tab. A.l. 
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u BH j V 

BH w BH 


-1 

r^(^2*^3qD, 

.hi B§i h2 

5§2 ^3 °^3-^ 


hih^hj 

L 5§i 


i 4 b , j 

. i 

i 

pSOi^h^qj^ P 

3§2 

S§3 

J 

hih2h 

3 ^ 

2^ 


.1 

4. 

,-B {h2h2(TiiU+T2iV+T3iW)} 










+ 


bT:^ ^ -I 


(A. 58) 


The components of the shear stress have been defined in the discus- 
sion of the momentum equation. Therefore only the components of the 
heat flux vectors are left to be defined to provide a complete state- 
ment of the energy equation. 

The heat flux vector as handled previously will be described 

as the sum of the diffusional and radiative heat flux vectors. The 

diffusional heat flux vector can be expressed as a function of the 

mass flux vector by simple manipulation of the equation given by 

T 

= - k'vi + > h.J. - Nk T ) V (— ) (A. 59a) 

cZ-.N,m. N 

i i ^ ^ 

where is not the ordinary thermal conductivity coefficient. The 

usual form of the diffusional heat flux vector is written in terms 

of diffusion velocities or mass flux vectors. This form eliminates 

N. 

V(^~) from the preceding equation and adds a term to k yielding the 
ordinary thermal conductivity. This step also introduces the binary 
diffusion coefficient into the Dufour effect term. Following 
Hirschfeldcr et. al., Ref. A. 8, and substituting for the diffusion 


velocities yields: 
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kVT f I h. J . 
i 


N 


I X tt: ( 7 ; 

i 


ij 


D./ 


(A. 59b) 


# 


where D. . is the binary diffusion coefficient 


3(M. + M.)P 
_ 1 1 

1 .] i-j 


©. . = 


(A. 60) 


The diffusional heat flux vector contains terms which respectively 
represent conductive energy flux, diffusional energy flux, and 
diffusion- thermo (Uufour) energy flux. The Dufour effect is kept in 
the above equation to be consistent with keeping the Soret effect 
in the species continuity equation. At this point it is appropriate 
to point out that the thermal conductivity in the conductive flux term 
is in general a tensor. For the case of an isentropic fluid, the 
conductivity reduces to a scalar. This is the form used in the 

diffusional energy flux vector above. 

Having stated the vector form of the diffusional heat flux vector, 
the components needed in the energy equation can be expressed. 

k 5T , T . T 

*^D,1 ' h^ u i 


T T 

N_. d/ f 


N 


L L m D,. , V 

i ^ 


T-) 


(A. 61) 


J-J 


# The perfect gas equation of state has been used to replace k^T in 
these equations from Ref. A. 8 with P/N. 
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k r\T 

^ ° + > h.J. 


h,2 ~ ^ 


T J . rr 

)>. , j,§9 ^>^2^ 


,£_\ \_-i_i-( ^ 

2 L U, 'T, V p . Pi 

i jf 1 


(A.62) 


%,3 


t^+Yh.j. , 

S§ 3 /_j ^ ^’^3 


^ ,Y "i “iJ ^ "j 


(A. 63) 


i j?^i 

where the components of the mass flux vector used in the above express roi 
are defined in the discussion of the species continuity equation. 

To calculate the components of the radiative flux vector 
where §. is an orthogonal coordinate, let us integrate Eq. A. 30. 


- 




.dr 






(A. 64) 


Note that V • q^isasca 


lar independent of coordinate system. The 


flux cornponents rnay be written. 


(A. 65) 


or by substituting from Eq. A. 30 


*R. 


f;(r ) ® / n^tT _ V 

^ = r 1 I ^4nB^ - j I^(r)dQydvh^dCj^ 


§(ro) 0 


(A. 66) 
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‘*R,2 

‘^R,3 



§(r 


O' 



^( r ^)0 


|.§(ri) 



_ \ 

I^(r)dQ jdvh2d§2 

0 

4rr _ \ 

I^(r)dQydvh^d§2 

0 



(A. 67) 


(A. 68) 


In addition to the general conservation equations an equation 
state is needed to specify the relationship between pressure and 
tenq^erature. A reasonable approxiaiation for the thermal behavior of 
a gaseous mixture is the ideal gas equation of state. 


P = p RT 

Another form of the ideal gas equation of state is 


(A. 69) 


P = Nk T 
c 


where 



N 


(A. 70) 


This last expression has been used previously to state Eq 


s A. 59 and 


A. 60. 



CONSERVATION EQUATIONS IN BODY ORIENTED COORDINATES 

In order to describe the flow over blunt bodies moving at hyper- 
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sonic velociLies, it is fomiu convenient to solve the conservation 
equations in orthogonal body oriented coordinate systems. The type 
of body under consideration, i.e. three-dimensional, axisymmetr ic or 
two dimensional, thus determine the stretching functions, h^^, h^, h^, 
discussed in the previous section. The class of bodies considered in 
this development are axisymmetric or two-dimensional and have the 
following stretching functions, see Tab. A. 2: 


= X, h^ = 1 + ny 


§2 ^ ^ 


I 

^2 


r 


3 


9’ 


h3 = r 


AXISYMMETRIC 


§1 = X, = 1 + Ky 


^2 


= y. 



1 


^3 


z. 



1 


TWO-DIMENSIONAL 


(A. 71) 


(A. 72) 


where k is the local body curvature and r is defined in Fig. A. 3. 
Using Fig. A. 3 the following relationship may be found 


r = r + y sin 0 (A. 73) 

w 

dr = sin 0 dy + k cos 6 dx (A. 74) 


where 


K = 1 + ny 


(A. 75) 
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(dC)^ = K^(dx)^ + (dy)"- + (r“dz) 


A, .2 


(A. 76) 


For the axisytnmetric case z has been substituted for cp. Thus by 
noting that the superscript A can be either 1 or 0 both the respective 


axisymme 


trie and two dimensional cases can be represented by one set 


of equations. 

Substituting the stretching functions A. 71 and A. 72 and relation- 
ships A. 73 and A. 75 into the general conservation equations for a 
tnulticomponent continuum gas in general orthogonal coordinates given 
in the second section yields the following equations. 


Continuity : 


9(pur 


Ur py)_ ^ 0 

dy 


Species Continuity : 

biv%C.u) hd rVi^) ' 

dx ^ dy dx 


- SCk r J ) _ ^ (A. 78) 

Sy ^ 

where J. and J . are the mass flux components of species i m the x 
L ,x 3. , y 

and y direction respectively. The mass flux vector is the sum of two 
vectors neglecting force diffusion and pressure diffusion. 

_ (A 79 

= h ^ 


The components are 
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concentration diffusion: 


(D) 

i,x 


2 


n 


Jt 

p (RT 


n 


I 



h W, /T,-V,Y 


1 ^\l 


Sx 


J 


K 






p IRT 


n 


\ — ' 

i 


r 

M.M.D. .1 Y. 
1 j ijL 1 




(A. 80) 


(A.81) 


thermal diffusion: 


(T) 

T 

®i S^n T 

(A. 82) 

i,x 

” Sx 



H 


(T) 

T S£n T 

(A. 83) 

i.y ~ 

” i By 



The two momentura equations can be expressed in the following 


manner . 

X - momentum 


pr\ + pv;r“v tr, “ 


~ A ^ 

Sy 

A . ^ k 


+ r 


A9P 


Sx 


5x 


Sy 


-V ** 

A . Sr _ 
xy zz 3x 


y - momentum 
A dv 


~ A Bv „ 2 

pr u -^ + p,,r v-^ - pnr u 


+ H r 


A ap 


Sy 


Sx 


Sy 


A ^ ~ - n 

+ V" aT ■ 


(A. 84) 


(A. 85) 
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where the components of the stress tensor are 


^XX ^ 


^ LSx 

K 


X rs(/ul + + 2p 1^ 

V " ^ L Sx ^ Sy -1 5y 

yX 

vX 


T = T* = M- 

xy yx 


[■J: ^ ^ ul 

L~ Sx Sy ~ -I 

w X 


The above stress components are also used in the energy 
energy: 

S(rV ) o(f^r%p ) 

A 5H ~ A BH _ P _ — - ■ ’ ^ 

P’^ ” ^ Sy " Sx Sy 

^ A j 


B ~ A 


^ .2— I Kr uT + X r vt 
By L xy 


yyj 


The components of the diffusional heat flux vector are 


^D,x 


- + ) h.J. 

~ Bx Z-i 1 ijX 

^ i 

T 

_ ^ N. D. 

“li 


J . 

...liiS 


J- \ 

P.. 1 


(A. 86) 

(A. 87) 

(A. 88) 
(A. 89) 

equation. 

(A. 90) 


(A.91) 
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- + 
^ 3y ^ 


JL. V V ^ 

" / ^ i "i “^ij ^ pj 'i ^ 

3 - 2T-^ 

From Eq. A. 66 and Eq. A. 67 the components of the radiative flux 
vector are: 

n = r ^ f O' (4ttB - I (r)df2)dvxdx 
x(ro) 

y(rj.“ p4rr _ 

1 //. .1 T ^ I \ 1 \r 


(A.92) 


q = ^ L Ck' (4 ttB - I (r)dD)dvdy 

%.y '' ''•'o'' 


(A. 93) 


(A, 94) 


The statement of these vector components completes the set of 
conservation equations expressed in body oriented orthogonal 
coordinates. By the use of the stretching functions listed in Tab. 

A. 2, the conservation equations can be written in the coordinate 
system desired by following the method used for the case under consid- 
eration in this section. Furthermore, it is noted that the conservation 
equations were obtained using orthogonal versions of the tensor terms. 
Consequently, if a nonorthogonal transformation is desired the equations 
may be used without fear of neglecting terms. 

Subsequent transformation of independent variables using Dorodnitsyn, 
Von Mises, Lees or one of many other transformations may be made in 
order to simplify the form of the conservation equations. The 
selection and use of these transformations will not be discussed here. 

The reader is referred to Dorrance, Ref. A. 9, and Hansen, Ref. A. 10, 
for suitable discussion and listing of similarity transformations. 
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APPENDIX B 


RADIATIVE TRANSPORT EQUATIONS 

The purpose of this appendix is to provide details of the 
radiative properties and transport equations used in this work which 
were developed by Wilson (Ref. B.l). Accordingly, the following 
development provides a synopsis of the radiative transport section 
vjf Ref. B.l. The synopsis is complete in itself and provides definition 
of symbols and terms not found in the main text. 

BASIC TRANSPORT EQUATIONS 

In Chapter 2 the radiative transport equation was solved in 
physical space for a one-dimensional slab to yield an equation for 
intensity. This equation, Eq, 2.49, and the equations for the flux 
Eq. 2.50, and flux divergence, Eq. 2.51, are the starting place for 
this development. 

In order to calculate intensities the absorption coefficient is 
separated into a line and continuum contribution. 

C , I" (B.l) 

O' = a + a ^ 

V V V 

Correspondingly, the flux normal to the body is expressed as a sum 
of the continuum only process and the line process corrected for 
continuum attenuation. 

■Jr ' \ * 

* The y subscript on the flux term which was used in the main text 
has been dropped in Appendix B for convenience. 
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The continuum flux is then expressed as 
z.” r ^E. (0,y) 




(T 


1 u 


B dE - V 
V V 


^ B dE Vdv 
0 ^ 


(B.3) 


= I - exp[- j J 


1 

where the emissive function, E^, is 

r-y 
y 

The line flux is expressed as 

;.w,(0,y) 

■*R 


(B.4) 


all lines 0 ® 


where the frequency integration has been carried out analytically and 
incorporated into an equivalent width variable, w. (y^y) , defined as 


W^(y' ,y) = (y.y')^ = 


1- e ^ 


f 


a dy 

/ V 




Jdv (B.6) 


In this definition it has been assumed that the continuum absorption 
coefficient and Planck function are frequency independent over the 

interval Av. 

The flux divergence is expressed as a sum of four terms. 


= Q = 

hy 

These four terms are defined as; 

(1) the energy emitted and absorbed by the continuum; 


(B.7) 


* Throughout this appendix the parenthesis sj'mbol ( ) is principly 
used to denote functional dependencies. 
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■>c,c _ 


c E (0,y) E (y,6) 

2n M B^(y^)dE (y,y') + B (y')dE (y ,y) 

On V LO A V V V ^ V 


0 '-"0 


- 2 B 


^(y) jdv 


(B.8) 


(2) the energy emitted by the continuum and absorbed by the lines; 

Q ’ = 2ti ^ J B^(y )[_S^(y) - A^(y ,y)jdE^(y ,y) 

all lines 0 

X 


r^E (y,6) 

B^(y')^S^(y) - A^(y,y')jdE^(y,y') 


(B. 9 ) 


( 3 ) the energy emitted by lines and absorbed by the continuum; 


Q^><= - 2„ > 


\ ( y ) J B^(y')dw^(y ,y) 

all lines 0 

i 


,^w.(y,6) 

+ J B^(y )dvi(y,y^) 

0 


( 4 ) the energy absorbed and emitted by lines: 

L L r- (^i<o.y) , - 7! a,, (y)<Jy 

^ I B (y )e ■^''y i dA (y,y^) 

Jq 1 1 


Q * = 2 n 2 . 

all lines 
i 


+j' 


1 |.y' c 

- 7 „ O' (y)dy 

j B (y )e ^ ^ 1 dA (y,y') 

> A J * i 


(B.IO) 


- 2 B.(y)S.(y) 


(B.ll) 


In Eqs, (B. 9 ) and (B.ll) the frequency integration is performed in 
terms of an absorption equivalent width variable A^(y^,y) defined as 



(B.12) 


A . ,y) 


A.[-(y,y^)] 




The line strength variable, S., is defined in section on transport by 
a single line. 

transport y.Y the COy^INUUM PROCESS 

Since the continuum absorption coefficient varies slowly as a 
tunction of frequency, except at photoionization edges, a monochromatic 
evaluation of the flux and flux divergence at selected frequencies 
followed by numerical quadrature over frequencies is realistic. 

Four atomic species are considered, H, C, 0, N. For these species 
and shock layer conditions of interest, the shock layer is optically 
thill below the first major photoionization edge. Furthermore, at 
frequencies above the first photoionization edge the absorption 
coefficient is essentially constant. Consequently the frequency 
dependent absorption coefficient can be represented by a series of 
grey absorption coefficients. For the optically thin groups the grey 
absorption coefficient is the Planck mean coefficient. 

Df B dv 

V V 

(B.13) 

B dv 



where 


a 

V 



j 


I 

j 


N.cr.ri - 

J JL 


-hvATI 

e 


(B.14) 
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and 


e)j(v,T) for j 


species H, C, 0, N 


The absorption cross section has the same form for the atomic species 

2 

considered (all cross sections are given in cm /particle). 


0 hv s hv^^^ 


a. = 7.26 X 10 r. kT e 
J 1 


16 „ - [ej 


- hv - A]/kT 

[hv]- 


hv > hv„ 


16 - 1 ^"- 
a. = 7.26 X lO" r. kT e ^ 


hv^J]/kT g. 


[hv]' 


v;here 


(B.15) 


(B.16) 


Fj “ species statistical weight factor 
= species ionization energy 
A = photoionization edge shift 
hv^ = species merged energy level limit 
g. = species nonhydrogenic correction factor 


Reference B,1 did not include molecular absorption or emission. How- 
ever, continuum cross section for the molecules N 2 , 0^, C 2 > C^, and CO 
were more recently included by Wilson in his Fortran subroutine TRANS 
which was obtained from Wilson and forms the basis of the radiative 
calculation method of the present work (Ref. B.2). In the present 
work the cross sections for and C^H were added. The addition of 
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molecules to the radiative calculators dictated several changes in the 
contjnurm calculation procedure described in Ref. B.l. The following 
development includes these changes. 

The continuum transport is calculated using a 12-band model (7- 
bands were used in Ref. B.l). In each group the absorption 
coefficient 


^ " Z. ^ ^3^ ^2^ ^2^ ^2’ ^2^^2^ (^-17) 

j 

is calculated. The equations for the band averaged absorption cross 
sections for each band group is listed below. 

Group 1: 0 ^ hv ^ 5.0 ev 

Atomic cross sections 

In one band the partial Planck mean is used in which a ^ for each 
element has the form 


CT . 

J 


5.04 X 1o\t r .K -e./'kl 


j O J 

£ J 


A/kT, jr. 

e hVj 





r 71 ■[‘"’r - 

+ kT je^ + 2pIkT] - kT e -* \0j + Pj 


(hv 


- hv^J) } 


B 

2"^ 


-{hvg- hv^^}{ 23 .[hvg- hv^JjkT + 23 .[kT] 2 } 

K i- e 


(B.18) 


^ , -16 2 2 
Kq = 7.26 X 10 cm ev 

hv„ = 5.0 ev 

D 


where 
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lA^n - hvw/kTp/^^B\^ 

w+ X 10 [kT] - e j 


^ & - K^) - ^]} 


For each atomic species the required parameters are; 



hv^'^(ev) 

Cj (ev) 

CD 

C.I. 


H: 

2.40 

13.56 

1.00 

0.0 

C: 

3.78 

11.26 

0.30 

0.0488 

N: 

4.22 

14.54 

0.24 

0.0426 

0: 

4.22 

13.51 

0.24 

0.0426 

Molecular 

cross sections 



O 

b 

2 2 

^2 S 

^CO ^ 

= 0.0 
2 " 

Group 2: 

5.0 5 hv ^ 

6.0 ev 



Atomic cross sections 




a . = 
J 

Eq. (B.16) 

j = N, 0, 

H, C 


Molecular 

cross sections 




= 3. X 10“^® 

-18 

+8. X 10 

- .5/kT 
e 



= 4. X lO"’-^ 




Q 

O 


^CO " 

= 0.0 
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Group 3: 6.0shv^ 7.0 ev 

A p. ,• 1 - r r> cti ons 


= Eq. (B.16) 


j = N, 0, C, H 


Molecular cross section^ 


CT = 1. X 10 
^2 


18 


^ , -18 - .7/kT 

°C0 - 3- 10 

"<>2 “ "“2 ■ ‘’•*2 ' "S ' ° 

Group 4: 7.0 ^ hv ^ 8.0 ev 

Atomic cross sections 


= Eq. (B.16) j = N, 0, H 

, -17 - 4.18/kT 

= Eq. (B.16) + 5. X 10 e /T^ 


Molecular cross sections 

, -17 - .5/kT 
a = 1.9 X 10 e 
*^C0 


CT^ = 6.0 X 10 

“2 


“Ic H " “ 


-19 


■18 


C2 N2 H2 C3 


= 0.0 
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Group 5; 8.0 ^ hv ^ 9.0 ev 


!3 sections 


(j = Eq. (B.16) j = N, 0, H 

j 


CT = Eq. (B.16) + 5. X 10 e 

c 


17 - 4.18AT 


-17 - 2.68/kT 
+ 2.2 X 10 e 


Molecular cross sections 




Q = 2.0 X 10 

°2 


Q = 8.5 X 10 


“"2 ■ ■ "S ' "S 


= 0.0 


Group 6: 9.0 5^ hv ^ 10.0 ev 


Atomic cross sections 


= Eq. (B.16) j = N, 0, H 


, -17 - 4.18/kT, 
a = Eq. (B.16) + 5. X 10 e 

c 

-17 - 2.68/kT. 

+ 2.2 X 10 e /5^ 


Molecular cross sections 


°co = “ 


Q - 1.0 X 10 
^2 
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, - CTj^ - cTh cr^ <^C„H 

'2 2 2 3 2 

! 0.0 ^ hv ;S 10.8 ev 


= 0.0 


Atomic cross sections 




a 


N 


= Eq. (B.16) j = 0, N 

, - 18 ^ - 10 . 2 /kT 

= 3.2 X 10 T e 




Or 


To ^ 10-17 - 1 . 26 /kT . 2 2 X lo"^^ e“ 2 - 75 /kT 

^ 8.5 X 10 e + 2.2 X lU e 


+ 5.0 X 10 


-17 - 4 . 18 /kT"',,^ 


Molecular cross sections 


_ ^ A 1 A 


19 


0 . 


L ^2 2 


= 0.0 


Group 8 : 10.8 ^ hv ^ 11.1 ev 


A t ora 1 c c ros s sections 


a = Eq. (B.16) j = 0, N 

J 

. 10-17 - 3 . 5 /kT 

= 5.16 X 10 e /Z^ 


'N 


'’c ■ L 


= l 8 . 5 xl 0 -^'e- + 2.2 X 10 - 1 ' 


CO 10-17 - 4.i8/kr^.,„ 

+ 5.0 X 10 e J'^C 


Molecular cross sections 




= 0.0 



349 


Group 9 : 11*1 ^ hv ^ IZ^O ev 


-sections 




= Eq. (B. 16 ) j = 0 , H 


N 


, , -17 - 3 . 5 /kT. 

5. 16 X 10 e /T^ 


-[B. 


5 X 10"^\“ + 2.2 X lO'^^e" ^ 


, -17 - 4 . 18 /kT-],_ 
+ 5.0 X 10 e 1 / 2 ^ 


Molecular cross sections 


CTj^ = 1.0 X 10 


18 


^2 ^2 ^3 

Group 10 : 12.0 ^ hv ^ 13.4 ev 


CO 


u 

V* 


= 0.0 


Atomic cross sections 


fj = Eq. (B. 16 ) j = 0 , H 
J 


a = [6.4 X 10‘^^ e" + 5.16 x 10‘^^ e‘ 

N L 

CT^ = [ 9.9 X lo"^^ + 8.5 X lo"^^ e" 

+ 2.2 X lO'^^ e" 2 - 75 /kT ^ ^ ^q -17 ^ 


Molecular cross sections 

,-17 


CT = 1.0 X 10 
2 


a„ = 2.7 X 10 

n.. 


-17 




. 75 /kT 







350 


Group 11: 13.4 ^ hv £ 14.3 ev 


At 07 T)ic cross s ections 

ajj = 1.18 X 


CTjj = 3.6 X 10‘^^/E(. 


" [ 


, , , -17 - 2.3/kT , . If,- 17 - 3.5/kT“j. 

6.4 X 10 e + 5.15 x 10 e 


_ ; „ Q ,f,-17 ^ . ,„-17 - 1.26/kT 

= I 9.9 X 10 + 8.5 X 10 e 

. _ , -17 - 2.75/kT ^ _ ,f,-17 - 4.18/kT‘;. 

+ 2.2 X 10 e + 5.0 x 10 e j/I^ 


Molecular cross sections 


a„ 


= 1. X 10 


-17 




= 2.7 X 10 


-17 


’0^ ■ ■ '’C3 “ "CO ■ "c^H ■ 


Group 12: 14.3 ^ hv ^ 20.0 


Atomic cross sections 


"h " "o “ 


(j^ = j^l.O X 10 + 6.4 X 10 e + 5.16 x 10 e J'^N 


.-18 


-17 - 2.3/kT 


-17 -3.5/kT‘ 


To Q 10-17 ^ , 10-17 - 1.26/kT 

(7^ = I 9.9 X 10 + 8.5 X 10 e 


^ - , -17 - 2.5/kT , . f, ^ 10-17 - 4.18/kT 

+2.2x10 e +5. Ox 10 e 
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Molecular cross sections 


- 1.0 X lO"^® 



In the above expressions E. 
species , 



0.0 


is the partition function of the j~th 


The flux and flux divergence equations in terms of the grey 
absorption coefficients for each group k are: 


\,k = ^ J - J Bj^(y')dEj^Cy,y') 

0 0 

CC . 

-\(y,6)_ _ , 

\(y')dEj^(y,y') - 2 B^j 

0 

where 


(B.20) 


(B.21) 


Bj. - j B (T)dv (B.22) 

and Ej^ is determined using Eq. B.4 and the grey absorption coefficient 

C 

“k • 

TRANSPORT BY A SINGLE LI>?E 

Consider the transport by system of nonoverlapping, electron- 
impact broaden lines. An isolated line has a Lorentzian shape 
characterized by a strength S and (half) half-width y, neglecting 
line shifts, 
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L Sy . 1 1 

V Tl L , N + 

(v - Vq) + Y 


(B.23) 


W i i X. C 


s 



O' dv ~ 
V 


ne 

me 


N f / 
n nn 



-hv^/kT- 


(B.24) 


and where N is the lower state number density and f /, the f-number. 
n ““ 

The f-number represents the transition probability strength and is the 
number appropriate for a single line, or multiplet, or whatever 
collection of line transitions is represented by the lower state 
number density, N . 

The spatial dependence of the half-width y in the denominator of 
Eq. (B.23) precludes a closed- form evaluation of the frequency 
integration required by the equivalent v/idth variables W and A. For 
computational expediency, the actual half-width in the denominator 
of Eq. B.23 is replaced with a spatial average y . When the line 
center is optically thick in near constant temperature regions the 
effect of the half-width y in the denominator of Eq. B. 23 is negligible. 
Thus, the spatially averaged value for y is defined such that in the 
optically thin limit the correct flux is obtained. To this end, 
consider the spatial integral appearing in W(y,y') and A(y,y') of 
Eq s . B . 6 and B . 1 2 . 


Z 



L 

a 

V 


(y)dy 


Define a transport variable 


z 


/TT 


py 

^ s(y) y(y)dy 


(B.25) 


(B.26) 
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then using Eq. B.23 and approximating y(.y) in the denominator by y(z) , 
Eq. B,25 is rewritten as 


Z = z! 


Lr -i2 - 2 , , J 

[v - Vq] + y (z) 


(B.27) 


- 2 , 


where y determined. Note y(z) is constant over the 

interval y to but not over the entire shock layer. Using this half- 
^^idth approximation the following expressions for W(z) and A(z) are 
obtained. 


W(z) = 2n Y t e ^[l^Ct) - ][continuum attenuation] (B.28) 

A(z) = S(y)[l - e^’^IgCt)] (B.29) 

. —2 

where t = z/2 y 

and Iq and are modified Bessel functions of the first kind. For 
the optically thin limit t « 1 


W 



S(y)d; 


(B.30) 


requires 

y 

J ^ S(y)y(y)dy 

Y(z) = 

J ^ S(y)dy 


(B.31) 


For the flux divergence equivalent width calculation the appropriate 
value for y is the value at the local point at which Q(y) is being 
calculated, i.e. y = Y^y) • 
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TRANSPORT BY A COLLECTION OF ISOLATED LINES 

To calculate the flux and flux divergence, line groups are used 
where a ^'group*^ of lines is defined as a collection of adjacent lines 
within a spectral interval across which both the Planck function and 
continuum absorption coefficient can be approximated as being 
independent of frequency. The contribution from only neutral atom; 

H, 0, N and C; transitions are considered since the contribution of 
ionic lines are negligible. All lines are considered non- overlapping. 
The overlapping corrections are developed in the next section. 

Eq. B.5 yields the total flux from all lines in a group if 
W(y, y') is interperted as the sum of the equivalent width from each 
individual line, 

n 

W(y,y') = w (y,y') (B.32) 

< -- i lii 

m=l 

where n is the total number of lines in the group. The right hand side 

of Eq. B.32 is approximated with a single expression of the form 
n 

W^(y,y") = nW(S ,Y ) (B.33) 

nF= 1 

where S* and y* are line parameters averaged over all lines in a 
group. Assuming that all lines in a group are either optically thin 
or optically thick expressions for S* and y* are derived. 

For the optically thin limit Eq. B.33 reduces to 

I J / (B.34) 

y y 


m 
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which requires 


~ ; s 

n Z_i in 


(B.35) 


m 


For the optically thick limit the tvjo sides of Eq. B.33 reduce to 


T- r 'y ^ ^y ^ . 

Zj J / \ J , S*(y)y*(y)dy^ 


(B.36) 


m 


Due to the square root operator a further approximation is made. 


m 


V 1 


{ I [ J } ” J y- {I 


(B.37) 


m 


This approximation is exact if the spatial variation of the 

same for all lines or differs by a constant factor only. 

The half-width v is proportional to the electron number density 


v„(y) ' 


(B.38) 


where B is a normalized half-width and is constant. The line strength 
m 

is proportional to the lower state number density and the induced 
emission factor, 


s (y) = 1. X 10' ^^f N (y)[l - 
m ^ mm 


(B.39) 


When all lines within a group have a common lower state Eq. B.37 is 
exact. Using this as a basis Eqs. B.36 and B.37 yield 

2 


n 

m 


(B.40) 
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which requires 


I ' 


„ I 


m 


(B.35) 


m 


For the optically thick limit the two sides of Eq. B.33 reduce to 


y y 


(B.36) 


m 


Due to the square root operator a further approximation is made* 

{ I L J } ~ J { £ [\(y>v„(y)] } 


(B.37) 


m 


m 


This approximation is exact if the spatial variation of 
STima for all lines or differs by a constant factor only. 

The half-width is proportional to the electron number density 




(B,38) 


where B is a normalized half-width and is constant. The line strength 
m 

is proportional to the lower state number density and the induced 
emission factor, 


S (y) = 1. X N (y)[l - 

m mm 


(B.39) 


When all lines within a group have a common lower state Eq. B.37 is 
exact. Using this as a basis Eqs. B.36 and B.37 yield 

n S* 


V* “ 4- { Z 


(B.40) 


m 
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A similar treatment o£ the flu. equivalent width function A yields an 

• ^ ^ expression for y*- 

For each atomic specie Individually Eq. B. 37 is reasonable. Thus 
for each line group the equivalent width for all lines of a given atom 
are computed. The total equivalent width for all lines of the group 


becomes 

n 

V' 




(B.41) 


where H. is the lumped equivalent width for each atom and where n . , 

S* andy* refer to the effective line parameters for that atom. The 
parameters and y* are calculated from Eqs. B.35 and B.40 where 

the summation includes only those lines for a given atomic species. 
Using Eqs. B.33 and B.39 the line width and line strength can be 

expressed as 


-16 r:* Pi -hv/kT-i 

s * = 1 . X 10 N.f . U - e J 

j J J 


V.* = N B.‘ 

e J 


(B.A2) 

(B.43) 


The £.* and p.* terms are 

J d , 


n . 
-J 

) 

• 

f g e 
mm 

m 



n. 


4 L 

J J ,,, 


-j rfn&£sL 


■ e ykT 
m 


The line spectrum for H, C, N. 0. atoms was collected into nine groups, 
Hithin each group four "effective lines" were considered, one for each 
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atomic species. The f-number and half-width y* are calculated 

by rown'tting Eqs. B.44 and B.45. 


J 1 


p.*=— I — y 

J J 


(B.A7) 


where all transitions energies are calculated with respect to a common 
lower state Z and where is the fractional population of state Z 

in species j. 


j 


-e^/kT 




(B.48) 


f J = Y f (B.49) 

Z A in 

m in X 

m in jJ 

The states Z considered for the H, C, N, and 0 species are listed in 
Tab • B . 1 . 

The data for the spectral lines considered is given in Tab. B.2. 
For each line group, its spectral location and interval is listed. 

For each element the number of lines n^ in the group and for each 
state Z of that element the parameters f^ and are listed. 

TRANSPORT BY OVERLAPPING LINES 

This section considers the correction to the group equivalent 
width which accounts for overlapping lines. 
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> f TP 1 


The correct e>.presslon for the equivalent »idth for a collection 
linea vithin a group is (on-itting the continuum attenunation 


factor for simplicity) 


group J 


- I - -K- I f 


(B.51) 


where the summation is over all lines and D is the spectral interval 
covered by the group. For the case of no lines overlapping in a group. 

V- (B.52) 

W = W = ) W 

group i m 

m 

Hovever. when lines do overlap, an analytical result for the frequency 
integration in Eq. B.51 is not available. To avoid prohibitive calcu- 
lation time required by a frequency Integration an empirical 
correlation for the line group equivalent oldth for overlapping lines 


is used. 


The spectral interval D defining a particular line group is a 
fixed interval. When the lines are optically thick and strongly 
overlapping within the interval, the line group equivalent width 
approaches the value of D. Thus a means of measuring the amount of 
overlapping within a group is to compare the isolated line value W^ 
with D. By comparison with exact calculations Wilson showed that the 
group equivalent width was correlated quite v^ell by 


1C- 

2 ^ - 1'. E ; 

= tan L 2 D J 
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For the flux divergence calculation, the parameter effected by 
is 

r(y,y') = S(y) - A(y,y') = 


ff *^(y)e "y 
V 


J^,aJ^(y)dy 

\T ^ r 


(B.54) 


which, for a group, is 


Tgro^p ■ J„{ A 

in 

For non- over lapping lines this parameter becomes 


group 


= ^1 = I 


Numerical results indicated that this parameter also could be 
correlated vjith the ratio W^/D by the expression 



The two correlations, Eqs. B.53 and B.57, were found sufficient to 
account for line overlapping in the transport process. 


local solutions for the flux D IVERGENCE 

The numerical evaluation of the line flux term, Eq. B.5, and the 
flux divergence terms, Eqs. B.8, B.9 and B.IO, present no particular 
problems. However, as noted by Wilson the numerical evaluation of 
of Eq. B.ll presents difficulties for the case of optically 
thick lines. To eliminate this problem Eq. B.ll is written as 
(omitting the continuum attenuation for convenience) 
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q'”'- . 2 n J 


Atoy-) ’ 


L,L,- 

■*■ ‘^local ^local 


(B.58) 


where 


= 2ni J 

local 


^ '*^^^^'\(y')dACy,y') - S(y)B(y)] 


(B.59) 


= 2n' 

local 


_ s(y)B(y)] 

^ -^0 


(B.60) 


Difficulties occur when, for a line 
AA” -♦ S(y). The solution is found by eva 
or B .60 by parts. For example, 


which is very optically thick. 


luating the integral in B.59 




dB 


qL,L,- ^ B(At')A(At ) - J dt 

local ^ 


where 


t = z/ 2 y^. Replacing dB/dt by constant 


(B.61) 


^ ^ BfAt'^ - B(01 
dt " At" 

Eq. B.61 becomes 


= _ B(At")r(A 

local 


t')r(At") + [B(At") - B(0)] A(At") 


(B. 62 ) 


(B.63) 


where 


r(At") = S(0) - A(At") = /2) 


(B.64) 
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A(At") = j A(t)dt = 


-At*‘/2 


[iQ(At m 


+ I^(At"/2)] 


(B.65) 


EquaLxmi provides 


the correct limiting 


form for the flux diver- 


gence for the case of very large optical depth, At »1 

equation for is used, 

local 


An analogous 



362 


table B.l (from Ref. B.l) 
fractional population data 



363 


1 . 

TABLE B.2 (from Ref. B.l) 


line group data 



Note; 
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table B.2 (Ref. B.1-) 
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